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Recent Investigation of Temperature 
Recovery and Heat Transmission on Cones 
and Cylinders in Axial Flow in the N.O.L. 

Aeroballistics Wind Tunnel 


G. R. EBER* 
U.S. Ordnance Laboratory 


ABSTRACT 


The Naval Ordnance Laboratory research program on the 
investigation of aerodynamic heating effects at supersonic veloc- 
ities is outlined briefly. Results of recent experimental investi- 
gations on temperature recovery and heat transmission along 
the surface of cones and cylindrical bodies in axial flow are pre- 
sented. The experiments were conducted in the N.O.L. 40 by 40 
em. Aeroballistics Tunnel at Mach Numbers between 1.5 
and 5.0, inclusive, and cover the regions of laminar, transition, 
and fully developed turbulent boundary-layer flow. The ex- 
perimental results are compared with analytical solutions for 
temperature recovery and heat transfer with laminar and tur- 
bulent boundary layer. The qualitative effect of the direction of 
the heat flow and the surface temperature on the limits of the 


transition region is discussed in detail. 


INTRODUCTION 


7 NEED FOR ACCURATE HEAT-TRANSFER DATA iS 
apparent if one considers that the surface of a 
supersonic vehicle flying in the atmosphere with 2,500 
m.p.h. or a Mach Number of about 3.4 will adopt an 
equilibrium temperature of about 540°C. above that of 
the ambient air. This rise in temperature is caused by 
conversion of kinetic energy of the flying body into heat 
by compression and friction in the boundary layer, 
which increases the temperature of the air in proportion 
to the square of the speed of the flying body. 

How much of the heat produced will be transferred to 
the surface depends greatly on the surface coefficient of 
This coefficient in turn depends on the 
Our knowledge 


heat transfer. 
flow conditions in the boundary layer. 
of the aerothermodynamic characteristics in supersonic 
Presented at the Hypersonic Aerodynamics Session, Nine- 
teenth Annual Meeting, I.A.S., New York, January 29—February 
1,1951. 
* Chief of Aeroballistic Division. 


flow is rather incomplete. It is, therefore, of great 


importance to investigate surface temperatures ann 
heat transfer in this range in order to provide the de- 
signer with the data that are needed so urgently. 
Considerable theoretical and experimental work has 
been expended in recent years on the solution of the 


supersonic heat-transfer problem. Despite many valu- 
able results that have been obtained, however, we are 
still far from a solution of this problem, especially in the 
region of turbulent-flow boundary layers. More theo- 
retical and experimental work must be done before we 
will be able to predict supersonic heat transfer reason- 
ably accurately. 

Accordingly, the Naval 
(N.O.L.) has started a comprehensive experimental and 


Ordnance Laboratory 
theoretical research program on supersonic heat trans- 
fer. The experimental program is being conducted in the 
White Oak Aeroballistic Research Facility; it consists 
of the following phases: (1) temperature recovery fac- 
tors, (2) heat-transfer correlations, and (3) bound- 
ary-layer transition. 

Although this program is far from completion, some 
significan® results have already been obtained and are 


presented here. 
EXPERIMENTAL FACILITIES AND INSTRUMENTATION 
The principal experimental facility is a 40 by 40 cm. 


The wind tunnel operates 
A range of Mach 


aeroballistics wind tunnel. 
with atmospheric supply pressure. 
Numbers between 1.2 and 6.5 can be covered in this 
tunnel by a set of fixed nozzles. (The tunnel, at present 
intermittent, is described in detail in reference 1.) 
Different models have been employed to produce 
laminar, transitional, and turbulent boundary-layer 
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20° DEGREE CONE MODELS 
Fic. 1. Three 20° cone models for recovery-factor and heat 
transfer measurements: (a) thin-walled cone, (b) thick-walled 
cone, and (c) cone with insulated frustra. 
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Fic. 2. Cone-cylinder model. 
flow. For laminar boundary-layer flow, small cones 


have been used. Three typical 20° cone-models are 


presented in Fig. 1. The first cone (a) is a thin-walled 
copper cone of 0.025-in. wall thickness with several 
thermocouples inserted along the skin of the cone for 
surface temperature measurements. The second cone 
(b) has a thick wall of copper, the heat capacity of 
is used for the determination of the 


which average 


heat-transfer coefficients. The third cone (c) is a cop- 
per cone subdivided into frustra. The latter are ther- 
mally insulated from each other to avoid heat flow in 
the axial direction and to allow the determination of 
local heat-transfer coefficients. Except for the sub- 
divided cone, cones with 10, 40, 60, and SO deg. total 
angle were used throughout the investigation. 

For transitional and turbulent boundary-layer flow, 
there were employed cone-cylinder models, as shown in 
Fig. 2. These models are composed of a 40° cone and 
a number of cylindrical sections 2 in. in diameter, up to 
2 in. longest. The cylinder can be built up to about 40 
in. total length. For insulated surface-temperature 
measurements, the sections are made of thermally in- 


sulating material; a surface temperature element as 
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shown in Fig. 2a can be inserted at any desired position 
of the cylinder. Four cu-const. thermocouples of low 
heat capacity are inserted equally spaced around the 
circumference of the measuring element. 

For heat-transfer measurements 2-in. sections of 
copper rings of '/s-in. wall thickness are used, as shown 
in Fig. 2b. These sections are mounted on a heat- 
insulating body. One to four thermocouples are in- 
serted in the copper rings; the heat capacity of the 
copper ring is used for the determination of the heat- 
transfer coefficients in the same way as done with the 


cones. 


TEMPERATURE RECOVERY FACTORS 


The surface of a body in supersonic flight will assume 
its highest temperature after it has been heated up to 
the temperature of the inner edge of the boundary layer, 
Then the transfer of heat has become zero. This maxi- 
mum temperature depends on the energy distribution 
in the compressible boundary layer and is usually ex- 
pressed in terms of the temperature recovery factor , 
This factor is defined as the ratio of the actual tempera- 
ture rise across the boundary layer to the adiabatic tem- 
perature rise. The recovery factor, in general, is a 
function of the similarity parameters, Reynolds Num- 
ber Re, Prandtl Number Pr, and Mach Number J, 
Thus, 
at — 5). fits, Py, M) 

Io- I, 

where 7; is the temperature of the insulated surface, 
7; is the ambient air temperature outside of the bound 
ary layer, and 7} is the stagnation temperature. 

According to the different energy distribution, we 
must expect different recovery factors on bodies with 
laminar, transitional, and turbulent boundary layer. 


{) Laminar Boundary-Layer Flow 


The results of the temperature measurements with 
the thin-walled cones are presented in Fig. 3. The local 
recovery factor determined with the local temperature 
at the inner edge of the boundary layer is plotted versus 
The 
local Mach Number was determined from the M.I.T. 


local Mach Number /, for the various cones. 


Tables* by disregarding the existence of the boundary 
layer itself, which was considered immaterial for the 
small cones. The measurements show that the tem- 
perature-recovery factor is 0.845 + 1 per cent independ 
ent of the Mach Number in the entire range and within 
the measurement accuracy in agreement with the the- 
oretical value r = Pr, provided that the Prandtl 
Number is evaluated with the physical properties of 
the air at equilibrium temperature, 7°. 

The obvious independence of Mach Number leads to 
the conclusion that the compressibility effect is the 
same at all Mach Numbers in the investigated range. 
It appears, therefore, that the above simple relation 
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TEMPERATURE 


can be applied with confidence in predicting insulated 
surface temperatures for laminar-flow boundary layer 
at least up to JJ = 5. It must be considered, however, 
that the recovery factors are determined at low temper- 
atures. It will be necessary to verify the results for 
actual flight temperatures. 

The Reynolds Numbers, encountered in the cone 
measurements with laminar boundary layer, presented 
The 


1.€., 


in Fig. 4, are ranging from 6,000 to 500,000. 
Reynolds Numbers are evaluated on 7, basis 
from the physical properties of the air at equilibrium 
temperature, the velocity v outside of the boundary 
layer, and the length of the cone along an edge. The 
investigations with the cones are reported in more detail 
in reference 3. 


B) Transitional and Turbulent Boundary Layer 


For higher Reynolds Numbers, surface temperatures 
were measured along the cylindrical part of the cone- 
cylinder model. The local flow conditions along the 
cylinder outside of a boundary layer were taken from 
computations the Ballistic 
Laboratories, Aberdeen, according to a method devel- 
Static pressure taken along the 


provided by Research 
oped by Clippinger.* 
cylinder was found to be within 2 per cent in agreement 
with the data given by Clippinger’s method. 

The temperature recovery factors determined from 
the surface Mach 
Numbers 2.87 and 4.25 are presented as a function of 
The experiments show 


temperature measurements for 


the cylinder length in Fig. 5. 
(a) that there is a separation of the data according to 
Mach Number for a cylinder length larger than 7 in. 
and (b) that for each Mach Number there are generally 
two distinct regions—one in which a rather large varia- 
tion of the recovery factor occurs and one in which 
the recovery factor is essentially constant. Schlieren 
photographs, taken with 0.5 microsec. flashlight (Fig. 
6), show that the boundary layer is fully turbulent in 
the region of constant recovery factors. Hence, the 
front portion of the cylinder can be considered as having 
transitional flow in the boundary layer and the rear 
portion as having turbulent flow. No means of de- 
termining precisely the beginning of the transition 
region on wind-tunnel models has yet been found. The 
value of the recovery factor | in. behind the leading 
This is about 4.5 per cent 
This high 


edge of the cylinder is 0.89. 
higher than the laminar value for the cones. 
value indicates that some turbulent motion exists al- 
As the air moves along the cylin 
more 


ready at this point. 
the turbulent 
pronounced, and the recovery factor increases until it 


der, motion becomes more and 


reaches a distinct maximum in the rear portion of the 
transition zone. The highest measured values are 0.96 
at Mach Number 2.87 and 0.98 at Mach Number 4.25. 
The recovery factor then drops to its turbulent values 
of 0.92 at M = 2.87 and 0.97 at M = 4.25. The data 


plotted against Reynolds Number are presented in 


RECOVERY 


AND HEAT TRANSMISSION 3 
The Reynolds Numbers are determined on 7, 
The character- 


Fig. 7. 
basis in the same way as defined above. 
istic length is the length of the model along the surface 
from the tip of the cone to the respective measuring 
point. The recovery factor results from the measure- 
ments on the cylinder may be summarized as follows: 

(a) In the transition region the recovery factor in 
creases to a value that is above that of the turbulent re- 
covery factor (an observation that was also made re- 
cently by Stalder’). 

(b) The recovery factor is essentially constant once 
turbulent flow in the boundary layer is established. 

(c) The turbulent recovery factor at JJ = 4.25 is 5 
per cent higher than at JJ = 2.87. 

(d) The Reynolds Number of transition to turbulent 
boundary-layer flow is lower for the higher Mach Num- 
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These experimental results differ from the usual con- 
cept of turbulent recovery factors deduced from ap- 
These analyses give the recovery 
i.€., 


proximate analyses. 
factor either as a function of the Prandtl Number 
r = Pr’ (as, for instance, Ackermann and Squire’)—or 
as a function of Prandtl Number and Reynolds Number 
(as, for instance, Seban*) but not as a function including 
the Mach Number. More data at low and high stagna- 
tion temperatures are needed to evaluate precisely the 
influence of the si:nilarity parameters Re, Pr, and M. 
However, it appears that the recovery factor is consid- 
erably higher than widely accepted from various analy- 
ses and the few experimental data from cones with 
artificially induced turbulence. The insulated surface 
temperature of supersonic vehicles with turbulent 
boundary layer must be expected to be near the stag- 
nation temperature. 


HEAT TRANSMISSION 


The next question is: How much of the heat pro- 
duced in a boundary layer will be transferred to the 
solid surface of a body at supersonic speed? The 
measured temperature recovery factors were used to 
determine the surface coefficient of heat transfer / and 
the nondimensional heat-transfer coefficient Nusselt 
Number Nu. The Nusselt Number is, in general, a 
function of the Reynolds Number, Prandtl Number, 
and Mach Number; thus, 


Nu = hl/k = f(Re, Pr, M) 


where / is the surface coefficient of heat transfer, / is the 
characteristic length, and & is the thermal conductivity 
of the air. However, in compressible flow the similarity 
parameters referred to free-stream conditions are not 
adequate for boundary-layer phenomena. The de- 
pendence of viscosity, thermal conductivity, and specific 
heat on the temperature suggests new variables that 
make the experimental investigations complicated and 
the comparison of the available data difficult. 





AT M=287 


40° cone-cylinder model at M = 2.87 in the N.O.L. 
40- by 40-cm. Aeroballistics Tunnel. 


Fic. 6. 
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(4) Laminar Boundary-Layer Flow 


Two sets of experiments have been made with laminar 
boundary-layer flow on cones: 
heat transfer over the entire cone was determined with 
thick-walled copper cones; in the other set, the local 


In one set, the average 


heat transfer for each section was determined with the 
subdivided 20° cone. 

Fig. 8 shows the heat-transfer data from the thick- 
walled copper cones in nondimensional form as Nusselt 
Number versus Reynolds Number. 
erties of the air are again evaluated at equilibrium 


The physical prop- 


temperature and the velocity at the outer edge of the 
The data are obtained for heat flow 
that is, into and out of the wall of 


boundary layer. 

in both directions 
the cone. The wall temperature, on the average, was 
either 25°C. below or above the equilibrium tempera- 
Local Mach Numbers at these measurements 
The data show 


ture. 
ranged from 0.88 to 4.65, inclusive. 
that the difference between heating and cooling is small. 
As in the case of the recovery factor, no influence of 
Mach Number was found. The data can be well rep- 
resented by the analytical solution for average heat 
transfer over cones, as given in reference 9, 


Nutay. = 0.767 Re®-®Pr'” 


In a similar way, the local values of the heat transfer, 
shown in Fig. 9, have been determined for the sub- 
divided cone in a range of local Mach Numbers from 
1.42 to 4.22. 
separation with respect to heat-flow direction. 


These data show even less evidence of a 
Again, 
the data can be well represented by the corresponding 


analytical solution for local heat transfer—viz., 
Nu = 0.575 Re®*®Pr” 
From these measurements we may conclude that 


application of the above heat-transfer correlations to 
predict supersonic heat transfer in laminar boundary- 
layer flow with reasonable accuracy is justified up to 
Mach Numbers of about 5. However, data at higher 
wall temperatures and higher supply air temperatures 
are still needed to verify the result in the high tempera- 
ture range. 


(B) Transitional and Turbulent Boundary Layer 


The 
boundary layer was determined with the cone-cylinder 


heat transfer for transitional and turbulent 


models. The recovery factor for each location along 
the cylinder observed under identical flow conditions 
was applied for computing the local heat-transfer co- 
the wall-temperature measurements. 


efficients from 


These heat-transfer coefficients were converted into 
Nusselt Numbers, in the same way as for the laminar 
boundary layer, by use of the equilibrium temperature 
as reference temperature for the thermal conductivity. 
The data, at J = 2.87, were taken with a model made 
up of 2-in. copper sections, thermally insulated from 


eachother. With this arrangement heat transfer occurs 
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TEMPERATURE RECOVERY 
over the full model length, the surface being at uniform 
temperature within +5°C. Again the heat transfer is 
measured for heat flow into and out of the cylinder wall 
below and above 
Fig. 10 shows 


at a wall temperature about 25°C. 
equilibrium temperature, respectively. 
the relation between Nusselt Number and Reynolds 
Number. Contrary to the result for cones with laminar 
boundary-layer flow, the change in the direction of the 
heat flow produces roughly 100 per cent change in the 
amount of heat transferred in the transitional, as well as 
in the turbulent, flow region. The question whether the 
separation of the data is due to the direction of heat 
flow or due to the difference in wall temperature cannot 
be answered from these experiments. Data in a larger 
range of wall temperatures are necessary to answer this 
question. 

The data from the starting section of the cylinder for 
the cold model seem to follow the equation for laminar 
flow over a flat plate as given by Crocco."© Taking into 
account the ratio of wall temperature and ambient tem- 
perature for the flow conditions prevailing on the model 
in the wind tunnel, the equation reads 


Nu = 0.295 Re®®Pr' 


The data for the heated wall indicate an appreciable 
amount of turbulent motion in the boundary layer. In 
the turbulent-flow region, the heat transfer for heat flow 
from the wall to the air follows closely the turbulent 
correlation of heat transfer as derived from the momen- 
tum-heat transfer analogy by Colburn'! for subsonic 


velocities—viz., 


Nu = 0.029 Re®8Pr'” 


The data compiled so far are not sufficient to make any 
final conclusions. The data, however, indicate that 
small changes of temperature distributions in the bound- 
ary layer near the wall are likely to cause large changes 
This is one reason for the nonuni- 
From the 


in heat transfer. 
formity of the available heat-transfer data. 
data that we have we must conclude that we cannot ex- 
pect a high accuracy in predicting turbulent heat trans- 
fer. The application of the well-known Colburn corre- 
lation for practical design purposes will give fairly ac- 
curate data up to Mach Numbers of about 5, as was 
found at free-flight tests by Fischer and Norris.'* In the 
higher Mach Number region, however, predictions are 
uncertain, as recently pointed out by Kaye,'* who at- 
tempted to determine temperature distributions in a 
wing flying at supersonic speeds. 


TRANSITION OF THE BOUNDARY LAYER 


Considerations of the stability of a laminar boundary 
layer are of considerable practical interest because of 
the large change in equilibrium temperature, as well as 
in heat transfer, associated with the change of bound- 
ary-layer flow. Boundary-layer theory shows that 
there is a principal difference between subsonic and 
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supersonic flow. Whereas in subsonic flow instability 
of a boundary layer will occur in all cases for sufficiently 
high Reynolds Numbers, in supersonic flow there exists 
a critical Reynolds Number above which a boundary 
layer will stay laminar. If the prediction as given in a 
theory by Lees and Lin" " is correct, radiative cooling 
might be sufficient in particular cases to keep a bound- 
ary layer laminar. 

In Fig. 11, the Reynolds Number of transition which 
can be observed in a wind-tunnel test is plotted as a 
function of the ratio of wall temperature, 7), to ambient 
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air temperature, 7), for ratios 7),/7 between 2.3 and 
3.0. The transition Reynolds Number increases by a 
factor of about 3 when the wall temperature is lowered 
over about 100°C. A wall temperature for which a 
boundary layer stays laminar over the entire length of 
the body has not yet been established. 


CONCLUSIONS 


The experimental investigations of temperature re- 
covery, heat transmission, and boundary-layer transi- 
tion in the N.O.L. Aeroballistics Wind Tunnel have 
shown the following: 

(1) Temperature-recovery factors in a laminar-flow 
boundary layer, measured on small cones with total 
angles between 10° and 80°, have been found to be inde- 
pendent of Reynolds Number and Mach Number in a 
range between JJ = 0.SS and 4.65 and can be repre- 
sented by the square root of the Prandtl Number for a 
Prandtl Number evaluated at wall conditions. 

(2) The transitional factor rises with 
Reynolds Number to a maximum that is higher than the 
turbulent recovery factor for the same Mach Number. 


recovery 


‘at Mach 


JANUARY, 1952 


SCIENCES 


(3) The turbulent recovery factor obtained from 
temperature measurements on cylinders at \/ 2.87 
and 4.25 is found to be independent of Reynolds Num. 
ber but 5 per cent higher for 17 = 4.25 than it is for 
M = 2.87. ° 

(4) The Nusselt-Reynolds relation for heat transfer 
with laminar boundary layers, determined from cones 
0.SS 


and 4.25, can be 


No influence of 


Numbers between 
represented by known correlations. 
Mach Number and little influence of the direction cf 
heat flow has been found. 

(5) Transitional and turbulent heat transfer for 4 
cone-cvlinder model at J = 2.87 shows separation of 
the data according to the direction of the heat flow, 
There are not sufficient data yet to approve or dis- 
approve existing approximate analyses. 

(6) All the experiniental data on temperature re 
covery and heat transfer are taken — 
+50°C. Verification of the data at higher tem- 
peratures will be indispensable tor practical applica- 


between 


and 


tions. 

(7) The Reynolds Number for transition to turbulent 
boundary-layer flow becomes lower for higher Mach 
Numbers. For a Mach Number, _ this 


Reynclds Number increases as the model wall becomes 


constant 


cooler. 
REFERENCES 
! Lightfoot, J. R., The Naval Ordnance Laboratory Aeroballistic 
Research Facility, Naval Ordnance Laboratory Report (N.O.L.R 
No. 1079, 1950. 


2 Tables of Supersonic 
Institute of Technology, Technical Report No. 1, 1947. 


Flow Around Cones, Massachusetts 


3’ Eber, G. R., Determination of Temperature Recovery Factors 
on Cones in the NOL 40 X 40 cm Supersonic Wind Tunnel No. 2, 
Naval Ordnance Laboratory Memorandum (N.O.L.M.) 10107, 
1950. 

‘ Personal communication with R. F. Clippinger and J. H 
Giese of the Ballistic Research Laboratories, Aberdeen, Md., 
1950. 

5 Stalder, J. R., Rubesin, M. W., and Tendeland, T., A 
Determination of the Transitional-, and Turbulent 
Boundary-Layer Temperature-Recovery Factors on a Flat Plate in 
Supersonic Flow, N.A.C.A. T.N. No. 2077, 1950. 

6 Ackermann, G., Plattenthermometer in Stroemung mit grosser 
Geschwindigkeit und turbulenter Grenzschicht, Forschung auf dem 


Laminar-, 


Gebiete des Ingenieurwesens, Band 13, p. 226, 1942. 

7 Squire, H. B., Heat Transfer 
British Air Ministry R. & M. No. 1986, 1942. 

*Seban, R. A., Analysis for the Heat Transfer to Turbulent 
Boundary Layers in High Veiocity Flow, Ph.D. Thesis, University 
of California, Berkeley, Calif., 1948. 

® Johnson, H. A., and Rubesin, M. W., Aerodynamic Heating 
and Convective Heat Summary of Literature Survey, 
Trans. A.S.M.E., Vol. 71, p. 447, 1949. 

Crocco, L., Laminar Boundary Layer in Gases, Translation, 
North American Aviation, Inc., Report CF-1038, 1948 

11 Colburn, A. P., A Method of Correlating Forced Convection 
Heat Transfer Data and Comparison with Fluid Friction, Trans 
A.I.C.E., Vol. 29, p. 174, 1933 

(Continued on page 14) 


Calculations for Airfoils, 


Transfer 





Am 
of win 
that li 
by sin 
tions ‘ 
wings 
lar to 
rect W 


R« 


gram 
nolog 
Dep: 

+} 





‘| from 


2.87 
» Num- 


t is for 


ransfer 

cones 
‘an be 
nee of 


ion cf 


for a 
ion of 
flow. 
r dis- 


re re- 
— 50 
tem- 


Dlica- 
ulent 
ach 


this 


ymes 


llistic 
L.R 


isetts 


ictors 


v0.2 


107, 


H 
Md., 


sser 


lem 
nls, 


lent 
sity 


ing 


mn, 


on 





The Lift, Rolling Moment, and 
Pitching Moment on Wings in 
Nonuniform Supersonic Flow 


HENRY L. ALDENT ann LEON H. SCHINDEL? 
Massachusetts Institute of Technology | 


SUMMARY 


A method is developed for computing aerodynamic coefficients 
of wings in arbitrary nonuniform supersonic flows. It is shown 
that lift, pitching moments, and rolling moments can be obtained 
by single integrations without calculation of pressure distribu 
tions as an intermediate stage. Application of the method to 
wings with supersonic leading and trailing edges is formally simi 
lar to conventional strip methods but gives results that are cor 
rect within the limits of three-dimensional linearized theory. 


NOMENCLATURE 


a = tangent of sweepback angle of leading edge 


b = wing span 
( = wing chord at root 
d = tangent of sweepforward angle of trailing edge 
é = intersection of Mach line and leading edge 
= rolling moment 
mi = pitching moment 
] = dynamic pressure (free stream 
t = parameter locating spanwise position of upwash pulse 
u = velocity perturbation in x-direction 
= velocity perturbation in y-direction 
w = upwash velocity perturbation in z-direction 
Wp = upwash due to pulse 
x = chordwise coordinate on wing 
y = spanwise coordinate on wing 
= vertical coordinate 
A.R. = aspect ratio 
C1 = lift coefficient 
a = pressure coeflicient 
F(t) = influence function for lift 
L = lift 
WV = Mach Number (free stream 
Pit pitching moment influence function 
R(t rolling moment influence function 
= free-stream velocity 
Gy local stream angle at wing tip 
6 = VM-1 
n = spanwise coordinate on wing 
é chordwise coordinate on wing 
¢ = perturbation velocity potential 


INTRODUCTION 


COEFFICIENTS OF LIFT AND PITCHING MOMENT 


HI 
I of supersonic aircraft are known to be affected ap 
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preciably by downwash flows over tails and aft wings. 
In cases where there are stringent requirements on con- 
trol and stability of the air frame, these effects must be 
included in any accurate performance estimate of wing 
tail combinations. Much work, both theoretical and 
experimental, has been carried out on various phases of 
the problem, including studies of pressure interference 
between wings and bodies and between vertical and 
horizontal surfaces of cruciforms, and investigations 
into the structure of flow fields behind lifting surfaces.! 
If a flow field is known, the problem of downwash effects 
is resolved into the calculation of aerodynamic coef- 
ficients of tail surfaces that operate within the nonuni- 
form supersonic flow. In principle, this calculation 
can be carried out by integrating fundamental solutions 
of the wave equation’ to obtain pressure distributions 
and, hence, the aerodynamic coefficients; in practice, 
the labor required is enormous. This report deals with 
the problem by developing a simple method for the 
direct computation of lift, pitching moment, and roll- 
ing moment without recourse to the intermediate stage 
of pressure distributions. The general approach em 
ploys the idea of lifting strips, as presented by Lager- 
strom and Van Dyke?* for the case of straight trailing 


edges. 
Basic EQUATIONS 


The theoretical derivation is subject to the restric- 


tions of nonviscous linearized theory and is based on the 
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2. Areas of integration for determining velocity potential. 











(a) (b) 


Fic. 3. Upwash pulse on trapezoidal wing. 
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Fic. 4. Limits of integration for pressure coefficient at points 


P\(x,y) and P.(x,y). 
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Areas of integration for determining F(t) at 4; and fe. 


standard form of the wave equation for the velocity 
potential of small disturbances, 
,7¢e We de 


— = ( (1) 


6 ox* oF 3s" 


where 0¢/0x = u, 0¢/Oy = v, and 0¢g/0z = w. The 
flow is represented as a stream of velocity U’ in the y- 
direction (see Fig. 1) with perturbations v and w. The 
absence of a perturbation in the x-direction (~) implies 
a constant-pressure field and no variation of v or w in 
the flow direction, a situation that physically corre. 
sponds to a downwash field several chord lengths be- 
hind a wing or at least far enough downstream for the 
flow pattern to have stabilized. Associated with the 
field is a velocity potential that satisfies Eq. (1). In 
this field is placed a thin airfoil that will disturb the 
field in a manner determined mathematically by a 
second potential also satisfying Eq. (1). It is this 
second potential that will describe the pressures and 
aerodynamic behavior of the airfoil and which will be 
termed the wing potential ¢ for the purpose of theoreti- 
cal development. 

The airfoil surfaces under consideration are taken as 
flat planar surfaces with supersonic leading and trailing 
edges, as well as tip edges that are aligned with the 
stream. Within these limitations the wing may have 
arbitrary angles of leading-edge sweep, trailing-edge 
sweep, and aspect ratio. 

The boundary conditions are thus set up for flat 
plates, since the increments due to thickness are addi- 
tive and can be treated independently. In the plane 
of the wing, the flow field has an upward velocity w(y 
and a lateral velocity v(y) which are constant in the 
chordwise direction; to ensure that flow over the wing 
is tangent to the surface, it is sufficient to specify that 
the vertical velocity of the wing potential cancels the 
vertical stream velocity or 


w(y) + (O¢g/dz) = 0 (2 


on the wing in the plane z = 0. Off the wing in the 
z = 0 plane, the pressure will be zero.2* Since C, = 
—2u/U = 0, 

u = O¢/dx = 0 (3) 


off the wing in the plane z = 0. Provided there are 
no subsonic edges, Eq. (1), with the boundary condi- 
tions expressed in Eqs. (2) and (3), determines a wing 
potential that will vanish at an infinite distance from 
the plane of the wing. For the present purpose, the 
argument is carried through for wings with supersonic 
leading and trailing edges, as indicated in the typical 
plan form in Fig. 1. The problem, as stated here, is 
completely equivalent to a warped wing in a uniform 
flow if small chordwise strips are inclined to the stream 
at various angles equal to w(y)/U. 

By superimposing fundamental source solutions to 
Eq. (1), several authors have been able to demonstrate 
that full solutions to the above problem can be expressed 
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LIFT, ROLLING 
formally as integrals. The lift, for example, can be 
written as a triple integral with variable limits. The 
wing potential, a double integral over the fore-cones “a 
Mach lines, is given as? 


=. [f g(é, n)dédn 4 
(x, ¥) = ) 
4 tJ JA V(x — 8)? — BXy — 0)? 


where g(£, ) = upwash distribution function. Sample 


MOMENT, 
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Fic.7. Linear symmetric upwash distribution w(t) and influence 


function F(t) for rectangular wing of effective aspect ratio 4. 
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AND PITCHING MOMENT 9 
areas of integration which would give the potential on 
the trailing edge are indicated in Fig. 2, showing that, 
for points (x, y) not influenced by tip Mach lines, the 
region A is the Mach fore-cone, whereas in tip zones 
the region is modified as shown. To evaluate the lift 
in a field of arbitrary vertical velocities by pursuing 
this method is generally impractical. A different ap- 


proach is desirable. 
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Fic. 9. Influence functions for rectangular wing of effective 
aspect ratio 4(8 A.R. = 4). (a) F(t) = lift influence function. 
(b) R(t) = rolling moment influence function (about the wing 
centerline). (c) P(t) = pitching moment influence function 
(about the mid-chord ). 
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It is assumed that the upwash can be divided into 
elements, as indicated in Fig. 3, with the objective of 
computing the increment of lift from a representative 
The total lift results from the proper addi- 
Specifically, an upwash element or 


element. 

tion of increments. 

yulse of unit strength, w,édn = 1, is placed at spanwise 
pen 

The 


pressure disturbance will spread over the Mach cone 


station, » = ¢, while everywhere else w = 0. 


from the leading edge and will produce, upon integra- 
tion, an increment of lift which evidently will be a func- 
tion of the pulse position ¢ and the geometry of the wing. 
The methods of solving for the influence of a single 
pulse are relatively simple and stem from application of 
the formal solution of Eq. (4) and formulas derivable 
from it. For given plan forms, the lift due to a unit 
pulse at ¢ can be expressed as F(t), the integral of pres- 
sure over the disturbed region. The total downwash 
is a continuous distribution of pulses of various strengths 
w(t)dt at various positions along the span, with the re- 
sult that, if the lift increment is written 


dL = [w(t) dt] F(t) 


it follows that the total lift is 


3/2 
) / F(t)w(t) dt (5) 
(6/2) 


The function F(t) is the influence function for lift; 
once it is known for a given plan form, the computation 
proceeds by determining the area under a single curve, 
a process similar to approximate strip theories. 


INFLUENCE FUNCTION FOR LIFT 


This computation is most easily carried out by use of 
an expression for pressure coefficient derived from Eq. 
(4). When the downwash is constant in the chordwise 
direction, as assumed here, Mirels® demonstrates that 
the difference in pressure coefficients between the top 
and bottom surfaces at a point (x, y) is given by a single 
integration 


- 


| I h(n)dn 
; (6) 
Un Je V(x — §)? — By — 7)? 


w, and where the 


AC. (x;y) 


where / (yn) = upwash velocity = 
typical paths are along the leading edge between e and 
fon Fig. 4. The unit pulse, w,éy = 1, results in zero 
upwash except for an extremely narrow zone 6n at 7 = 
t; at this point, £ = an = at. The radical is substan- 
tially constant over the range of ¢ to ¢ + én, and there- 


fore Eq. (6) can be evaluated at once as 


| a 
. (4) 
Ur LV (x — at)? — By — t)? 


The regions of influence of the pulse in the central and 
As a result of the spe- 


AC, (x, vy) = 


tip regions are shown in Fig. 5. 
cial integration limits in the tip area, there is a section 
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where the pressure cancels and the value of AC, is zero 


not that given by Eq. (7). The resultant lift is the 


influence function and is given by 


= ¢g / J AC,dxdy = 
eB 
4g si dxdy 
. (S) 
Un BV (x — at)? — B(y — t)? 


The plan 


F(t) 


Fig. 5 shows sample areas of integration. 


form characteristics—chord, span, and leading- and 


trailing-edge sweep will define the exact areas B for a 
given Mach Number and will thus determine the shape 
of the curve F(t). However, it is significant to note 
that in any case the areas B are simply Mach cones 
starting at the leading edge and spreading to the trail- 
ing edge; near the tip, the Mach line reflects back from 
the tip edge. Thus, if the flow were coming from the 
opposite direction, these areas would represent the 
This 
suggests that F(t) is also the solution to a reverse flow 


fore-cones of points on the downstream edge. 


Let it be assumed that there is a reverse flow over 
Refer- 
ence to the integral of potential of Eq. (4) and its areas 


the wing of Fig. 5 at the same Mach Number. 


of integration in Fig. 2 shows that the areas B of Fig. 5 
are fore-cone areas of integration for calculating a re- 
verse flow potential at the downstream edge. It is con- 
venient to note in this connection that spanwise loading 
is related to the wing potential at the trailing edge 
by 


dL “TE. 
= g AC,dx = 
dy L.E 


tg ' K Og / tq 
x*= P( Trailing Edge) 
: : |? mo ( i > FG l Ed 


A comparison may be made between the influence 
function, Eq. (8), and the spanwise loading for a uniform 
reverse flow with unit constant upwash. Using Eq. 
(4) with g(&, ») = 1 and using the coordinate system of 
Fig. 5 in which the downstream edge for reverse flow 1s 
defined by x = ay, the resulting span load is 


dL tq didn 
) 


oy ae (9 
dy iJ J, V (é — ay)? — Bn — y)? 


By comparison, the function dL./dy(y) is identical with 
F(t). 
loading of the reversed wing flying at a constant angle 
of attack. 
plan forms are available,’ and they can be used directly 


The influence function for lift is, then, the span 
The span loadings for wings of various 
to obtain the influence function, F(t). For generalized 
trapezoidal wings it is given below. The wing is split 
into three spanwise zones, as shown in Fig. 6, having 
associated with them, respectively, three formulas for 


the influence function. 
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LIFT, ROLLING MOMENT, AND PITCHING MOMENT Li 


For Zone I: 


1g j d (B? — d*)t 
F(t) = 1e- « — dt) cos! a “ 


Urv B? -— dad B(c — ai — dt) 
d (Bp? — d*)t ‘l) 


(10) 
8 B(c — at 4 it) |f 


(c — at + dt) cos | 
For Zone II: 
| 
: (11) 


Fy,(t) = (c — at — dt) 
UV B? -— @ 


For Zone IIT: 


ig f¢ — al — dt (8 + d)(b — 2t) 2(b — 2t)(c + dt -— db — at) \ 
cos~!| 1 — — (b- (12) 


= ~ 2t)*? 
Fun!) Ux V B? — ad? c—at—dt i-<¢ 


If the Mach lines from the tip and the root overlap, the influence function in the overlapping region is 


F(t) = Fy(t) + Fiyy(t) _ Fy (t) 


EXAMPLE OF LirT CALCULATION 
In general application, the lift is computed from Eq. (5) by multiplying the curve F(¢) by w(t) and measuring 


the included area graphically. However, for illustration, an example is chosen which can also be calculated 


analytically. Consider a rectangular wing of aspect ratio 2 or greater [for which F(t) is found by letting a = d = 0 
in Eqs. (10), (11), and (12)] in a field of linear symmetric upwash (Fig. 7). Thus, for the right half-wing, 


F(t) = Fy (0) = 4gc Ups 


: tgc {1 Bb Bt l b t Bb Bt\*) 
Fin (t) = — 4-— cos! |] 1 — — 2 + 2(8 — 28 — ( — 2 
UB lar C Cc 7 re Cc Cc t f 


and w(t) = [Uar/(b/2) |t 
The product of F(t) and w(t) is shown by the dotted curve in Fig. 8 and encloses an area that can be evaluated 


analytically to give 


9 *b/2 
C, = / F(t)w(t)dt 
gbc Jo 
lGar ((0/2-(/8 Gar (%/2 ty Bb Bt 
= tdt + m <cos~!| 1 — —2 + 
Bb Jo Bb? Jiws2) -(c/g) © | c Cc 
2(é - ws 7 - a Ay a 
c Cc Cc Cc { 


= (2a7/B8){1 — (1/BA.R.) + [1/2(BA.R.)?]} 


Lagerstrom and Graham’ treated this case, with the same result, by integrating the pressure coefficients. In 
addition, the actual span loading that they obtained is plotted in Fig. 8 for comparison with the product F(f)w(?). 
Although both yield the correct lift coefficient upon integration, only the true loading curve could be used to com- 
pute the bending moment about the wing root chord, since an attempt to treat F(/)w(t) as a load distribution curve 
Similarly, treating F(/)w(t) as a load 


obviously would produce an erroneous bending moment and is unjustified. 
To com- 


distribution curve in general will not correctly give the roll moment on a wing in a nonuniform stream. 
pute roll moments in nonuniform fields it is necessary to use a roll moment influence function analogous to the 


lift function. 
INFLUENCE FUNCTIONS FOR ROLL MOMENTS 


By similar methods it is possible to derive a rolling moment influence function R(¢) such that the total rolling 


b/2 
l= / w(t) R(t)dt (13) 
—b/2 


moment around the root chord of a wing is 
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Consider the contribution to rolling moment due to the unit pulse of Fig. 3. The pressure at each point within 
the disturbed region must be multiplied by its moment arm y and then integrated over the area. The pressures 
are given by Eq. (7). The contribution of a pulse at station / to the rolling moment is now 


R(t) = = ydxdy 
Ur B Vv ( (x — at)? — By — t)? U4) 


The integration areas are again those of Fig. 5. As before, this expression is to be compared with the span k vading 
of a reverse flow. This time the reverse flow contains a linear variation of vertical velocity along the span and is 
zero at the root with a slope of one—that is, g(&, ») = n in Eq. (4). Again using the coordinate system of Fig. 5, 
the span loading in reverse flow is 


dL Pa ndédn 
ay B ‘(— — ay)? — Bn — y)? 


Since the limits of integration are those used previously, Eq. (15) is the same function as Eq. (14). The flow de. 
scribed by Eq. (15), or the roll influence function, is equivalent to the span loading of a reverse uniform flow in 
which the wing is rolling at 1 rad. per sec. The span loadings on rolling trapezoidal wings have been calculated 
by Lagerstrom and Graham.’ 

Referring to the zones of Fig. 6 and the work of Lagerstrom and Graham, the influence functions for roll can be 


obtained. 


For Zone I: 


t Bt td d d Bt ad 
1—-(a+d)]}2 + (a —d) — - (1 + ) 
Ri(?) 2qbc Cc Cc Bc B B c B? 4 
— cos 
. UrB (Bb/c)[1 — (d?/B?)|” 1 — (t/c)(a + d) 
t Bt td ‘| d Bt ( 4 dt 2at t? | 
— —d _ 1 1 — 1 — — 2 — gq? 
E ; (a —¢ )|[2 8c (a +d) + 8 |e 2 ; 3? 3(b/2) " 2 (B a’) 
Lt 2 cos = 
(Bb c) [1 — (d?/B?)]” 1 — (t/c)(a — d) 1 — (d?/B?) | 
(16) 
For Zone IT: 
Ru(t) = 2 tle — t(a + d)] _ (d/B)[c — - sd | (17) 
UB\ V1 — (d?/p?) Bll — (d?/p?)]}'* § 
For Zone III: 
l Bt td d Z b 
1—-(a+d)}}2 + (a —d) — B+d —t 
lik an 2qbc c c Bc 8 oe) an c 2 4 
- Urp (Bb/c) [1 — (d2/B2) |” 1 — (t/c)(a + d) 
b 2 td d bd ] 
joe (8 —a) - (8+d)+1 (8+d)— (d +a) + _ (6B + d) 
2c Cc cB B Bc 
)(B — d) 1 — @/f*) J 
‘(6 are P 
2Bb I t Naa a 2c | a died (18) 
3¢ b/2 (b/2c)(B — d) 





If the Mach line from the tip and the root overlap, the influence function in the overlapping region is 


R(t) = Ri(t) + R(t) — Ru(t) 


INFLUENCE FUNCTIONS FOR PITCHING MOMENTS 


Proceeding in the same manner, the pitching moment about the 7 axis of Fig. 5 is expressed in terms of an in- 


fluence function P(t) by 


(b/2} 
m= / w(t)P(t)dt (19) 
~—(b/2) 





In 2 
the 


Thi 
moi 


7 


I 


Py 





1t within 


ressures 


(14) 


loading 
n and is 
f Fig. 5, 


(15) 


low de. 
flow in 
culated 


can be 


(18) 


n in- 


(19) 





LIFT, ROLLING MOMENT, AND PITCHING MOMENT 13 


Computing the pitching moment of a unit pulse gives 


tg an a x dx dy 
Py) = — I = (20) 
Ux B V(x — at)? — By — t)? 


1a uniform reverse flow the function P(t) compares with the span load of wings pitching at 1 rad. per sec. about 
I 


the downstream 7 axis. For this case g(£, ») = & in Eq. (4), and 


dL 1g a x tdédn 
inn / / (21) 
dy Ur J Je V(é — ay)? — Bn — y)? 


Thus Eqs. (20) and (21) are identified; this relationship is preserved for arbitrary choices of the axis of pitching 
moment. 
The influence function curves for pitching trapezoidal wings are: 


For Zone I: 


29c2 t t 2t t 
P(t) = — = ; aE —a--—d k (d-—a)+d@ (2 — i) x 
Un(B? — d?)’? | c cic c 


t t 
d{1l1—a — ~p 
c c t t | ft t 
cos~! +il—a-c+d — B?-(d+a)+d? {2a — | +4 
t ) rt C 4 C 


t 
B(1-a —@- 


Cc c 


t ai 
a(1—a )+ 6 
c 


- c l fic - t*) 
cos~! + 2d (1 —a ) V B? — ad? (1 —a ) -— 2 — (22) 
t Cc c cf 


l 
g(1-< +d ) 
‘ 


( 


For Zone IT: 


, t i}? 
7 a? (1 —a--d ) 
2qc"* c c (23) 


t t 
Py(t) = (a —a )(: —d-—-a ) _ - 2 
UV Bp? — ad? c Cc c c/ p? — d 


For Zone III: 


Qe? f . Cee f ; —— 
Putt) = ——— —, {jl —a--—d (B? — d*) |d- —a —d’?{l—a--d x 
Un( Bp? — d?) Cc Cc c c Cc c 


b t 
(8+ d) ( —2 . 
Cc Cc 1 /b t l ; , t t 
cos~!| 1 — + - —? (8 + d)|3{1— 2a (Bp? — d*) —({l—a--d > 4 
3 \c Cc Cc : 


t t > \e r ri 
l|—a--d 
ri Cc 
2 
— | 
’ t | (6 +d) b ,f 
(Bp? + 4dB) — B(B + d)(B — 24) ( —2 )| j ' a (24) 
Cc c} } 


The influence function in the overlapping region is 
P(t) = Py(t) + Pm(t) — Put) 


The axis of moments for these.formulas goes through the middle of the wing root chord. 
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DISCUSSION 


Wing lift, moments of roll, and moments of pitch are 
readily computed in downwash fields by Eqs. (5), 
(13), and (19), respectively. For any wing with super- 
sonic leading and trailing edges, the appropriate in- 
fluence functions have been identified with load dis- 
tributions in reverse uniform flows. The problem of 
thickness ratio, not considered in this report, can be 
treated independently by the linear theory and is the 
same as in the case of uniform flow; previous results 
for wave drag due to thickness are valid for nonuniform 


flows. 
The influence functions may be computed by applica- 


tion of Eqs. (8), (14), and (20). 
trapezoidal wings. 


They are given here 
for generalized For rectangular 
wings, the functions are illustrated in Fig. 9. 
they show, as expected, the effects of three-dimensional 
flows in the vicinity of roots and wing tips. 
parison with two-dimensional strip theories, the three- 
dimensional treatment shows marked differences in 
these root and tip regions, confirming the idea that low 


In general 


In com- 


aspect ratio wings that consist mostly of root and tip 
cannot be treated with any success by two-dimensional 


concepts. 


Inspection of influence functions for lift and pitching 
moment shows that it is possible for different distribu- 
tions of downwash to produce the same total lift but 
different moments. Thus, the spanwise distribution 
of downwash will have an effect on the chordwise center 
of pressure. From the point of view of tail efficiency, 
then, the reduction in lift due to downwash might be 
expressed in terms of one efficiency, while the moment 
efficiency would have a different value because of the 
shift in the center of pressure. This fact, which is al- 
ready known from experiment, is shown by analysis to 
be an inherent feature of nonuniform downwash dis- 


tributions flowing over planar tails. 
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The analysis, however, is restricted by the assump. 
tions of the linearized theory. If, as would be expected, 
the importance of wing thickness is similar to that jy 
the case of uniform flows, the pitching moment should 
be modified by a correction such as Busemann’s higher 
order approximation, whereas the lift and rolling mo- 
ment would be properly predicted for closed profiles by 
the linear formulas. 

The extension of influence functions to include sub. 
sonic leading and trailing edges has not been attempted, 
Such extension would introduce questions concerning 
the Kutta condition at trailing edges and the proper 
application of leading-edge suction and would make the 
problem in general more complicated. The results of 
an investigation, however, would be of interest to the 
designers of airfoils having highly swept plan forms. 
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The Interpretation of Failure Loads in the 
Plastic Theory of Continuous Beams and 
Frames 


P. S. SYMONDS* anp B. G. NEALt 


Brown University and Cambridge University, Respectively 


SUMMARY 


Plastic failure loads on continuous beams and rigid frames of 
ductile metal can be simply computed when certain assumptions 
concerning the moment-curvature relations are made and when it 
is assumed that individual loads remain in constant ratio to each 
other during the loading process, so that the load value is de- 
fined by a single parameter. The simple theory is mainly con- 
cerned with the prediction of the maximum load a structure can 
carry before ‘‘plastic collapse’ occurs. However, in some cases it 
will be desirable to estimate some of the deflections, and a case of 
a continuous beam is described herein in which the calculation of 
plastic failure load may be misleading unless accompanied by a de 
flection analysis. The deflection values that make the plastic 
analysis meaningful are calculated for this case by a simple ap 
proximate method that requires little additional labor after the 
plastic failure load has been determined. The results are com 
pared with test results previously published by Sttissi and Koll- 
brunner!! and by Maier-Leibnitz."' '? The method may be used 
for general framed structures, and an analysis of a two-story 
portal frame is given as a further illustration of its use 


(1) INTRODUCTION 


: ea PLASTIC METHODS OF ANALYSIS and design for 
beams and frames of ductile metal'~’ are mainly 
concerned with the load at which failure by “‘plastic 
collapse’ would occur. The basic assumption of these 
methods is that the bending moment at any cross sec- 
tion cannot exceed a value termed “‘the fully plastic 
moment.’’ This value corresponds to the spread of 
plastic zones across the entire cross section, and, as it is 
approached, large changes in curvature can occur, with 
negligibly small changes in bending moment. A “‘plastic 
hinge’’ is said to occur at a cross section where the fully 
plastic moment is developed, since it is assumed that 
hinge action can then take place at this cross section 
while the bending moment remains constant. 

Under these hypotheses, when the loads on a framed 
structure are steadily increased in fixed proportion to 
each other, a critical load is finally reached at which the 
structure would become a mechanism owing to the 
formation of plastic hinges at a sufficient number of 
sections. It has been shown*~ how this critical load can 
be easily calculated as compared with elastic analyses 


for highly redundant structures. 
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Tests,' * ° as well as theoretical considerations, show 
that in actual structures a certain amount of strain- 
hardening always takes place at the cross sections where 
plastic hinges form in the theory, so that the load can 
be carried slightly above the theoretical collapse value. 
In practice, complete collapse, with a drop or total loss 
of load-carrying capacity, occurs only when additional 
effects arise, such as lateral instability of beams, rupture 
of welds, additional moments due to large deflections, 
etc. These effects ordinarily appear only after large 
plastic deformations have been produced, at loads ex 
ceeding the theoretical collapse value. The collapse 
load of the plastic theory then gives a good estimate, 
generally conservative, of the observed load at which 
small load increases produce much larger deflections 


than at previous load levels. 


Although most of the experimental investigations 
have dealt with mild steel structures, the methods 
should be useful in cases of structures of other metals of 
sufficient ductility. The effects of strain-hardening 
should be qualitatively the same: The actual structure 
does not collapse at the critical load of the theory but 
carries loads more or less higher than this, depending on 
the rate of strain-hardening. Further experimental in- 
vestigations are needed; questions of failure by rupture 
or by buckling, among others, must be investigated for 
the metals that have not yet been studied. 


The plastic methods are directly concerned with 
loads rather than with deflections, although the criterion 
of failure is the imminence of large plastic deflections. 
In cases where the precise values of the permissible de- 
flections themselves are the primary quantities that 
control a design, the plastic theory will probably not be 
applicable. The exact calculation of deflections in a 
partially yielded redundant structure is a formidable 
task, especially if measured stress-strain characteristics 
are used. Even for simple beams, fairly lengthy calcu- 
lations are required,’ while, for redundant structures, 
trial-and-error procedures must be carried out based on 
extensive preliminary computations.’ Except for re- 
search purposes, such calculations would never be justi- 
fied, however. They are necessarily based on assump- 
tions as to constancy of physical properties, such as the 
yield stress, the strain at which strain-hardening begins, 
etc., and as to conditions of end-fixing, rigidity of sup- 


15 
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P (1/4)PL2 = Mn — (1/2)Mz — (1/2)Mc = Mn — Me 
A B z' c D (1) 





(a) PP >> —S«s Ms > ke 


(b) 


(c) 





Type of continuous beam tested by Stiissi and Koll- 
brunner!! and Maier-Leibnitz.!? 


Fic. 1. 


ports, absence of residual stresses, and the like. Many 
tests have shown that in actual structures these 
conditions may vary by large amounts quite unpre- 
dictably. Moreover, if elaborate calculations to de- 
termine deflections had to be made, the great ad- 
vantages of the simplicity and directness of the plastic 
methods would be lost, and conventional elastic methods 
would be better employed. 


In many cases, however, fairly rough estimates of de- 
flections either at the critical collapse load or at the 
normal working loads may be desirable, or necessary, 
and it would be useful to have for these cases a method 
of estimating important deflections which would require 
only a small amount of work beyond that needed for the 
usual plastic analysis. The present paper shows how 
such deflection estimates can be easily obtained, with 
some indications as to their probable accuracy. 


(2) ConTINUOUS BEAMS 


As an example in which deflection estimates are essen- 
tial in order to make the plastic analysis fully meaning- 
ful, consider the continuous beam and loading shown in 
Fig. la. Tests on this type of structure were reported 
by Stiissi and Kollbrunner,'' whose main purpose was 
to investigate the effect of changing the ratio L/L. on 
the observed failure load. The concept of failure load 
used by Stiissi and Kollbrunner was not the plastic col- 
lapse load derived from the fully plastic moment values, 
which is the basis of present-day methods. They were 
concerned with loads at which complete collapse oc- 
curred rather than with the load at which a given load 
increment first begins to produce much larger plastic 
deformations than previously. 


Let the fully plastic moment of the beam (assumed of 
uniform section and material) be 1/,. Then the load P 
has reached the plastic collapse value P, when the bend- 
ing moments at the midpoint m and at supports B and 
C have reached the magnitude /,. To calculate P,, 
consider the moment equilibrium equation for the 
middle span under the given symmetrical loading, 


Now, since M, < M,, M,. > — M,, it follows that 


PL, < 4(M, + M,) = 8M, (2) 
The critical value P, is therefore given by 
P, = 8(M,/L2) (3) 


This is the highest load that can be applied without 
violating either the equilibrium conditions or the re- 
quirement that the bending moment shall not exceed 
the fully plastic moment in magnitude at any section. 
If the outer spans were removed, leaving the central 
span simply supported, the collapse load P.’ would be 


P = 4M, Le (4) 


Apparently, therefore, the effect of attaching the outer 
spans is to double the load at which plastic collapse 
occurs according to the plastic theory, and this result is 
true independently of the length L;. However, if L; is 
made extremely large compared to Ll», the middle span 
becomes effectively simply supported. There is ap- 
parently an inconsistency in the plastic theory, which in 
the one case predicts a failure load 4M,/L»2 and in the 
other case predicts twice this value for effectively the 
same simply supported beam. This was pointed out by 
Stiissi and Kollbrunner, who cited their test results to 
show that the failure load of the continuous beam with 
large ratios of L,/L» was less than double the failure load 
of a beam without outer spans. They concluded, there- 
fore, that the “‘moment equalization’? method, as they 
called it, was not sound, since it led to an overestimate 
of the failure load in this case. 

There is, in fact, no inconsistency in the plastic 
theory. According to this, if there are outer spans of 
any finite length, then the rate of increase of deflection 
with load is finite, since it is limited by the elastic con- 
straint of the outer spans. This is true until the 
moments at the supports B and C reach the negative 
fully plastic moment, when theoretically the deflection 
begins to increase at an infinite rate. 

The difficulty in the present example is that if the 
ratio L/L. is made extremely large, then a large deflec- 
tion of the central span will occur before the moments 
Mpzand M_, reach the fully plastic value. Normally, the 
plastic failure effect is observed as a transition from 
moderate to much higher rates of deflection with load. 
In the present case, if L,/L2 is made large, this would 
become increasingly difficult to observe. It is evident 
that the theoretical collapse load may not, in this case, 
give a reliable indication of the limiting useful load 
that the structure can carry; for, if L;/L2 is large 
enough, tolerable deflection values would probably be 
exceeded before the theoretical failure load is reached. 
The usual plastic analysis must therefore, in this case, 
be accompanied by an estimate of important deflections. 
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INTERPRETATION 


To estimate the central deflection 6,, in the present 
case at the instant when the critical load has just been 
reached, we note that continuity of slope still exists 
across any section until the full plastic moment is de- 
veloped there. Therefore, at a load infinitesimally be- 
low the critical value there is continuity of slope across 
the section at C, since it is obvious in this case that the 
fully plastic moment develops first at the midpoint m and, 
finally, at C, leading to plastic collapse. Free-body dia- 
grams of the lengths mC and CD are shown in Fig. Ic, 
with the moments and deflections as they would be just 
as the failure load is reached. 

The following two assumptions will be made: 

(1) Each span retains its original flexural rigidity // 
except at the cross sections where fully plastic moments 
are reached. 

(2) Rotations occur freely at constant moment values 
at the cross sections where the fully plastic moments are 
developed. 

The first of these amounts to neglecting the additional 
flexibility caused by the spread of plastic zones along the 
length of the members. The second states that the 
strain-hardening at sections where plastic hinges are 
formed is to be ignored. The errors caused by these 
two assumptions tend to cancel each other, so that the 
calculation of deflections based on them may often be 
accurate. This calculation can be made directly from 
the collapse analysis without requiring successive elastic 
solutions for the beam as various plastic hinges have 
formed. 

The last plastic hinge having formed at C, continuity 
of slope just as the collapse load is reached requires 
that 


Bem = Ben (5) 


these being the clockwise angles of rotation at the ends 
at Cof the spans Cm and CD, respectively. The slope- 
deflection equations* for these spans are 

SEL bm 


M, —- (6) 
es L.? 


2EI l l 
Le bcm = 5 M, —- 


(EI/11)0cn = —(1/3)M, (7) 


After equating the expressions for 6¢, and @¢p and 
solving for 6,,, we obtain 


5, = (M,L2?/24ET) (1 + 4(L1/L2) | (8) 


This result shows at once that the plastic collapse 
load calculated for the extreme case in which L;/L. = © 
has no practical significance, since the deflection 4,, 


* The general form for a span PQ is as follows: 
(EI 1)(Op@ om Opq’ — 6 1) = (1/3) Mpg — (1 6) Mop 


where EJ = flexural rigidity (constant), / = length of span, Ope 
= slope angle of P, @pg’ = slope angle at P which given span 
loads would produce if the span were simply supported at P and 
VY, 6 = relative displacement of ends, and Mpg, Mgp are end 
moments at ?, Q, respectively. Signs of all slope angles, 6//, and 
all end moments are defined by the positive-clockwise convention 
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Kollbrunner!! on continuous beam of dimensions shown. P- 
is collapse load and 6, is approximate deflection at collapse, 
determined by simple plastic theory. 
would also become infinite. This limiting case, how- 
ever, is of interest mostly as a mathematical curiosity. 
For spans of reasonable length ratios the failure load 
and central deflection values of the simple plastic 
theory have their usual significance of marking the 
transition from moderate rates of deflection with load 
to much higher deflection rates, with useful accuracy. 

Stiissi and Kollbrunner"! tested steel beams with span 
ratios L;/L: equal to 3, 2, 1, and '/2, in addition to a 
simply supported beam. In all cases Ll. was 60 cm. 
However, the only curve for load plotted against central 
deflection reproduced in their paper is that for the case 
L,/L. = 2. The measured points given for this case 
(beams 532/6 and 534/8) are shown in Fig. 2a, together 
with the point corresponding to the calculated collapse 
load and central deflection according to the simple 
plastic theory. No value of M, is given in the paper; 
however, a sufficiently close estimate is obtained by 
multiplying the yield moment M, by a shape factor 
1.10 appropriate to the rather light /-section used. 
(The dimensions given are shown in Fig. 2a.) 


M, = 1.10M, = 1.100,Z (9) 


Fig. 2a shows how P, and 6,, given by the simple theory 
appear in relation to the curve of load vs. midpoint de- 
flection measured by Stiissi and Kollbrunner. The 
point with coordinates (P., 5) does not fall on the 


measured curve. However, the value P, gives a fair 
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estimate of the plastic failure load, and 6,, indicates 
roughly the magnitude of the deflection at the load at 
which large plastic deformations begin to occur. Stiissi 
and Kollbrunner gave the total failure load for this 
beam as P,, = 3.90 tons, without stating precisely what 
happened at this load except that “deformation = ©.” 
Since the usual effect of strain-hardening is to cause the 
load-deflection curve to rise steadily until an additional 
effect, such as rupture or buckling, causes a drop in the 
load-carrying capacity, it is difficult to interpret the 
meaning of P, without supposing that one of these 
effects occurred. The failure load P,, with which the 
plastic theory is properly concerned, is conservative, 
since in almost all cases slightly higher loads can be sup- 
ported. 

Further tests on a beam loaded as shown in Fig. la 
and with L,/L, = 2 were made by Maier-Leibnitz.!” 
The beam tested was of steel, with area A = 25.8 cm.?, 
section modulus Z = 90cm.', J = 453 cm.!, Le = 60cm., 
E = 2,035 tons per cm.?, a, = 2.455 tons per cm.”. The 
latter is the lower yield stress obtained by averaging 
values from four specimens cut from flanges. The 
value of M, calculated from the cross-section dimen- 
sions (see Fig. 2b), and the above value of o, is 248 ton- 
cm., corresponding to a ‘‘shape factor” of 1.12. With 
these values the Eqs. (3) and (8) yield the results 


P. = 16.5 tons, 6 = 1.45: em. (10) 


This point is plotted in Fig. 2b for comparison with the 
measured deflection curve. The beam used in these 
tests was strengthened at the supports and load point, 
and Fig. 2b shows that this had considerable stiffening 
effect. If the value for /, of 262.5 cm. is used, as sug- 
gested by Maier-Leibnitz on the basis of tests ona 
simply supported and centrally loaded beam, the calcu- 
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lated values are P, = 17.5 tons, 6, = 1.54 cm. This 
point lies close to the experimental curve. However, it 
seems advisable to calculate /, from the measured 
yield stress and section dimensions in assessing the 
simple plastic theory. Maier-Leibnitz did not carry the 
test to complete failure and therefore did not obtain a 
value of P,, corresponding to that of Stiissi and Koll- 


brunner. 


As a second example of the estimation of deflections 
at the collapse load of the plastic theory, consider 
another continuous beam for which Maier-Leibnitz 
gave test results.'_ In this reference tests are described 
on the beam shown in Fig. 3a. The supports are at 
equal heights. For this steel /-section, the author gave 
the value of /, as 660 tons per cm., and further data as 
L = 80 cm., J = 1727 cm.‘, E = 2100 tons per em.‘ 
According to the plastic theory collapse occurs when 
the magnitudes of the moments at A and / reach the 


value M,. The moment equilibrium equations are 


M, = PL + (1/3)M, (11) 
M, = PL + (2/3)M, (12) 


From Eq. (11), we obtain 
PL = M, — (1/3)Ma4 < (4/3)M, (13) 
Corresponding to the upper limit 
P, = (4/3)(M,/L) (14) 
the set of bending moments are 


M,\M;| Mz Mn 
— M,|M,|(2/3)M,|(5/6) My 


(15) 


The theoretical collapse load is, therefore, as given by 
Eq. (14). 


To estimate the deflection at the collapse load under the same simplifying assumptions made previously, the 


free-body diagrams of Fig. 3b are used. 
tinuity of slope across the middle section requires that 


Ain B 


It is assumed that the last plastic hinge to form is at section 1. Con- 


4 (16 


mA 


These angles are given by the appropriate slope-deflection equations, in view of the assumed continuity of slope at 


section |, as 


2 EI P I (° Vu ) 5 Pl kl : 2) 
m is sul p > cls Om \ 
“Eis 3\6° 7") 54 ((3/2)L}? 
2 El P 1/5 vf ) n 5 PI kl ; ' 
m aaa 4 ? cle ™ On, (IS 
sit” se 54 [(3/2)L]? 
The value of 6,, at collapse is found to be Inserting numerical values, the results are 
bm = (103/144) (7,L7/ET) (19) P, = 11.0 tons; bm = 0.833 cm. (20 
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INTERPRETATION 


The point with these coordinates is plotted in Fig. 4 for 
comparison with the measured load-deflection curve ob- 
tained by Maier-Leibnitz. It will be seen that in this 
case the calculated point locates accurately the load and 
deflection values when large plastic deflections are im- 


minent. 
(3) RiGiIp FRAMES 


The method of estimating deflections used in the 
above examples may be applied to portal frames and 
other types of rigid frames such as that shown in Fig. 
ja. It is especially advantageous in these cases because 
of the high degree of redundancy usually involved. 

The method used depended on an assumption as to 
which plastic hinge of all the plastic hinges that occur 
in the collapse mode of the structure was the last to form. 
In the beam examples taken, the choice of the last plas- 
tic hinge to form was obvious. However, in more com- 
plicated frames, it is rarely obvious which hinge will 
form last, although usually it is easy to guess where the 
first few plastic hinges will occur. The procedure is to 
assume, in turn, that each of the hinges is the last to 
occur and, by applying the slope-continuity condition 
to this hinge, to calculate the corresponding value of 
some particular one of the deflections of the frame. 
The choice of hinge that yields the largest value of that 
deflection is the correct one. The reason for this is that 
imposing a false continuity condition amounts to re- 
moving a ‘“‘kink’’—-i.e., a discontinuity in slope which 
was produced at that section in the course of bringing 
the loads up to the collapse value. However, the only 
way such a kink can be removed, while maintaining the 
values of the moments that actually occur at the col- 
lapse load, is by a motion of the frame as a mechanism 
in the reverse direction to that which would occur in the 
actual mode of collapse. Thus, when a condition of 
slope continuity is assumed across any hinge other than 
the one actually formed last, the resulting value of the 
deflection may be considered as the true value dimin- 
ished by a certain amount due to a backward motion of 
the frame as a mechanism. No such backward mecha- 
nism motion is imposed when the correct choice of hinge 
is made, and this corresponds, therefore, to the largest 
value determined by the several choices. 

To illustrate by a simple example, consider again the 
continuous beam of Fig. 1. If the last hinge to form is 
assumed (erroneously) to be the one at midpoint, then 


the condition 6,,, = Onc 1s sufficient to determine a 
value 6,,’.. The pertinent slope-deflection equations 
are 
2E1O nn ti] : 
ows —_ 5 (21) 
Lo 6 L.? 


216 nc | 


( *) 
=-M,-—- * (22) 
Ly 6 L,* 


Thus, the assumption of continuity of slope at m yields 


the result 
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Fic. 2b. Curve shows deflections measured by Maier-Leib- 
nitz on beam of dimensions shown. /, is collapse load and 6,» is 
approximate deflection at collapse, from simple plastic theory 
(ignoring stiffeners actually present at load points). 
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Fic. 3. Type of continuous beam tested by Maier-Leibnitz.' 


6, = (1/24)(L.2M,/ED) (23) 


Comparing this with the result for 6,, given by Eq. (8), 
the above value of 6,,’ is seen to be smaller than 6,, for 
any value of 1; greater than zero. The agreement be- 
tween the two results for the case L; = 0 indicates that, 
when the beam is clamped at B and C, both plastic 
hinges form simultaneously. 
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Fic. 4. Comparison of collapse load and deflection at col- 
lapse according to simple plastic theory with measured load-de- 
flection curve obtained by Maier-Leibnitz.! 
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Rigid frame and loads used to illustrate deflection 
calculation. Fic. 5b Mode of collapse; black circles show plastic 
hinge locations. Fic. 5c Free-body diagrams used in calculation 
of deflection (shear reactions omitted for simplicity). 


Fic. 5a. 


If the result Eq. (23) is used to compute the discon- 
tinuity in slope across the joint C (or B), we find 


Aen — bom — (1 3) '(L,M, /FI) (24) 


The negative sign indicates slope angle values at section 
C corresponding to a concave upward shape of the beam. 
This absurd result is due to having made a false assump- 
tion as to continuity of slope at the midpoint. As 
pointed out above, this false assumption corresponds to 
removal of a “kink” that actually existed at m by 
means of a ‘“‘backward’’ mechanism motion of the span 
BmC. 


sult from Eq. (24) by imposing a “‘forward’’ mechanism 


This suggests that we can obtain the correct re- 
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motion of magnitude just sufficient to cancel the dis- 
continuity given by Eq. (24). In this motion the outer 
spans are held rigid, all motion occurring at the hinges 
at B, m, and C. The resulting angle change at C, 
— Adem, is set equal to (1/3)(11M,/EI). 
sponding increment of deflection A6,, is given by 


i | LiL2M, 
( = Adem ) = = 
2 6 El 


The corre- 


Abn = (25) 
The result of adding this to the value of 6,’ given by 
Eq. (23) is 


(26) 


: Le? Lyle\ M, 
bn + 46, = | —- + — 
24 6 EI 
This is identical with the correct result for 6, given by 
Eq. (8). The above process provides a check on the 
value obtained for 6, which is useful in more compli- 
cated problems. 


As an example of deflection calculations for a more 
complicated frame, the two-story portal frame shown 
in Fig. 5a will be considered. The vertical loads are 
taken as W and the horizontal loads as 0.90W each. 
All members are taken to have the same fully plastic 
moment M,. The mode of collapse under these circum- 
stances is as indicated by Fig. 5b, the black circles 
representing plastic hinges. There are six redundant 
moments in this structure, and the mode of collapse in- 
dicated in Fig. 5b has only six plastic hinges. This is a 
case of incomplete or local collapse,*: * [see Section (5) of 
this paper] in which the elastic moments are not com- 
pletely determined by the values of the moments at 
plastic hinges and conditions of statical equilibrium. 
However, the elastic moments can readily be calcu 
lated, the method of least work being especially con- 
venient for this purpose. It is found that there are four 
undetermined elastic moments that are related by three 
equations of equilibrium. The values of these moments 
san be obtained by minimizing the strain energy of the 
frame with respect to one of the four moments, taking 
proper account of the equilibrium equations but re- 
garding the remaining eight moments as constant. The 
complete set of moments thus obtained is given in 
Table 1. 
is indicated by the dashed lines in Fig. 5a; 


The moment sign convention used in Table | 
positive 
moments are ones that cause tensile stress in the side of 
the beam adjacent to this line. 

The collapse load W, is found to be 


W,. = 2.13(M,/L) (27) 
TABLE | 

Section A F G I K i 
Bending Moment 
ei Sim : 4 = Ere 

M, 
Section B Cc D E H J 
Bending Moment a. ' - - ” 

7; 0.90 O87 0.038 0.87 0.07 —0.93 

M> 
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INTERPRETATION 


The details of the analysis to determine the collapse load 
and the bending moments at collapse need not be given 
here; references 3-6 may be consulted if necessary. 

The first step in the deflection analysis is to draw the 
free body diagrams of the members of the frame, as 
shown in Fig. 5c. Using these, one can then calculate 
any one of the deflections, such as the side-sway 6, at 
joint B, from an assumption of slope continuity across 
any one of the plastic hinges in the collapse mode. 
Suppose, for example, that the hinge at L is assumed 
to be the last one to form—i.e., that @,¢ = 6,;.. Writing 
the slope-deflection relations for members CL and L/, we 


obtain at once 


l 0.87 l | 
——- M, — —M,+ 4, == M,—-=M, — A 
3 6 3 6 , 
(28) 
whence 
A, = (EI/L*) 6, = 0.323M, (29) 


Now, since sections B and C are not plastic hinges, we 
know also that 6¢, = 024 and can write, therefore, 


0.90 l 


0.87 
M, + ; M, + Ap 
) 


I 
rm M, + 6 M, + AL = 
0 ) 


(30) 


whence 
Ap = (El/L?*)bg = 0.913M, (31) 


An assumption of slope continuity at each of the other 
five hinge locations leads, after a similar calculation, toa 
value of A, [defined as in Eq. (31)]. Some of the cal- 
culations are shorter than that just given; others are 
somewhat longer, but none takes more than a few 
It is unnecessary to give the full 
The choice 


minutes to work out. 
details. The results are set out in Table 2. 
of F as the last plastic hinge to form leads to the largest 
value of 6,. It is therefore concluded that this is the 
correct choice and that the value of the side-sway at 
joint B when the collapse load is just reached is given by 


bp = 1.363(L2M,/EI) (32) 


As a check, the alternate procedure used in the pre- 
vious beam example may be applied. The (incorrect) 
choice of L as the final plastic hinge allows the slope dis- 
continuities to be calculated at all of the hinge loca- 
tions. For example, the rotation angles at F are found 


to be given by 


Orn = — 0.605(LM,/ET) (33) 
Oro = — 0.294(LM,/ED) (34) 
TABLE 2 


Hinge where slope 
continuity assumed 


F L G I A K 


EI . = 
Result for 5 bz 1.363 0.913 0.908 0.450 0.183 0.178 


2Mp 
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From these, the slope discontinuity (in the sense indi- 
cated in Fig. 5b) is found to be 


Ore — Ove = —0.899(LM,/ET) (35) 


The negative value of this means, of course, that the 
choice of L as the location of the last hinge to form is 
wrong. However, we can now imagine a motion of the 
frame as a mechanism in the actual mode of collapse of 
just such an amount as to remove this slope discontinu- 
ity at F. Let the clockwise rotation of leg AE be y; 
then the relative rotation at F is 2y, and this is chosen 


so that 
2y = 0.899(LM,/E1) (36) 


Corresponding to this mechanism motion, there is a 
side-sway displacement at B equal to 


YL = (0.899/2)(L?M,/EI) (37) 
Adding this to the result for 6, given by Eq. (31), we 
obtain 
4 ow 0.899\ L?M, ong LM, 
62 = (0.913 + 9 EI 1.363 EI 


(38) 


which agrees with the result obtained previously by 
assuming the last plastic hinge to be that at F. It may 
be checked that all of the slope discontinuities are now 
positive—i.e., in the sense shown in Fig. 5b. 

It may be noted that a consideration of deformations 
in the elastic structure will usually indicate which are 
the most severely stressed sections. In the present exam- 
ple the joints at A, AK, J are obviously highly stressed 
initially. It would have been a safe guess to assume at 
the start that none of these is the last plastic hinge to 
form; this would have eliminated about half of the work 


involved in the calculation. 


(4) AccURACY 


As pointed out already the present method of esti- 
mating deflections depends on assumptions that amount 
to neglecting (1) additional flexibility due to spread of 
plastic zones along the lengths of members and (2) 
stiffening due to strain-hardening at joints where plastic 
hinges are assumed to occur. 

The accuracy of the method depends on the extent to 
which these two effects cancel each other. In the cases 
of continuous beams used as examples in this paper, the 
bending moments changed in sign at successive joints. 
This meant that the plastic zones were rather localized, 
and the additional flexibility due to them was nearly 
canceled by the stiffening effect of strain-hardening. In 
other cases some members may be in pure bending, or 
nearly so, and, if plastic flow occurs in these, the actual 
deflections may be considerably larger than those pre- 
dicted by the present method. 

The comparisons shown in Figs. 2 and 4 give some in- 
dication of the accuracy of the method as applied to 
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continuous beams. Comparisons between deflections 
calculated by the simple theory and those observed in 
tests on small-scale portal frames have been given in 
another place.'* The test results were obtained by 
Baker and Heyman* using portal frames of 4-in. span 
and 2-in. height and loaded by a horizontal load at the 
beam height and a vertical load at the midpoint of the 
span. Since the load ratios in these tests were chosen to 
cause collapse in all the possible modes, they indicate 
what may be expected for portal frames of rectangular 
cross section under a wide range of conditions. The 
comparison between observed deflection curves and cal- 
culated points at collapse shows that in many cases the 
agreement is extremely close. In other cases the calcu- 
lated deflection is considerably too large and in others 
too small. However, the agreement in general is good 
enough so that, in view of its simplicity, the »»ethod may 
be considered a useful means of estimating deflections 
at the theoretical collapse load of a structure. Such 
tests as these on model frames must, however, eventu- 
ally be supplemented by similar tests on full-scale frames 
made from standard commercial sections using normal 
fabrication procedures. 


(5) SPECIAL CASES 


Three types of collapse under proportional loading 
can occur, depending on the structure and the ratios of 
Let 1 be the number of redundant moments 
The type of collapse which most 


the loads. 
in a beam or frame. 
investigators have considered'~‘ is that in which » + | 
plastic hinges have formed just as the collapse load is 
reached; this may be called ‘‘complete”’ collapse. The 
beams used as examples in this paper failed with com- 
plete collapse; in both cases n = 1, and there were two 
However, collapse with fewer than 
This 


hinges at collapse. 
n + | hinges will often occur in practical cases. 
has been termed “‘local’’ collapse or preferably ‘‘incom- 
plete’’ collapse. For example, the two-story frame of 
Fig. 5a fails, under the loading shown, with six plastic 
hinges as in Fig. 5b. In this structure » = 6, so this is a 
case of incomplete collapse (but the term “‘local’’ seems 
inappropriate). Finally, a third type of collapse is 
possible, although in ordinary structures under normal 
loadings it is less likely than the other types. This is the 
type of collapse in which several hinges form simul- 
taneously (in theory) just as the critical failure load is 
reached, so that in the collapse mode there are more 
than n + 1 plastic hinges. As an example, consider 
the two-story frame of Fig. 5a, but let the two trans- 
verse loads each be JV instead of 0.901V. Then the 
mode of collapse is actually one with ten plastic hinges, 
only the moments J/p and 7, being elastic. This sort 
of failure might be termed ‘‘overcomplete”’ collapse. 
(Alternatively, the terms “indeterminate,” 
nate,’ and ‘“‘overdeterminate’’ might be used to specify 
the three types of plastic collapse.) 

In the rather special cases of overcomplete collapse, 
the present method for estimating deflections may be 


““‘determi- 
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difficult to apply, since in order to calculate one of 
the deflections it will often be necessary to assume Slope 
continuity at a group of hinge locations rather than at 
merely one. Obviously, there may be many more pos. 
sible groups than there are plastic hinges. In some 
cases, however, such as the portal frames whose test re- 
sults and calculated failure points are discussed in 
reference 13, the occurrence of more than n + 1 hinges 
does not present any difficulties that cannot be easily 
resolved by the use of a little judgment concerning the 
probable behavior of the frame. However, in a more 
complicated case, such as the two-story frame cited, 
the situation may be awkward. Probably the simplest 
expedient is to alter slightly the ratios of the applied 
loads so that the mode of collapse will be different and 
of either the complete or incomplete type. In either of 
these cases there will be no difficulty in making the de- 
flection analysis as described in this paper. The result- 
ing deflection values may either be regarded as suf- 
ficiently accurate in themselves or may be taken as a 
guide in investigating deflections under the original 
loading. 
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Simplified Treatment of the Turbulent 
Boundary Layer Along a Cylinder in 
Compressible Flow 


HANS U. ECKERT* 
Wright Air Development Center 


(1) SUMMARY 


An attempt is made to determine the effect of lateral surface 
curvature upon the turbulent boundary layer of a circular cyl- 
inder immersed in a compressible fluid with its axis parallel to the 
direction of flow. The calculation is carried out by means of the 
momentum theorem with the assumptions that the velocity pro- 
file (1/7-power law) and an empirical relation that puts the wall 
shear stress inversely proportional to the fourth root of the bound- 
ary-layer thickness (Blasius Law) are unaffected by the curva- 
Further assumptions are constancy of static pressure and 
On this basis the following results are obtained: 


ture. 
total energy. 
Full thickness and displacement thickness are less than on the 
flat plate, while momentum thickness and local and mean fric- 
tion coefficients are higher. For the displacement thickness the 
deviation from the flat plate value may become appreciable al- 
ready if the ratio of full boundary-layer thickness to cylinder 
radius is less than 0.1 and for the full thickness if this ratio is 
somewhat above 0.1. The friction coefficients on the other side 
are rather insensitive to curvature and deviate appreciably only 
if the ratio becomes as large as one. Compressibility in general 
somewhat increases the effect of curvature except for the displace- 
ment thickness where a slight decrease is found. The reliability 
of the basic assumptions and effects of their variation are dis- 
cussed 


(2 


<) 


INTRODUCTION 


¢ A PRECEDING PAPER,' an approximation method 
has been presented for treatment of the turbulent 
boundary layer on a flat plate in compressible flow 
based on von Karman’s momentum theorem and the 
assumptions of constant static pressure, constant total 
energy, and a velocity distribution according to the 
\/7-power law. A shear stress relation involving the 
same assumptions has been derived from tests on com- 
pressible flow in pipes and applied to the case of the 
plate in a manner that is consistent with earlier success- 
ful procedures in the incompressible régime. 

In the same way as it has been shown in that paper 
already on some parameters for flow in pipes, it can be 
proved generally for cases of rotational symmetry that, 
if the velocity profile follows a power law, boundary- 
layer parameters (as, for instance, mean values for ve- 
locity density and temperature or momentum and dis- 
placement thickness) can be expressed by corresponding 
parameters for two-dimensional flow. These consist 
of the ordinary flat plate term based on the same ve- 
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locity profile as it occurs on the rotational body and a 
second term based on a profile with half the exponent 
of the common profile, which may be called second-order 
The second term takes account of the lateral 
It is negative for axial flow inside 


profile. 
surface curvature. 
cylindrical surfaces (concave curvature) and positive 
for flow outside such surfaces (convex curvature). 

In the present paper the influence of convex surface 
curvature will be investigated. This case has some 
practical significance for boundary-layer experiments 
in wind tunnels where, though ultimately data for two- 
dimensional flow are desired, cylindrical surfaces are 
often preferred to flat plates in order to avoid the margin 
effects. In reverse, the problem emerges if two-di- 
mensional data are to be applied to missiles of large 
fineness ratio. In cases of model tests with large dif- 
ferences in Reynolds Number an additional scale effect 
produced by the lateral curvature may have to be 
taken into consideration. Finally, the problem is of 
some importance for determination of the base pres- 
sure on projectiles since a close relationship between 
this quantity and the boundary-layer flow has been dis- 
covered in recent years.* * 

It can be anticipated that the behavior of the bound- 
ary layer along a cylinder will deviate appreciably 
from that on a flat plate if the boundary-layer thick- 
ness becomes comparable in size to the cylinder radius. 
The question is, however, in which way the boundary 
layer adapts itself to the changing condition. 

Without a detailed knowledge of the mechanism of 
turbulence one can suppose merely that the following 
two cases form limits within which the actual behavior 
may occur. One of these probable limits is that the 
shear stress at the surface, and thus according to the 
momentum theorem also the momentum defect of the 
fluid, remains unaffected by the lateral curvature. 
For equal momentum thickness it follows, out of the 
cylinder geometry, that the full thickness of the 
boundary layer will be smaller on the cylinder than on 
the plate if the simplifying assumption is made that 
the velocity distributions within the boundary layer 
are the same for cylinder and plate. If this limiting 
+ The existence of this relation has been indicated in an earlier 


paper by Tetervin.? 
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Fic. 1. Nondimensional momentum thickness of second-order 
velocity profile related to same value of first-order profile vs. free- 
stream Mach Number for various power law exponents of the 
first-order profile. Values refer to flat plate flow. 





case would be true exact application of friction (and 
heat transfer) coefficients from flat plates to cylinders 
and vice versa would be possible. 

The other limit may be that the full thickness of the 
boundary layer remains constant. For same velocity 
distributions this means higher momentum defect 
and thus higher shear stress on the cylinder. On this 
basis Jakob and Dow’ made an estimation for incom- 
pressible laminar and turbulent boundary layers to 
find out how closely their measurements of heat trans- 
fer from cylinders may represent flat plate values. For 
turbulent boundary layers they assumed a velocity 
distribution according to the 1/7-power law and ob- 
tained for ratios of boundary-layer thickness to cylinder 
radius between 0.249 and 0.654 ratios of cylinder to 
plate coefficients between 1.075 and 1.196. The au- 
thors find no unequivocal support for the existence 
of increments of that order by comparing their test 
results with flat plate values from other experimentators 
and suggest that in the case of the cylinder the higher 
“wiping capacity’ of the flow will produce a smaller 
boundary-layer thickness. 

The method proposed in the present paper leads to 
results that are intermediate but closer to the first one 
of the extreme cases. Instead of keeping the friction 
coefficient or the boundary-layer thickness constant, it 
is assumed that the relation between friction coefficient 
and boundary-layer thickness applied to the flat plate 
in referetice | remains applicable to the cylinder. The 
assumptions regarding velocity profile, static pressure, 
and total energy are also the same as made for the flat 
plate in that reference. 


(3) THICKNESS OF THE BOUNDARY LAYER ON A 
CYLINDER 


To simplify the calculations it is required that the 
direction of the free stream just ahead of the cylinder 


leading edge is parallel to the surface, that the boundary 
layer is turbulent from the beginning, and that it starts 
at the leading edge with zero thickness. 

Under actual conditions these requirements might 
be best approximated by equipping the cylinder with 
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a suitably shaped nosepiece and sucking off the bound- 
ary layer formed along the nose. 

If no pressure gradient has to be taken into account, 
the momentum theorem reads, for this case, 


Py =6§ 
d ‘ 
/ pu(U — u)2x(r + y) dy = 2nrr (1) 
dx Jy=0 

where x is the distance from the leading edge along the 
surface in the flow direction, y is the normal distance 
from the surface, and 6 is the boundary-layer thick- 

, y-lay c 

p and u are density and velocity inside, and U is 
r is the cyl- 


ness. 
the velocity outside, the boundary layer. 
inder radius, and r is the shear stress at the surface, 
All flow characteristics are time-averaged. 

The momentum thickness 0 is defined by 


! ci y 
6= if pu = w) (1 +*) ay (2) 
pul . y=0 r 7 


where py is the density outside the boundary layer. 
Introducing the ratio 0/6 = Q and splitting the in- 
tegral, the following expression is obtained: 


7] *] 
pu u y *f pu 
Q = 1 — 1 {- 
I = ( r)a(z) +5 0 oe 
(: = A y d (7) (2a) 
Us 6 5 


The first integral apparently represents the ratio of 
the momentum thickness to the full boundary-layer 
thickness for a flat plate. By applying the same pro- 
cedure as described in Part (III) of reference 1, it can 
be shown that the second integral represents half the 
value of the same ratio for a plate with a velocity pro- 
file of half the exponent. 

If the profile is denoted as 


u/U = (y/6)*"™ (3) 


this relation can be indicated by writing Eq. (2a) 


plate plate 
Q, = Q, + (6/2r) Qon (4) 
plate 


Numerical values of 2, as a function of Mach Number 
up to Ma = 10 are given in Fig. 3 of reference | with the 
parameter m varying between 5 and 10. The corre- 

plate 
sponding values of Q:, can also be determined from the 
data of Fig. 1 and Eq. (26) of that reference. (NOTE: 
the term (2 is designated by 0 in reference 1.) 

To eliminate the wall shear stress 7 from Eq. (1), use 
is made of Eq. (34) of reference 1, which is the above- 
mentioned relation derived from pipe-flow experiments. 
After replacement of the pipe radius d/2 by 64, this rela- 
tion becomes 


r br 1/4 i [(w+ 1)/4] 1 
~ = 0.0233 ( ) (= ) (5) 
pul Upvd Ty 


where yy is the free stream viscosity, 7'/7y is the tem- 
perature of the mean state across a pipe section related 
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TURBULENT BOUNDARY 
to temperature at the axis of that section evaluated 
for 1/7-power law velocity distribution, and w is the ex- 
ponent of the viscosity variation power law, in reference 
1 taken as 0.8. 

Using the identity 9 = 26 and Eqs. (2), (4), and 


(5), we obtain, from Eq. (1), 


d a 6° a Ki vi 
62, + Qe, |} = 0.0233 : x 
dx 2r l pu 
T —(3—w)/4 
(. ) gis (6) 
7 


The superscript “‘plate’’ on Q has been omitted for 
simplicity, since only flat plate values of this term appear 
in the following. Eq. (6) can be readily integrated 


to give 


1/4 ‘i —(3—w)/4 
0.023: (i ) (|) x (7) 
Upr 7 


or, if solved for 6/x and with U’pyx/py = Re, 


i (3—w)/5 
0.059 2, ~'“* (- ) Re~"/* 
U 


= eos (7a) 


mo 


For 6/r = 0, Eq. (7a) becomes identical with Eq. (37) 
of reference | for the plate. Dividing Eq. (7a) by the 
plate equation, the effect of convex lateral wall curva- 
ture upon the boundary-layer thickness is obtained. 


6 5 Q», 6\" (8) 

’ (1 + tg *) 

9 Q, 

Q,/Q, has been plotted versus Mach Number in Fig. 
lforn = 5,7, and 9. This ratio increases with Mach 
Number and decreases with m. Both effects are mod- 
For Ma = 0, Q2, /2, = (n + 2)/(2n + 1), which 
(S) thus reduces for the in- 


erate. 
for n = 7 yields 0.6. Eq. 


compressible case to 


(2) I 
On i (1+ (6/3r)]”* 


inc 


(Sa) 


The ratio 6/6, is plotted vs. 6/r in Fig. 2a with 
Mach Number as parameter using » = 7 which is con- 
sistent with the applied shear stress relation. For 
Ma = Oand 6/r = 0.1 the boundary-layer thickness is 
about 3 per cent less than on a flat plate under same 
1 the decrement is about 20 per 
effect 


conditions; for 6/r = 


cent. Compressibility slightly increases the 
of curvature. 

The fact that Eqs. (8) and (Sa) and their graphical 
representation in Fig. 2a are implicit with respect to 
6 involves no serious complication for practical pur- 
poses. First determine 6,, and use 6,)/r as independent 


variable to obtain a preliminary value for 6/6,, from 


LAYER 


to 
ou 


ALONG A CYLINDER 


Eq. (8) or Fig. 2a. With this preliminary value, cor- 
rect that initially used 6,,/r to obtain an improved 
value for 6/6,,. If 6/r > 0.1 it may be necessary to 
repeat the operation. 


(4) FRICTION COEFFICIENTS 
With 6/r known, the coefficients of local and mean 
friction can be found in the following way: 
The /oca/ friction coefficient is defined as 
C, = 21r/py? (9) 
With the aid of Eq. (5) and by relating the skin friction 
value for the cylinder to that for the plate, 


(10) 


Crp 


Combining Eq. (10) with Eq. (8) yields 


Cy (: ‘ 5 Qe, *" 
Crypt 92, 7 
Numerical values for Eq. (11) are plotted versus 6/r in 
Fig. 2b for Mach Numbers 0 and 10. 
The effect of curvature on c,; is small and becomes 
appreciable only if 6/7 approaches unity. 


is obtained. 


(11) 
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Characteristics of the turbulent boundary layer along 


Fic. 2 
a cylinder in compressible flow related to flat plate values vs. the 
ratio boundary-layer thickness to radius of the cylinder 1/7 


power law velocity distribution. Parameter: Mach Number 


(2a) Full thickness of the boundary layer; (2b) local friction 
coefficient; (2c) mean friction coefficient. 
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The mean value of the friction coefficient over the 


C= (1 yf cy, dx (12) 


From the integrated momentum Eq. (1) and with 
the aid of Eqs. (2), (2a), and (4), the ratio of cylinder 


distance x is 


friction to plate friction is 


26 6 
Q, + — Qe, 
Cy x 2r (13 
— ») 


Cipt (26, eile 
Using Eq. (8) to eliminate 6/6,, from Eq. (12) yields 


6 Qon 


Cy or Q, 


= = re? 14) 
CIpi 56 Qon 7 
aul Or <. 


The ratio ¢;/C;,, is plotted versus 6/r in Fig. 2c for 
Mach Numbers 0 and 10. The deviation from the 
plate value is still smaller than for the local friction 
coefficient, which is obvious from the fact that the mean 
friction coefficient includes all values from the leading 
edge where 6/r, and thus the curvature effect, is smaller 
than at the downstream end position. Eq. (14) and 
Fig. 2c also represent the ratio of momentum thickness 
on the cylinder to that on the flat plate, as is evident 
from Eqs. (4) and (13). 

If in Eq. (13) is set 6 = 6,, and Q4/Q2; = 0.6, then 


Cy/Cpp, = 1 + (3/10) (8/r) (15) 


which is the result obtained by Jakob and Dow.*® 


(5) DISPLACEMENT THICKNESS 


In determining the influence of the cylinder radius 
upon the displacement thickness it must be remembered 
that this quantity for higher Mach Numbers becomes 
comparable to the full thickness of the boundary layer 
itself, and it is therefore not permissible to disregard the 
quadratic term as in the case of the momentum thick- 
ness. 

From the continuity of mass flow, it follows that 


puUnl(r + 6 — 6*)? r?] = 


wy=65 
/ pu2r(r + y) dy (16) 
y=0 


where 6* denotes the amount of radial displacement a 
free-stream line experiences up to its intersection with 
the “edge’’ of the boundary layer. This definition 
differs from that used conventionally for cylindrical 
surfaces. It is felt, however, to be preferable because 
it represents the correct influence upon the potential 


flow. 
After rearrangement and division by 2mrépyl/, Eq. 
(16) can be written in nondimensional form 
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6 g* 5* 5* 


2r 26r 6 r 
1 1 
pu y 6 pu y y 
f (5) + f a (%) (17) 
0 pul 6 rJo0 py 6 6 
In the same way as for the momentum thickness, it can 
be shown for the displacement thickness by the 


procedure in Part (III) of reference 1 that as 


plate 


1 
u ) 6* 
[2 (2) =1-( ) (18) 
0 pul 6 5/n 


there is also 


1 plate 
1 ! ! l 6* 
f ~(?) d (?) = [ - ( ) | (19) 
0 pul \6 6 2 Jon 


where again the superscript denotes plate values and the 
subscript velocity profile parameters as defined by Eq. 
(3). With Eqs. (18) and (19) and A = 6*/6, Eq. (17) 
becomes, after rearrangement and canceling out equal 
terms on both sides, 


) 


6 6 plate 6 plate 
A, (1 + ) — A,’ =A, + Aon (20) 
r 2r 2r 


or 
plate 
14 6 a) 
A, aay J : 
plate = plate (21 
An 6 Ay, Ay 
Tx ik ae 
r 9 plate 
By 
plate 


The ratio (A:,/A,) which is required for evaluation 
of Eq. (21) is plotted in Fig. 3 versus Mach Number 


plate 
forn = 5,7, and 9. As, has been determined with the 
aid of Eq. (24) and Fig. 1 of reference 1. The be- 
plate plate 


havior of (A2,/A,) is similar to that of (Qe,/Q,). For 
Ma = 0, 


plate 
(Ao, /A,) = (nm + 1)/(2n + 1) 
plate 
For values of 6/r up to 0.1, the ratio A,,/A,, can be calcu- 
lated conveniently and with an error of less than 0.5 
per cent in one step from Eq. (21) by setting on the 


plate 
right side A,/A, = 1. For larger values of 6/r, one 
has to use either successive approximation or the 


solution of the quadratic Eq. (20), which is 


~ 
~ 


) plate §2 plate 
jl + (1 — A,) + (1 — Ao,) (22) 
\ , - on) | ( 
The ratio displacement thickness on the cylinder to 
displacement thickness on the plate is 
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TURBULENT BOUNDARY LAYER ALONG 


plate 
6*, 5,,* ~~ (A, A,) (6 Opi) (23) 
or with Eqs. (8) and (22) 
6 | 26 plate 62 plate 
T si i] + Q — A.) + (1 = Ao, ) 
"dad i stelle 
s * 2 plate 
Opi gPlate 5 5/Q, 4 
By) ir ( 
r 9 r\Q, 
(24) 
§*/5,,* is plotted versus 6/r for n = 7 and Mach Num- 


bers 0, 2, 5, and 10 in Fig. 4. As can be seen, the in- 
fluence of lateral curvature on the displacement thick- 
ness is larger than on any of the other quantities 
considered. For incompressible flow, about 90 per 
cent of the flat plate value are obtained at 6/r = 0.1 
and only about 50 per cent at 6/r = 1. Itis noteworthy 
that in this case Mach Number somewhat decreases 


the effect of curvature. 


(6) SUPPLEMENT 


In the absence of experimental data for comparison 
with the results above, some considerations are re- 
quired as to the reliability of the basic assumptions and 
the effects of modifications to which they may be sub- 
jected. 

Measurements of turbulent boundary layers in com- 
pressible flow are at the present time available only for 
For flat plates the velocity 
profiles obtained from Pitot tube readings by Wilson, 
Young, and Thompson® at Mach Numbers from 1.6 to 
2.2 and Reynolds Numbers from 4 X 10® to 20 & 106 

Except 


two-dimensional cases. 


are in good agreement with the |, 7-power law. 
for a region in proximity of the wall, the same can be 
said about similar tests by Rubesin, Maydew, and 
Varga’ at Mach Number 2.5 and Reynolds Numbers 
varying from 2.7 X 10® to 6.2 & 10° Comparison of 
the values for 6*/x and cy, derived from these experi- 
mental data with the theory of reference 1 indicates 
their variation with Re~’’*, while the Mach Number 
effect appears to be underestimated by that theory. 
This does not, however, influence to a great extent the 
results of the present paper, since its purpose is mainly 
to show differences in the behavior between flow along 
cylinders and flat plates where only a second-order 
effect of Mach Number comes into play. 

At higher Reynolds Numbers it must be expected 
that, similarly as for incompressible flow, the 1/7- 
power law for the velocity profile and the Re~'”* law 
Yor the friction coefficient will no longer hold. The 
velocity profile will probably change to a fuller shape, 
and the effect of Reynolds Number on skin friction will 
become weaker. Both effects are interrelated. If -1/m 
denotes the Reynolds Number exponent, the relation 
between m and n can be obtained approximately from 


a dimensional analysis by von Karman." It reads 
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m = (n + 3)/2, which for n = 9,t for instance, pre- 
. - + . 1/, 
dicts an effect of Reynolds Number according to Re™ “*. 

Eq. (8), which in the general form reads 


6 1+ m Qo, *) (m—1)/m si 
= ) 
6 2m — 1, 4 


pl 


t+ Wilson® obtained from measurements at the wind-tunnel 
wall for Mach Number 2 a velocity profile close to a power law 
with » = 9. The equivalent flat plate Reynolds Number was 


determined to 50 106 
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Fic. 3. Nondimensional displacement thickness of second- 
order velocity profile related to same value of first-order profile 
vs. free-stream Mach Number for various power law exponents of 
the first-order profile. Values refer to flat plate flow 
































Displacement thickness of the turbulent boundary 
layer along a cylinder in compressible flow related to flat plate 
value vs. the ratio boundary-layer thickness to radius of the 


Fic. 4 


cylinder. 1/7-power law velocity distribution. Parameter: 


Mach Number 
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can be written with the relation above 


6 n +4 2 Qon 6 —(n+1)/(n+ 3) 
= ( + ) (Sc) 


Spt 2n+4Q, r 


For large values of n, Eq. (Sc) converges toward 
1 Qo, 6\~" 
~{1+ (n— o) (Sd) 
: r 


which for the incompressible case reduces to 


Re ae 8 
= n—-> © Sse) 
5p ine 4r ' , 


For 6/r = 1, the value obtained with Eq. (Se) differs 
by less than | per cent from the one obtained with Eq. 
(Sa) where m = 7. With respect to the small depend- 
ence of the ratio Q.,/2, upon Mach Number, it can be 
concluded that Eq. (8d) will give results that are simi- 
larly close to the values obtained with Eq. (8) for com- 
pressible flow. This means that for nm > 7 and 6/r < 1 
the ratio 6/6,, is practically independent of n. 

The ratios of the local and mean friction coefficients 
decrease in steps of about 0.5 per cent if ” is increased 


Spi 


in steps of unity between 5 and 9. Forn—~> o, 


C7/ Ci > €/ Cr —> | 


These results are valid for the assumption that the ve- 
locity distribution is not influenced by lateral wall curva- 
ture. There is reason to expect, however, that the 
profile will be altered by convex curvature toward 
larger n—that is, to a fuller shape compared to the flat 
plate—because the flux of momentum from the free 
stream toward the wall is concentric. This will lead 
to higher velocities within the boundary layer and a 
stronger velocity gradient near the wall. 

An effect in this direction has been noticed for in- 
compressible flow by Schultz-Grunow,’ who found that 
velocity profiles are fuller on flat plates than in pipes. ft 

In order to allow an estimation of the effect of differ- 
_ent velocity profiles, some calculations have been car- 
ried out assuming » = 9 and an accordingly altered 
shear stress relation for the cylinder while retaining 
n = 7 for the plate. For this condition the ratios of 
cylinder to plate values depend, besides on 6/r, directly 
on the Reynolds Number, and it is found that for Re > 
10’ the iriction coefficients may well exceed those on 
flat plates by 20 per cent. It can be supposed that 
the effect of curvature on the velocity profile will be- 
come significant only if 6 > r, but the indication is 
given that the results of this paper might be low. 

The assumptions of constant total energy and con- 
stant static pressure have been checked experimentally 
for two-dimensional flow at Mach Number 2.8 by Spi- 
vack."” Variations of these quantities with wall dis 

1 The fact that this difference is rather small in view of the 
differences in the flow geometry may be due to the pressure drop 
inherent to subsonic pipe flow which tends to compensate the 
effect of the concave wall curvature. 
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tance by about 5 per cent have been noticed. The in- 
fluence upon integral parameters that enter the present 
analysis is irrelevant. 

It may be added that for a laminar boundary layer 
and incompressible flow the problem has been treated 
more fundamentally by starting with the differential 
form of the boundary-layer equation in a paper by 
Seban and Bond.!! Results given for the local friction 
coefficient are similar to those of the present paper but 
indicate a stronger effect of lateral curvature. For the 
behavior of the displacement thickness, no clear result 
has been obtained because of the limited number of 
terms retained in the series expansion. 


(7) CONCLUDING REMARKS 


The results obtained by the approximation method 
above indicate that the growth of the turbulent bound- 
ary layer along a cylinder becomes considerably slower 
compared to the flat plate if the boundary-layer thick- 
ness approaches the dimension of the cylinder radius. 
The effect is still more pronounced on the displacement 
thickness 

More important for engineering purposes, however, 
might be that the influence upon the mean friction co- 
efficient comes out to be extremely small and hardly 
exceeds the uncertainty of present-time experimental 
determinations even if the boundary-layer thickness 
equals the cylinder radius. This would substantiate 
the use of flat plate values for cylinders and vice versa. 
The reliability of this conclusion depends, of course, 
upon that of the assumptions involved, mainly the in 
variance of velocity distribution and shear stress law 
with curvature radius. 
that the results will, by their presentation as fractions 
of corresponding plate values, not be sensitive to 


It can be expected, however, 


variation of parameters that simultaneously affect 
conditions at both cylinder and plate. This is evident 
from the small influence of Mach Number. 

If at higher Mach Numbers or higher Reynolds Num 
bers the actual behavior of the turbulent boundary layer 
on the plate does not follow the relations of reference 
1, then a similar deviation differing only by higher 
order effects will also exist for the cylinder and qualli- 
tatively useful results can be expected by basing the re- 
sults of this paper on improved plate values as they may 
be given, for instance, by the theories developed in 
references 7 or 8. 

Once dependable relations for the case of the flat 
plate are established, it will be desirable to check also 
in particular the effect of lateral curvature upon the 
boundary layer by appropriate experiments. 
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On the Flight Dynamics of Slender 
Special-Purpose Aircraft 


R. M. ROSENBERG* ann GEORGE STONERt 
Boeing Airplane Company 


SUMMARY 


In this paper, general equations of motion are derived for 
slender special-purpose aircraft of nearly circular cross section. 
It is shown that considerations of steady roll and lack of circular 
symmetry of cross section lead to Mathieu equations. Aero- 
dynamic force derivatives for flight not near the flutter speed are 
given which appear to be useful in that their application has re- 
sulted in theoretical predictions that agree satisfactorily with 
flight experience. Finally, some possible instabilities are dis- 
cussed. The method is illustrated by an example 


INTRODUCTION 


a THE STANDPOINT of the flight dynamics of 
modern special-purpose aircraft such as missiles, 
it is of interest to develop fairly general equations gov- 
erning the flight of a slender, essentially cylindrical 
The 
solution of these equations leads to the analysis of os- 
cillations that may occur in flight and which may be- 
come unstable under certain unfavorable circumstances. 
The results of the theory presented here have been com- 
pared with flight experience, and astonishingly good 


agreement was achieved. 


body that is elastic in a sense to be defined later. 


THE DEGREES OF FREEDOM 


The aircraft under consideration is assumed to be 
slender and to possess geometric, elastic, and inertial 
characteristics having nearly circular symmetry with 
respect to all points on a longitudinal axis through the 
body of the aircraft. 

The right-handed orthogonal xyz-system of coor- 
dinates is fixed in the craft; it coincides with the prin- 
cipal axes, and the x-axis is positive in a forward direc- 
tion. Inasmuch as the coordinate system is right- 
handed, the positive direction of all other quantities 
relating to the coordinate system are now defined. 

The nomenclature and sign convention of the rigid- 
body degrees of freedom and of certain kinematical and 
dynamic quantities relating to them are defined in 
Table 1. All velocities, forces, and moments are posi- 
tive in the direction of positive displacements. 

In addition to the degrees of freedom defined in Table 
1, the aircraft is considered to be elastic in the sense 
that it can execute beam bending vibrations. <A vibra- 
tion in the xy-plane leads to a displacement o(x, /); that 
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in the xz-plane leads to a displacement r(x, ¢). Both 
vibrations are restricted to occur in the fundamental 
modes only, and no elastic deformations, except for the 
two listed, are considered. The nomenclature regarding 
the kinematical and dynamic quantities referring to the 
elastic degrees of freedom are listed in Table 2 where 
d/dt = D. 

o is positive in the direction of positive y. 
All other directional 


rT is pos- 
itive in the direction of positive z. 
quantities in Table 2 are positive when they act in the 
direction of positive displacements o and r. 

The aircraft under consideration is assumed to have 
certain controls having principal actions as listed in 
Table 3. These are described with reference to an axis 
system fixed in the aircraft 
by an observer on the aircraft. 
instance, may give rise to a displacement that would ap- 
pear vertical to an observer on the ground if the appli- 
cation of that force had been precededby a roll displace- 
ment of 90°. The symbols in the right-hand column 
are regarded as indicating the degree of application of 
the controls leading to the responses indicated in the 
All control actions are 


i.e., they are those seen 
Thus, a side force, for 


left-hand column of Table 3. 
positive in that direction in which they give rise to re- 


sponses in a positive direction. 


THE BENDING VIBRATIONS 


In an aircraft capable of roll, the inclusion of the 
vibrational degrees of freedom presents considerable 
difficulties. 

Consider that the aircraft is not rolling and that 
bending vibrations are set up in the xz-plane, say. If 
the aircraft is now given‘a displacement in roll, the 
plane of the bending vibration does not rotate with the 
aircraft; Let the 
plane of bending vibrations be the sx-plane, and, before 
rolling, the s-axis coincides with the z-axis. When a 
rolling displacement ¢g(f) takes place, the s-axis re- 
mains fixed in space, the g-axis remains fixed in the air- 
The vi- 


it tends to remain fixed in space.! 


craft, and the angle between the two is ¢(f). 
brational displacement 7 now has a y-component equal 
to r(x, t) sin g(t) and a z-component equal to r(x, ¢) 
cos g(t). Obviously, these displacements introduce 
nonlinearities into the equations of motion which render 
analytical treatment extremely difficult. If the roll is 
steady at the angular velocity fo, the components be- 


come r(x, /) sin pot and r(x, t) cos pot. In this simpli- 
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TABLE | 
Coordinate -——— Parallel to Coordinate Axes——-—. - ——_——_ ———About Coordinate Axes - = 
Axes Velocity Acceleration Forces Displacement Velocity Acceleration Moment Inertias 
x u az X rr p deg i A 
y v dy r 6 q ae M B 
z w dz A y¥ r ay N c 


fied case, the vibrational displacements introduce 
periodic coefficients into the equations of motion which 
also make the analysis difficult. Finally, if the air- 
craft has unequal stiffness in mutually perpendicular 
directions of beam bending, then the spring restraint 
becomes, even in steady roll, a periodic function of time. 
Again, the system of equations will have periodic co- 
efficients leading to Mathieu equations. Because of 
the difficulties resulting from a consideration of roll, 
the analysis will be carried out for a nonrolling aircraft, 
or g(t) = p(t) =0. 

When no roll is permitted, the vibrational displace- 


ments are 


a(x, t) = oo(t) f(x) ) 


I= (1) 
tT(c,t) = To(t) fo(x) f 
where /;2(x) are the mode shapes satisfying 
fio(x) dm = 0 (2) 


l 


where dm is an elemental mass located at x and the 
symbol ff indicates integration along the entire 
length of the aircraft. 


PROCEDURE OF ANALYSIS 


The general procedure of analysis will be to write 
one equation of motion in each degree of freedom. 
Each equation will be of the form a = f/m, where a 
is an acceleration, linear or angular, f is a force or 
torque in the direction of the acceleration, and m is a 
mass or inertia. 

Because of the essentially aerodynamic nature of the 
problem in hand, the natural dependent variables are 
the velocities rather than the displacements. 
ception to this occurs in the vibrational degrees of free- 


An ex- 


dom because vibrational displacements, as well as ve- 


TABLE 2 
Natural 
Displacement Velocity Acceleration Forces Frequencies 
o o = Do ¢ = De ad ray 
T += Dr 7 = D*r Q ws 
TABLE 3 
Principal Action Symbol 
Thrust é 
Side force n 
Vertical force ¢ 
Roll torque p 
Pitch torque 8 
Yaw torque Y 


locities, produce changes in angle of attack. Thus, as 
far as the vibrational degrees of freedom are concerned, 
the displacements will be regarded as the dependent 
variables. 

Inasmuch as the coordinates move in space, the ac- 
celerations in the equations of motion are the absolute 
accelerations of a moving coordinate system. These 
will, however, differ from the well-known forms of 
these accelerations in that they contain contributions 
from the vibrational degrees of freedom. The forces 
(or torques) in the equations of motion are the summa- 
They 
arise from two sources: One type of force is that arising 
This type of force 


tions of all forces that are not inertial in origin. 


from the motion of the aircraft. 
consists of such quantities as lift, drag, pitching mo- 
ments, aerodynamic damping, and others, and these 
will be regarded as response forces; thus, they will 
appear on the left-hand side of the equations of motion. 
The second contribution is introduced by activation 
of the controls. The forces arising from this source 
will be regarded as exciting forces and will appear on the 
right-hand side of the equations. 

The solutions of the equations of motion give the 
time histories of the motions; thus, they must clearly 
depend on the manner of activation of the controls that 
cause the motion. The equations here derived are 
completely general in this respect, and the investigator 
may search for steady-state solutions by introducing 
periodic control activation, or he may obtain transient 
solutions to unit-step or arbitrary control activation. 
The theory here presented was actually applied to 
transient solutions resulting from unit-step control 
motion. 


THE ACCELERATIONS 


Derivation of the linear accelerations presents no 
great difficulty. They are easily shown to be, respec- 
tively (for p(t) =0), 

a, = (du/dt) + (w+r)q — (vt o)r) 
(dv/dt) + (da/dt) + ur (3 
a, = (dw/dt) + (dr/dt) — uq ) 


ag 
| 


Eqs. (3) are obtainable from purely kinematical con- 
siderations. 

The angular accelerations must be derived from dy- 
namic relations. They are obtained from Euler’s equa- 
tions of inertia moments of a body rotating about a fixed 
point coinciding with the center. of gravity, and their 
simple form used here requires that the coordinate axes 
coincide with the principal axes. 
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FLIGHT 


The derivation of the angular acceleration about the 
y-axis will be shown here. For the other accelerations, 
only the results are listed. 

The angular velocity vector @, whose components 
are the angular velocities about the xy and 2-axes, re- 
spectively, is 

@®=qj+rk 
The velocity of a point P located at 


F=xi+ yj + zk 


v7, =(2+r)¢g—(yto)r 
Vy = XP 
v7, = —xq 
Then, the angular momentum of a mass dm located at 
Pis 
= 7dmi 


= hii + Ayj + Ak 


> 


where 
Az = dm [(y + a)v, — (2 + 7)2,] 
hy = dm [(z + 1r)v, — xz] 
A, = dm [xv, — (y + a) 2;]} 
The inertia moment due to the mass dm about the x- 
axis 1S 
dL = (dd,/dt) + q rz — Wry 
dm} —xoq — x(y + a)g — xtr — x(a + r)* — 
g¢+roytog?t+ (@+7) yt a)r? + 
[(y + o)? — (2 + 7)* Ir} 


The total moment is 


L= -4 f x0 + «) dm — if x(z + r) dm — 
ef +17) (y+) dm + nfo +o) (2+ 7r)dm+ 
o| fo + a)?dm - fe + 7) dm |- 
qo f xfi(x) x dx — rto f xfo(x) x dx 


By making use of the following notations 


y) dm = m; dm 
1 
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where x = dm/dx is the mass per unit length and use 
has been made of Eqs. (1). It can be shown that, if the 
vibrational displacements o and 7 are sufficiently small, 
the above expression for the moment can be approxi- 


mated by 


L = (C — B)qr - ain f xfi(x)x dx — rv f xfo(x)x dx 
l l 


where higher terms of the small quantities and products 
of inertia with regard to the vibrations have been neg- 
lected. Finally, the acceleration about the x-axis is 
obtained by dividing the expression for the moment 
by the moment of inertia of the body about the x-axis, 


or 


C-—B l : 
a. = gr — - xfi(x)x dx | goo — 
¢ , 2 { f 1 
f Zé l 
l 
xfo(x)k dx | rt) (4) 
Asi ~ 


Because the rolling velocity vanishes identically, par 
ticularly simple expressions for the other two angular 
accelerations are obtained in an exactly analogous man- 


ner. They are 


ay = dq/dt\ 
a, = dr dt{ 


THE EQuaTIONS OF MOTION 


Assuming that the force per unit length X, parallel 
to the x-axis, is a linear function of all variables and of 
all controls, the equation of motion of an elemental 
length of the aircraft (of length dx and mass dm) may 
be written as 


(% — rv + qw — roofi + Qtof2) dm = 


Dae +O0j+ Bk + Zah+ 
‘ J 


X ,Tof2 + X tof ; dx (6) 
where 7 can take on the values u, v, w, g, 7; j can take 
on the values £, , ¢, p, 8, y; the subscripts denote par- 
tial differentiation; and the bar denotes forces per unit 


length. 


X,/m = X, 


Kk dx; 


y) X,dx = X;; d/dt = D 
l 


and of Eqs. (2), the integral of Eq. (6) over the aircraft length is 
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D(u) — rv + qu — 1D Xt + a0 fre dx + Dio) f Xf dx + nf X fo dx + p(n) f X,fe dx} Xj 
i l I Y j 


(7) 
Similar equations for the forces parallel to the y- and z-axes are found to be 
J a oe ra 
D(v) + ur —\i1 it + uf Y,f: dx + Do) f a1 fi dx + nf, fo dx + D(r0) f/f Y, fe dct = 2Y,j 
& « l 
(8) 
D(w) — uq — ) Zit + oo f Zofi dx + D(a) f Z3hidx + nf Z,f2dx + Dim f 2, fe dec} = > 2,j 
; l eee ‘il 
(9) 


The equations of torques about the coordinate axes are obtained similarly to Eqs. (7) to (9). Of these, the roll 
equation does not exist because roll is not an admissible motion.* 
With the aid of the following notation: 


M,/B = M,, | M,dx = M; 
! 


Ni, / N,dx = N; 
the equation of rotation about the y-axis is 


D(q) - 1D Mya + ov fa. dx + Dio) f M, fidx + wf M,fsdx + Dow f, M, fe axl = > M,j 
i l 


(10) 


N, 1C 


and the equation of rotation about the z-axis is 


D(r) - 1z Nia + n f.X ofidx + Do) f N,fidx + nf XI 2 dx + D( wf Ni fo dx | = WN, 4 
l l 


(11) 


The remaining equations to be derived are those in the vibrational degrees of freedom. These are most easily ob- 
tained by application of Lagrange’s equation. That for the o degree of freedom is 


D?(a0) + w1200 — )Di [rs dx + of Pf? dx + D(a) / Pi fi dx + vf P fife dx + D(r) X 
s « . l 
f Pas fa dt rif? jfidx (12) 
l 


where use has been made of the notation 


II 


P, 


P, if fiex dx 
f 

O, | fo*x dx = QO; 
1 


the equation in the 7 degree of freedom becomes 


D*( To) + Wwe — 1xi f O: fs dx of wf 0, fife dx + Do) f Osh dx + TO fous dx fe D(70) 4 
i l 
ff ote dxt =D f O,feav (13) 


* It is not an entirely trivial matter to show the nonexistence of the roll equation because Eq. (4) shows that a roll acceleration 


With the notation 


does exist even when p(/) 0, and there are undoubtedly nonvanishing forces due to many of the variables that contribute to the 
roll torque. 





V 





Xj 


“J 


(9) 


» roll 


10) 
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The assumption that all forces are linear in their arguments is less satisfactory for the force X than for the others. 
Therefore, Eq. (7) should be more properly written as 


Diu) — rv + qu — +34 | X; di +f / XY, fidoo dx +f / NX, fidoy dx +f [Xsan dx + 
1 l ve 
IS NX, fodto dxt = oy | > g dj 
l j 


This equation is seen to introduce nonlinear terms into the system of equations of motion. However, the diffi- 
culty is circumvented by the following well-known simplification : 

In the conventional analysis of airplane dynamic stability, it is a familiar assumption to consider the variation in 
forward velocity negligible compared to the mean value of that velocity, or u = U' + u’, where u’ < U and U is 
aconstant. In the problem here considered, it is sufficient to regard the maneuvers so small that their effect on 
the forward velocity need not be considered. For this special case, the first of the equations of motion—1.e., Eq. 
(7)—must necessarily express the equilibrium between thrust and drag. Thus, Eq. (7) simplifies to 


Kat fe X ifo (1 4) 


where lL’ is the constant forward velocity and £ may be regarded as the constant throttle position. Eq. (14) is an 
equation containing none of the variables under consideration; thus it separates from the system of equations. 
Those that remain are conveniently written as the matrix equation, 


[E] {V} = [Cl] {x} (15) 


(— f rnas)(— fnac)(— raids) (— Pas) 


(-foe)(-foana)(- fsa) 


0. 
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where 
if) = 
- | Y,fidx D / Y, fe dx D 
(D — Y,) (— Yy) (—¥,) (U — Y,) i fo 
= | Tf, dx “ | Y.fs dx 
l 
- | Z,/, dx D - | Z, 2 dx D 
(—Z,) (D — 2.) —(U+2Z,) (-Z,) : ,, 
- [2 fi dx - [2s dx 
l l 
-f Mh dx D - | M,f fo ¢ dx D 
(—M,) (—M,) (D—M,) (—M,) > A 
- | M, fi dx -/ Mf dx 
l l 
= | N,fidx D\[ — / N, fo dx D 
(jt (—N,) (—N,) (D—N,) ‘> A 
- fn chi dx -f V fe dx 
l l 
“fe | 
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y; Y Y, - Y, 


n $ ’ } 7 
z Z, a £. , a Z, 
M, M, M, M, Mz M, 
N, N, N, N, Ng N, 


[Cc] = ; ‘ , . a 
| Pf; dx Pi fi dx / Pf dx | Phi dx / P afi dx 4 P, fi dx 
l l l l l l 
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fo f, dx | QO, fe dx | O, fo dx fos dx f Ooh dx fo. fy dx 
ae ee l l l l J 
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THE AERODYNAMIC DERIVATIVES 


Eq. (15) and the equations defining the matrices of 
Eq. (15) demonstrate the necessity for determining a 
large number of aerodynamic derivatives. 

It is realized that, in the problem at hand, the aero- 
dynamic forces (whose derivatives with respect to the 
dependent variables are desired) are those arising from 
unsteady flow about airfoils. Inasmuch as these forces 
have been determined for both subsonic incompressible 
and supersonic flows, their substitution in the equa- 
tions of motion presents no insurmountable obstacle. 
However, their use complicates the problem so that a 
solution without high-speed computing devices seems 
forbidding. 

It is not unrealistic to suppose that the application 
of the air forces arising from steady flow will predict 
the dynamic behavior under investigation reasonably 
well, provided the velocity of translation of the airplane 
is not near the flutter speed. However, in applying the 
steady-flow air forces, the investigator deprives himself 
of the possibility of analyzing flutter. 

The analysis here presented was undertaken in an 
effort to duplicate analytically, in character as well as 
quantitatively, the appearance of oscillations observed 
in flight. As a first step, the forces appearing in the 
equations of motion were considered to be those of 
steady flow—i.e., the effect of vortices shed into the 
wake was neglected. The results of this simplified 
analysis agree well with observation; therefore, it was 
considered satisfactory to treat the air forces in this 
simplified manner. 

It will be noticed that many aerodynamic derivatives 
appear simply as partial derivatives of forces with re- 
spect to some of the dependent variables, while others 
are integrals in which the derivatives occur in the inte- 
grand. The former are considered here to be quanti- 
ties known either from theory or from wind-tunnel 
tests. The latter will be termed integral derivatives 
for convenience, and it is these that will be considered 
in this section. Specifically, the integral derivatives 
will be here computed from the “‘simple’’ derivatives. 

It is supposed that the simple derivatives (which are 
known) were obtained by superposition of contributions 
arising from certain components of the aircraft. Thus, 
to consider one of these as an example, it is assumed that 
Y, was obtained by the operation 
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y * 1 0Ff, 


k 

, OY, 

i Ee de = 1, Bin sw te 
n=1 e 


n=1m OV 
(16) 


where the k numbers » represent k stations along the 
x-axis of the aircraft where components (such as air- 
foils) are located, each giving rise to an aerodynamic 
contribution to Y,. If downwash or any other factor 
should introduce a phase lag between the k components 
to be added in the manner of Eq. (16), this effect can be 
taken into consideration by considering some of the 
OY,,/0v as complex quantities. In any case, it is as- 
sumed that all components of the simple derivatives 
and thus their sums are known. It will now be shown 
how the integral derivatives can be obtained from the 
“simple’’ derivatives. 

The integral derivatives will be considered in certain 
groups, the first of these being that in which the dif- 
ferentiation is carried out with respect to a vibrational 
velocity. 

Consider P, which is the rate of change of the vibra- 
tional force per unit length P with respect to the vibra- 
Both the force and the velocity are 
in the y-direction. Thus, with regard to an element of 
length of the aircraft, P, is indistinguishable from Y,, or 


tional velocity a. 


oY/dv = 0P/d«a 
But, since 


oF 1o0¥Y oP l oP 


ov md’ Os Sfx dx de 


it follows that 
oP/de = (m Uf dx) (OY /dv) 


In view of Eq. (16), it is clear that 


ga m oY, 
P,; fer dx = - - | fi(x,) |? - (17) 
I ; Si fit dx u ov ~ | 
where /i(x,) is the value of f(x) when x = x,—.e., it is 


the deflection in the f;(x)-mode at that value of x at 
which the mth station is located. In a similar manner, 
all other integral derivatives with respect to linear veloc- 
ities can be found. 

In determining integral derivatives with respect to 
angular velocities, use is made of the operations 
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ra) _ O Ov O ra) Ow 
dr wor dW dwg 


However, since the right-handed coordinate system 
implies v = rx and w = —qx, it follows that 


0/Or = x(0/dv) 
while 
0/0g = —x(0/dw) 


Thus, for instance, 


. m oY, 
Pf,dx = - X_ fi (Xn) (18) 
I - Sifi2x dx u Ov in 


But it is found that 


- m or, 
Fe Xnfilx,) (19) 
| qj 1 Shir dx u Den Ih 


The next group of integral derivatives to be considered 
is that of nonvibratory forces with respect to the vibra- 
tional velocities. 

In view of the fact that, for example, 


oO} l oY 0} 
= dx = dx 

Ov m v Ov 

oO} l OY oO} 
— = — dx = — dx 

Oa m J! Oa l Oo 


it follows, from reasoning similar to that applied pre- 


eb oY, . 
J Y,fi(x) dx = > F(Xn) (20) 
! 


. oOo 


viously, that 


It remains to evaluate the integral derivatives of forces 
and moments with respect to the vibrational displace- 
ments. The aerodynamic effect of a vibrational dis- 
placement on a component of the aircraft is the change 
in angle of attack which the component experiences. 

Consider as an example the quantity 0/00. This 
quantity will be first regarded in its relation to the yaw 
angle y. In fact, 


OV /de = (OF /dp) (Op /da) 
But the yaw angle may be written as 


da d p - fi’ (x) 
y = => (oof 1(X)) = dol (x) = oof (x) : = 
dx dx , Si(x) 
fi’ (x) 
eS 
h(x) 
where the prime denotes differentiation with respect to 


x. Then 


OV /00 = (OV /dy) [fi'(x)/fil(x)] 


However, y is expressible in terms of the sideways ve- 


locity v. In fact, 


tany =v7/\UDSy 
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where the approximation was made in view of the small 
maneuvers considered here. Then, 
(OY /dv) (Ov/OW) 


wv/dy =U, dY/dW= 


or, finally, 


OV /de0 = (OF/dv) U[fi'(x)/filx)] (21) 


Similar reasoning, together with the fact that w/l’ = 
— a, leads to 


OV/dOr = —(OV/dw) Ufo’ 


(x) /fo(x) | (22) 
From Eqs. (21) and (22), two typical integral deriva- 
tives with respect to vibrational displacements are 
evaluated by methods similar to those employed for 
Eqs. (17), (18), (19), and (20). They are 


. 


| Pf dx 
I 


Ps oF. . 
/ V.fedx = —U> fo! (Xn) (24) 
l 


. ow 


ll 


—— 
l u ap (x, ) (23) 


and 


A list of all required integral derivatives is given in 
Table 4. Their derivation is in all cases analogous to 
one of those shown in detail in this section. 
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formvarwe | __evawwavion ___Demivarive | Evauuation | 
| Seeex (n/['« xdz) 338 fee.) ie f,dx Listes 
| [Retde | (m/ffrndx)? Bt eax, | 1 f.dx | Ree fin 
| [Stax | (m/fer ds) 3 f crafaces | fp Piefde | ae 
| [@effax (n/[the ae) 285 stay | [Rede Titi 
| (tax | e/fe ca)S hey | [Agdx | Liem | 





| (n/[ thea) ete, |[Refax UR Fe feo 

ate) | [tax Uae he 
[zfax UL tie 

| [Exf.de une fix | 

[Bede | Ouleeady Fetmaties 
[@stdx | (w/ (feds) S356) |[Rtede | {rou /f @ cine fx fteo | 

| [&f,dx | ‘ogess gent ,) fer f,f,dx (muff ‘\dzL ov abel x rficx, 
[agar | eM nav nfs [ace (moult cay Bhertien 


| SRefax | 
ff fidx L(n/[*«d) eda) 2 
[Bde | yes satthoans 
[fax t/fcee) 5t fix) 
JQfidx (m/f *«dx)> See Serf (x,) 





| [Shee | aC eda) | [Pi fax Us Fees 
| Sifax | ERE [Pi gae us te ices 
| [Wide | 25 Betas J Rifdx UZFe tees 
te F | 
| J2,fidx Trt |[Rufedx | UL ae 0%) | 





JOURNAL 


OF 


THE 





AERONAUTICAL 





TIME SGALE 
1 2 3 4 5 6 




















™) 








ACCELERATION SCALE 


c 





























TIME SCALE 
3 4 


~ 




















ACCELERATION SCALE 
o 























Upper curve is the result of flight measurements for 
Isower curve is the result predicted by 


Fic. 1. 
one flight condition. 
theory. 


A SIMPLE APPLICATION 


The foregoing development will be here applied to 
the following example: 


A step-function type of control activation in ¢ will 
be applied at ¢ = 4; all other controls will be left un 
disturbed. Thus, the controls activity considered will 
be ¢ = lf, or 1¢(¢t) = Ofort < # and 1¢(¢#) = 1 fort bp. 
Furthermore, only “‘primary’’ aerodynamic forces will 
be considered to exist. Thus, the following derivatives 
vanish—.e., 


Nic Fa Va BV ed Ki ite te. Vs Be Ze Z., 
Si Zhi dx, Si Zs fr dx, Si Pu fy dx, Si Pf, dx, 
Si Pi fife dx, Ni Pifife dx, Si P-fy dx, Si On fe dx, 
Si Orfedx, S) Oufife dx, Si Os fife dx, St M, fr dx, 
Si Mfr dx, M., M,, Si N,fo dx, Si Ni fo dx, 

Weis IN gig LY 


t 


are all zero. In an aircraft of conventional design, 
these derivatives are actually sufficiently small com- 


pared to those retained that they may be neglected com- 


pared with all others. 


fq. (15) becomes 


For this system, I 
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Clearly, the equations giving rise to the first three rows of Eq. (25) are functions of those variables of which the 
equations corresponding to the last three are not. This observation is equivalent to noting that in a nonrolling 
aircraft the longitudinal and lateral motions may be discussed separately. 

Inasmuch as the proposed control activity is one giving rise to longitudinal motion, the lateral motion will not 


be further considered. The longitudinal motion is the solution of the following matrix equation: 
; —S{Z,f2dx D — Tf, 

j (D — Z,) SZ, hs —l s Ze 

—S\iZ, fo dx : 


oar D? — ft0, ferdx D - - | ' 
S Oufe dx (? SO, fe ie rN on? = Si Outs dx %. i= Si O¢fe dx 1¢ (26) 


— Sf M, fo dx D 
— oy erty _ 
My ere (I M,) | a | | M, | 


If the occurrence of the acceleration D7) in the central element of Eq. (26) is considered an inconvenience, it may 
= Dr. The augmented 




















be eliminated by augmenting the equations of motion by the additional equation 7 


matrix equation becomes 


p= £) ore 0 =i w Zz 


—S\iZ, fo dx 


— — NiO; f22 dx D — — 
ina Ce 2 Ix . ie 9 ) wa QV; 2 aX 0 QO fe x . 
Si Qufa dx & oger ' Si Qufe ds . SO; fads 1¢ (27) 


— Si M, fz dx D 
— ij Ww « iia ) ) ae - ; : 
: ere dx (L M,) T V, 
O —D 1 0 . 0 


Then, ¢ becomes a function of one or more of the de- 
6 

pendent variables. This can either be accomplished 
by a mass unbalance** incorporated in the control sys- 
tem or by an accelerometer that governs control action 


Eq. (27) has been applied in a practical case in order 
to attempt the duplication of observed flight experi- 
ence. The aerodynamic integral derivatives were com- 
puted by the theory presented here. The analytical 
results showed remarkably good agreement with ob- 
served performance. In particular, an oscillation of 
considerable amplitude was explained qualitatively, as 
well as quantitatively. The actual task of integrating Yo Sons 
the equations of motion was entrusted to an electronic ) , - + 
analog computer, and the results were obtained as trac- 
ings of an oscillograph. These results are shown in - 
Figs. 1 and 2 for two different flight conditions. In -0.5 aie tit de 


both illustrations, the upper curves are the results of 
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flight measurements and the lower are the results of the 
analysis presented in this paper. The scales are the 
same in the two illustrations—i.e., the peak acceler- 
ation in Fig. 1 is nearly eight times as large as that in 
Fig. 2. It will be noticed that the agreement between 
theory and experiment is extremely satisfactory in both 





cases. 
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It appears that the system will be stable if the char- 
acteristic equation of the system satisfies the stability 
criteria of the Routh discriminant. However, for cer- 
tain types of control motion the system may become 
unstable even if the Routh discriminant of the left- rob ‘fe oe ee 
hand side of Eq. (27) indicates stability. Thus, con- ; 
sider a control activation that is governed by the ac- _Fic. 2. Upper curve is the result of flight measurements for a 

. et 2 flight condition different from that of Fig. 1. Lower curve is the 
celeration as measured somewhere on the aircraft. result predicted by theory. 
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If the acceleration-controlled aero- 
dynamic surface is on one side of a node of the vibra- 


by electric means. 


tion mode and the acceleration-sensitive unit is on 
the other, then the accelerations of these two com- 
ponents will be 180° out of phase with respect to 
one another. Now assume that the controlling de- 
vice is adjusted so as to counteract any accelera- 
tion experienced by it. This will generally be the case 
if acceleration control is introduced to produce sta- 
bility. In that case and for the unfavorable location 
of the two units in question, the control action now ex- 
cites the vibration mode at resonance and the system 
becomes unstable. This same effect can be produced 
with favorably located components if there is a phase 
lag equal to an integer multiple of one-half the vibra- 
tional frequency between control command and con- 
trol response. This effect may be analyzed by the 
equations derived here; it is merely necessary to insert 
the proper control forces on the right-hand side of the 
equations of motion. 


CONCLUDING REMARKS 


A theory of the flight dynamics of cylindrical bodies, 
nearly circular in cross section and elastic in a limited 
sense, has been presented. A large number of restric- 
tions were introduced in order to render the treatment 
manageable. These might be grouped into those 
whose removal would render the treatment more tedious 
but not theoretically more difficult and into others that 
make any analytical treatment quite forbidding. Among 
the restrictions that could be removed without causing 
undue difficulty are those relating to the elasticity and 


to the air forces. 
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The admission of elastic deflections decomposable 
into m + 1 normal modes rather than only one would 
merely add m equations of motion to the system and y 
Thus, the 
consideration of three normal modes instead of one 


terms to the equations used in this analysis. 


should certainly be considered feasible. 


Consideration of the air forces arising from unsteady 
flow is possible since these air forces are known. Ip 
fact, it appears that even transient solutions could be 
obtained rather easily if use is made of analog or digital 
computing machines. However, when nonstationary 
air forces are introduced, the problem appears to be- 
come numerically too complex for solution by means 
of desk computers. 


All restrictions except the two just mentioned were 
introduced in order to retain differential equations that 
Theo- 


retical predictions of the flight dynamics of the air- 


are linear and have constant coefficients. 


craft under consideration do not appear hopeful if these 
restrictions are relaxed. However, an analysis in the 
neighborhood of the problem solved here can be carried 
out by perturbation methods. 


REFERENCES 


' Pringle, J. W. S., The Gyroscopic Mechanism of the Halteres 
of Dipteria, Phil. Trans. Roy. Soc. London, 233B, pp. 347-384, 
1948 


? Phillips, Wm. H., Effect of Spring and Gravity Moments on 
the Longitudinal Stability of the Brewster XSBA-1 Airplane, 
N.A.C.A. Wartime Report L-263. 


3 Rosenberg, R. M., Spring and Gravity Moments in the Con- 
trol System of an Airplane; presented at the Annual Summer 


Meeting of the Institute of the Aeronautical Sciences, July, 1944 





Simplified Treatment of the Turbulent Boundary Layer Along a Cylinder 
in Compressible Flow 


(Continued from page 28) 


? Tetervin, N., Approximate Formulas for the Computation of 
Turbulent Boundary-Layer Momentum Thicknesses in Compressible 
Flows, N.A.C.A. A.C.R. No. L6 A 22, 1946. 

3’ Kurzweg, H., A Few Aspects of Future Supersonic Windtunnel 
Design and Test Techniques, Symposium on Experimental Com- 
pressible Flow, U.S. N.O.L., White Oak, Md. NOLR 1133, pp. 
103-119, 1949. Also, Interrelationship Between Boundary Layer 
and Base Pressure, Journal of the Aeronautical Sciences, Vol. 18, 
No. 11, pp. 748-748, November, 1951. 

* Chapman, D. R., An Analysis of Base Pressure at Supersonic 
Velocities and Comparison with Experiment, N.A.C.A. T.N. No. 
2137, 1950. 

5 Jakob, M., and Dow, W. N., Heat Transfer from Cylindrical 
Surfaces to Airin Parallel Flow with and Without Unheated Starting 
Section, Transactions of the A.S.M.E., Vol. 68, No. 2, pp. 123 
134, 1946. 

6 Wilson, R. E., C., and Thompson, M. J., 2nd 
Interim Report on Experimental Determination of Turbulent 


Young, E. 


Boundary Layer Characteristics at Supersonic Speeds, Defense 


Research Lab., Univ. of Tex. Report CM-501, DRL-196, 1949. 


7 Rubesin, M. W., Maydew, R. C., and Varga, S. A., An Ana- 
lytical and Experimental Investigation of the Skin Friction of the 
Turbulent Boundary Layer on a Flat Plate at Supersonic Speeds, 
N.A.C.A. T.N. No. 2305, 1951 

* Wilson, R. E., Turbulent Boundary-Layer Characteristics at 
Supersonic Speeds—Theory and Experiment, Journal of the Aero 
nautical Sciences, Vol. 17, No. 9, pp. 585-594, September, 1950 

® Schultz-Grunow, F., Neues Reibungswiderstandsgesetz fuer 
glatte Platten, Luftfahrtforschung, Vol. 17, pp. 239-246, 1940 
(Translation, N.A.C.A. T.M. No. 986.) 

” Spivack, H. M., Experiments in the Turbulent Boundary 
Layer of a Supersonic Flow, North American Aviation Report 
CM-615, 1950. 

11 Seban, R. A., and Bond, R., Skin Friction and Heat Trans 


fer Characteristics of a Laminar Boundary Layer on a Cylinder in 


Axial Flow, Univ. of Calif. Contract No. W33(038)-ac-15229 
with the A.A.F., A.M.C., 1949. Also, Journal of the Aeronaut- 
ical Sciences, Vol. 18, No 10, pp. 671-675, October, 1951. 

12 yon Karman, Th., Ueber laminare und turbuiente Reibung, 
Zeitschrift fuer angewandte Mathematik und Mechanik 1, pp 
233-252, 1921. (Translation, N.A.C.A. T.M. No. 1092.) 





TI 
orde 
the | 
for t 
sure 
the | 
cont 
the | 
free- 

T 
Eve 
exte 
mat 
out 
that 

A 
tent 
the 
in t! 
The 
shaq 
Tab 
mer 


T 


thr 
floy 
VOI 


Fe 


sur 
str 


the 
de; 
de: 
sid 
TI 
tin 
mit 
19¢ 
shi 


ing 
To 





Sable 
would 
and 
S, the 
f one 


‘eady 

In 
ld be 
igital 
mary 
» be- 
jeans 


were 
that 
hec )- 

air- 
hese 

the 
ried 


5 on 


ane, 


“on- 
mer 
44 


Tr 


na- 
the 
ds, 


at 
ro 


ter 


+() 
ry 


rt 


in 


29 








Velocity and Temperature Distribution 
Through the Laminar Boundary Layer in 
Supersonic Flow’ 


CARLO FERRARI+ 
Cornell Aeronautical Laboratory, Inc. 


SUMMARY 


The laminar boundary layer in a compressible gas is studied in 
order to obtain the velocity and temperature distributions through 
the layer in the presence of any pressure distribution whatsoever 
for the external flow, under the usual assumptions that the pres- 
sure is constant through the layer and that the temperature of 
the constraining surface is constant throughout its extent. The 
contour of this surface is supposed to be so gently turning that 
the surface elements may always be considered aligned with the 
free-stream flow. 

The method described is well suited for routine computations. 
Even though it makes use of successive approximations, only one 
extensive recomputation is required; after the second approxi- 
mation is obtained, further approximations, if desired, are carried 
out by means of short cuts that reduce the labor involved to less 
than one-tenth of that required in the first recomputation. 

A significant effect of scale on the pressure changes in the po 
tential flow exterior to the boundary layer is observed, although 
the magnitude of tne changes in these external pressures is small 
in the example selected for numerical illustration of the theory. 
The numerical example shows the progressive change in profile 
shape up to the point where imminent separation can be expected 
Tables and aids to calculation are provided, along with the nu- 


merical example 
INTRODUCTION 


[ew REPORT DEALS WITH A NEW METHOD for the 
determination of the velocity and temperature 
through the laminar boundary layer in supersonic fluid 
flow. This method is related to one expounded by 
von Karman and Tsien! and to the ideas suggested by 
Feinsilber,? but the present treatment may be applied 
to a compressible gas, and it holds for any kind of pres- 
sure variation along the boundary in the external 
stream. 

The method discussed here is actually one involving 
the technique of approximations. The 
degree of approximation associated with the solution 


successive 


desired in this investigation, however, may be con- 
sidered as more than adequate for practical purposes. 
The procedure described is, besides, well suited for rou- 
tine evaluations because of the fact that the numerical 
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values of the functions entering the final formulas have 
already been tabulated in standard works. 

The assumptions upon which this method is based 
are the same as those that have been employed pre- 
For instance, it is 


‘ 


viously in similar investigations. 
taken to be permissible to use the “‘simplified’’ equa- 
tions of the boundary layer in order to describe the vis- 
cous flow therein, and, hence, it is implied that the 
static pressure is constant throughout every cross sec- 
tion. The experimental results of Ackeret, Feldmann, 
and Rott* controvert this assumption, of course, more 
and more as the particular cross section under con- 
sideration is selected closer and closer to the one at 
which the shock wave starts. 

Such an assumption of constant pressure, on the 
other hand, does not preclude the derivation of a for- 
mula from which it is possible to deduce how the pres- 
sure within the stream external to the boundary layer 
is modified by the presence of the boundary layer itself. 

In a preceding report of this series,‘ it has been 
shown how, in the case where the Prandtl Number is 
taken to be unity, the problem can be solved approxi- 
mately through use of a step-by-step procedure by 
means of which the pressure variation through the 
cross section is taken into account at each step of the 
calculation processes. 

Further, it is assumed here, even as was done in the 
above-mentioned report, that the thermal conductivity 
of the wetted wall is infinite (so that the temperature 
of the wall is constant over its whole length), and, in 
addition, the curvature of the wall is taken to be so 
small that its surface elements may be considered co- 
incident in direction with the x-axis and the outward 
normal to be in the direction of the y-axis. 

A numerical example is presented to illustrate the 
application of the suggested method, and the necessary 
tables employed in the computations are also appended. 


PRELIMINARY CONCEPTS 


(1) Symbols 


U = component of the velocity along the x-axis 

V = component of the velocity along the y-axis 

Vi = limit velocity attained by the free stream when ex- 
panded toa vacuum 

Ue = velocity of the undisturbed stream 
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velocity of the stream on the outer edge of the boundary 
(for all intents and purposes, this vector is taken to 
be parallel to the x-axis) 


U/V; v=V/Vy; u.=U./Vy uo =Ue/V; 


= density of the fluid; p, = density of the fluid in the ex- 
ternal stream; po = density of the fluid in the un- 
disturbed stream; pp = density of the fluid in con- 
tact with the wall 

= pressure 


absolute temperature 
viscosity coefficient 
kinematic viscosity coefficient 





The above symbols with subscripts e, ©, and p have analogous 
J g 


meanings to those given for p. 


L = length of the wall; & = x/L; n = y/L 
Rw = Reynolds Number of the asymptotic stream (undis- 
turbed flow at infinity) = Uo.L/vo; n* = nV Ra 


jf + you _ dp > ( 
pe” Ox F Oy r dx oy 


Energy equation: 


The eontinuity equation is satisfied by writing 


pu = pu. (Oy~i/On) | 


—patt, (Ox/de) f (3) 


B 
II 


Through the use of the transformation previously em- 
ployed by von Mises,’ and by letting 


Z = M/be (4) 
where 
gs = (7/T.)’; p/pp = T./T = 27%" (5) 
and 
*t we ’ 
ePe j , 
f= I Ik; wt =VRoh 6) 
0 bapa 
while also 
Z = (u,2/2u..2) — (u?/2u,2)) (6’ 
: , )’) 
=e] 24" f 


then from Eq. (1), one obtains 


OZ/de, = (u/u.) (0/OW*) [27 ?/" (0Z/Oyv*)] + 
(1 — 2'/")(dZ,/dé:)] (7) 
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Y = stream function (dimensions of a viscosity coeftlicient 

Wi = V/(poloL); fi* = Ww Ro 

2 = p/mes 2 =1—s 

Z = (te?/2u (u?/2uam*); Ze = te?/2um?; and 
$= Z -—-Z = u?/Quon? 

Cy» = specific heat at constant pressure 

y = adiabatic exponent (= 1.4 for air) 

P, = Prandtl Number (= 0.75 for air) 

n = a coefficient defined by u/ue = (7T/T.)" (n = 0.75 for 
air) 


Other symbols are specified in the text. 


(2) General Equations 


On the basis of the fundamental assumption con- 
cerning the constancy of the pressure as mentioned in 
the introduction, the equations defining the flow within 
the boundary layer are as follows: 

Equation of motion: 


ad _ oy, Ws 4 ad = P 
P oy/ aa” oy , oy " 


( i 4 v—) > yo 4 & > ( ~~) 4 (~) 

“— ae” oy dx P, oy “ Ov . Oy 
7 VU dU’, 1 & ° ( or 4 (S-) : 
on ae” oy ‘a ov ‘ oy ” 


Applying the von Mises transformation to Eq. (2), one 
obtains in an analogous manner: 


02/08 = (1/P,) (1 — m)"~?"" x 

(0/Oy~r*) [(u/u.,) (O2,/O~r*)] — °°" ~?)/" (nu. /u) X 
[2u..?/(1 — u,”)](OZ/Oy*)? (8) 

where 2; = 1 — sg. 

For Eq. (7): 


The boundary conditions are: 


Z = Z, tory," = 0; Z= Ofory,* = @~ (7’) 
For Eq. (8): 

2, = 1 — (7,/T,)" for y,* = Ol 

2: = Ofor yi* = of 

wherein 7), the constant temperature at the wall, is 

supposed to be known for this phase of the calculations. 


SOLUTIONS OF THE TRANSFORMED EQUATIONS 


(3) External Solution of the Equation Derived from the 
Equation of Motion Through Use of von Mises Trans- 
formation 
Within the boundary-layer region, but near the ex- 

ternal stream, it is permissible to write 
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VELOCITY AND 
0Z/d~i* = 3 = O and u/u,, = u,/u. (9) 
and letting 
TT I [ Me ae (10) 
| = &1 
. PJo uu, 
it follows from Eq. (8) that 
02,/0&,** = 072, /Oy*? (11) 


which is identical to the equation of heat conduction 


within a medium having unit diffusivity. Let 


a) 
m 


¢ eT eee **KM P 
(Siy,%=9 = Ay = > Afi (12) 


then the integral of Eq. (11) which satisfies condition 


(8’) is 


z= > Am2?m! §**" i” erfe [Yi*/(2V &**)] (13) 
m 


(4) External Solution of the Equation Derived from the 
Energy Relation Through Means of the von Mises 
Transformation 
In order to determine the external solution of Eq. (7), 

the result just obtained as Eq. (13) will be used di- 


rectly. Within the region of the boundary layer now 
n. pine ts 90” | 00 oot, 9 ; 
F = - )A,,2*"mié, i erfc = = YosB.E 
nm “Vv & s 
wherein 


C,. = (A,,/n) 


B,' = sB, 


Thus, it appears feasible to write 


9 
2v & 


m 


* 
alles Tat] CB ar + C,B,,' 7? erfe vs _ 4 + C2B',—1 14 erfe 


Now return to consideration of Eq. (15). Applying to 
it the Laplace transform with respect to the variable 
£,*, one obtains 


(d?Z/dy,**) — pZ = —F (20) 


= L(Z) = J e pei® Zz dt,*\ 
0 


= L(F) 


in which 


NI 


>» | 


The Laplace transform of such an expression as 
f*7" erfe [Wit P, /(20/&*)] is discussed in Appendix 
V of The Conduction of Heat in Solids by Carslaw and 
Jaeger. Using this result, it is noted that 
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(n I)/m (. 1/3 


7 


under consideration, one may take z = 2 
for air) to be constant, because it varies but slightly 
In fact, there can be no objection to 
In addi- 


within the layer. 
selection of this constant as precisely unity. 
tion, letting u/u., = u,/u., Eq. (7) may be written in 


the form 
OZ/dE, = (u,/u.) (0?Z/OW*2) + (1 — 2'/") (dZ,/dé) 


or also, if we set 


£,* - P,§,°* (14) 
then 
r or 9 1 y 
OZ o£,* = (0°?Z Oy,*?) + (l—3z 7 (dZ, dé,*) 
~ . l/n : a . 
Expanding | — z‘” in a power series in the neighbor- 
hood of z = 1| and retaining the first degree term only, 


one obtains 


OZ 0°?Z 1 dZ, 0°Z .F 1. 
Z = : (Lo) 
n ' dé* oy *? . 


o%* di*? 
Let 
Z. = > B,t** (16) 


then Eq. (13) yields 


1) _ . * 
= 2, CA” 


m 


* > 
Ei a. E | +B E*°-” (17) 


9 * 
2vV bi . 


22m! (1/P,™) 
l (18) 


/ 


Viv P, 


9 * 
2V fi 


+...+ 


‘i 
‘ ss ; P; 
C,, By?” erfe v1 v = | (19) 


») 
“V vu 


Lh e*i® erfe [Wt P,/(2V78*)]} = 

i(—ti” r 4'] (d™ "/dp™ " (e w*a/p' +r) (21) 
for m > r and wherein g = VV P,p. It follows that 
—F turns out to be a summation of terms ot the type 


. qvi*® 1+1/2 f Q«< 
Gc, @”’ /s yi? (22) 


where g and ¢ are positive integers. It follows that one 


may write 


Em 2 des (23) 
t 


g 


wherein Z, , is defined by 
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(d*Z,9/dyr*?) — pLig = Cig*hi™ (e~™"/p't*) (24) Z,,=0 (25) 
for ¥,* = 0, as well as for ¥,* = @ must be kept in 
where the conditions that mind. 


One deduces, consequently, that 


2 7 c. F 1 E _ > n(2) "8 h _ e V pPfi* g! 
ie QQ re™ aft — (1/4/P,)) i=o0 NA/ "1 —(1/V PP) gl — 1/VP,) 


C, .* K . *e—h e7% pix “ 
ji t, 9 — Ga? h! () V1 — , o! (96) 
2(a/ P, — 1)pr’* h=0 kh} (VP, — 1)"p"" —p’*(n/ PP, — 1)*° 


Then noting that 


L* 
{ oni (2+t+h)/2 -24t+h *./P, 
Lo" Wiest 5)) = (46*)* itil eet wiv 
\ pit arti / ie 
and 
L- (e~ Vine pg tietite a (4€,*)Ct+)/2 2+! +8 rfc [yi* (20/§,*)] (27) 
one obtains 
g *o—h 
; (fy ; 
“ C," . h jew) 2HIt+h)/2 2+ttHe ge VN Ge g! *\(2+1+ 2)/2 
Liz => 9 ; ) > ri $l ) 1 eric 9 — — = (4, ) 4 
2(\/P, — 1) Lixo (+/P, — 1) 2V i ivr, — 3) 
h 
c yi" g 
* > 4 m(£\( Je x * 
2+ite  u¢ v1 i Ci, p 2 h “Vii e &\(24+1+h)/2 (2+t+h ~V~, 
1 erfc ; = 7 —| h=0 y i (4,*) 1 erfc — 
2Vii* 2(VP, — 1) (VP, — 1) 2V b* 
z! ea) te ei eee _* 
g , (4b *) t+e2 a t ® erfc y (98) 
(VP, — 1) 2vV i:* ss 


The expression for Z which satisfies Eq. (15), and which, for ¥,* = 0, coincides with the Z, given by Eq. (16) is 
therefore 


* 
‘ i sf V1 , 
Z = >Bn2?"m!,*"12" erfc +> Z,, (29) 
94/¢.¥* , 
m “Vv $l 1.8 


Thus, for example, in the most simple case in which 


Le == By + B,&,* (16’) 
one has 
+ * *x\l+h 
P a - (4&,*) 
Z = By erfe - v1 = + B,dt:*i? erfe ~ v1 =+ DY C,.*—— 
2V & 2V & t=0,2,4,... AVP, — 1) 
* oo 
& “erfc es = ak lies “ | (29) 
: 2vV .* 2V i* 
wherein ¢; = ¢/2, and 
cy Pa = 0 for ¢ = B :: 5, etc. 
Ci, 0* = —(1/4") (A;,/m) (2 th! B,/P,") for t = 0, 2, 4, ete. (30) 


(5) Internal Solution of the Equation Derived from the Energy Relation Through Means of the von Mises 
Transformation 


In proximity to the solid wall (where ¥,* = 0), Eq. (7) may be written in the form 


vA , (n—1)/n]- 5 (n—1)/» oF p 
= = {(u/u.) [(n — 1)/n]2!"—?/")—"az, Ovi*) | y*-9(OZ/O¥i*) + (u Uw)Sy"~/"(0?Z/Oy,*2) + 
1 


(1 — ,'/")(dZ,/dé) (31) 
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But it is true that 


}(u/u.) [(n — 1) ‘n}e'"—Y/"I-1 (ag Ovi") f yi .o = [(n — 1)/n](p./pe) (02/0n*), «0 (32) 


Now this term vanishes for m = 1, whatever the heat exchange AQ through the wall happens to be, and it vanishes 
as well for the case in which P, = 1, regardless of the magnitude of n, if there is no heat exchange (AQ = 0). 
Therefore, it seems allowable, in the case in which AQ = 0, for fluids that, like air, have values of n and P, both 
close to unity, to disregard the term within brackets. In the general case, we may still consider this term to be 
so small as to permit, without appreciable error, the substitution of the constant value taken on at the wall by 





the derivative 0Z/Oy,* for the actually variable value of this quantity. 


Eq. (31) the expression 


= 


SS 


wherein 


Dy = [1 — n)/n] (p. 
—(0Z/0~1*)y.*-0 


II 


In the development given below, the 


Consequently, one may substitute into 


| yi * = DoVYo 
lyr ° oy," i L lyx=o 


pe) (02/On*),*-ol a 
j (33) 


( 


solution of Eq. (31) is obtained by following a procedure analogous to that 


applied by von Karman‘ to the study of the same problem in the case of incompressible flow. Let 


f=Z2,-Z2 (34) 


then Eq. (31) yields 


0% /dyi*? = —(1/¥ 2) (1/V 9) (z,'/" 


1 Z,\" 1) "\ (dZ, dé,*) \ ‘9 V/Z, 4 
[1/(+/22,"-?/)] (A/V) [(08/08:*)V/2Z + DoYo] (31’) 


Analogously to what was done by von Karman in the above-mentioned report,® the latter equation is replaced by 


O%¢/dpi*? = (A*/V/e) [1 — VE/tr 


wherein 


A* = Ao*(1 + Y,A;:*): Ao* = 


gyi/2qr l/n n 1)/n 
—(Z,'"LZe'S )/ Sp» ; 


(31°) 


A,* = Do(1/V2) [1/(2,'/"Z.’) ]; 
Z,’ = dZ./dt,* (35) 


while ¢, denotes the value of ¢ which corresponds to the inflection point in the graph of ¢ vs. ¥i*/(2+/ &*) obtained 


from the internal solution. 


¥ 


Jaresin (24*/V 4A *? + 2Y,2A* 


(6) Join Between the External and Internal Solution of 

Eq. (7) 

The connection between the solution given by Eq. 
(36) and the one represented by Eq. (29) is made by 
following the same procedure as that suggested by von 
Karman. From the graph of ¢ vs. wWi*/(2V&*), 
plotted through means of Eq. (29) for the various 


parametric values of £,*, the value of 
H = —}(02Z/d[wi*/(2V8*)) f U/(Wa*) Ds =e, (37) 


is determined. 

Hence, by imposing the condition that, for ¢ = ¢,, the 
value of Of/O0Y¥,* given by Eq. (36) must be exactly 
equal to 77, one obtains the condition that 


H? = Y,? + 2A*y/5, = 
Vo? + 20/f, Ao* (1 + YoAi*) (37’) 


Vt) — aresin [((2A* — 


From Eq. (31”), one deduces that 


h* = (Wf A*)[Yo — V Yo? + 4A*y dige 2A*(E/V ot) ] + (WV 1¢,2/VA*) X 


2A*V ¢/¢,)/V4A*? + 2V2A*/Ve]} (36) 


This relationship permits one to evaluate Yo. Then, 
for each £,*, Eq. (36) gives the value of ¥,* in terms of 
¢, provided that D)—that is, {(0z/0n*) |,«~9—is known. 


The complete diagram of ¢ in terms of ¥,*/(2+/¢,*), at 
£ b Vv $l 


every given cross section of the boundary, £,* = const., 
is now determined, because for the interval 0 Z ¢Z ¢,, 
it is given by Eq. (36), and for the interval ¢, Z ¢Z Z,, 


it is defined by Eq. (29). 
DETAILED PROCEDURES OF SOLUTION 


(7) Second Approximation for the Temperature Distribu- 
tion—i.e., the Law of Variation of the z Values 


In order to determine a second approximation for the 
distribution of z through the boundary layer, one may 


proceed as follows. - Rewrite Eq. (8) in the for ..1 
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(u/u.)o~?"—~Y/"(O%,/08:) = (1/P,) (u/s) (L — 21) ~~" (0/OP*) [(u/uw) (O%/O*)] — 


Then let 


F(z) = f? "dz 
0 


and so, Eq. (8”) becomes 


0?F/On*? = — 2nP, [2u..2/(1 — u.”)|z [(0/0n*) VZ, 


In addition, Eq. (3) may be recast as 


yi* T 9./e* 

pal. 2Vt 

¥ si f ; dy,* = a 
. pl V/2 


The forms of Eqs. (40) and (39) now suggest a pro- 
cedure of successive approximations for the determin- 
ation of the manner in which z varies within the layer. 
If one introduces into the second member of Eq. (40) 
the value of z obtained from the external solution, 
considered as a first approximation, then Eq. (40) per- 
mits one to obtain n* in terms of ¥,*/(2+/¢,*), for each 


value of £;*. Now evaluate the quantity 


Bin — (2) OZ 
aye ee - 2 V2 \b0/ O[vi*/(2V :*)] 2VE* 


by making use of the values of z given by the external 
One is then able to evaluate ¢ in terms of 
n*, and, hence, obtains 


n* n* 
F= f ant [ gdn* - $1, 0 n* + Fo (41) 
0 0 


Pf. o = (OF ‘On*),,* 9 = z,| "(Oz On*),« 0 


solution. 


wherein 


and 
Fy = (F),*=0 = [n/(n + 1)]2,"7?” 


Because Eq. (41) permits the evaluation of F, then one 
has simply that z= [(7 +1) /n]"/"*? F"’“"* Pasa second 
approximation for z. It is easy to see that this proced- 
ure may be extended in order to obtain further approxi- 


mations. 


(8) Practical Method for the Determination of the Mo- 
tion Within the Layer. Evaluation of the Heat Trans- 
mission 
The calculations are generally laborious if the amount 

of heat AQ transferred from the fluid to the wall is taken 

to be a given quantity; they are simplified if one fixes 

In fact, under 


7;, and consequently evaluates AQ. 
these conditions, a first approximation to the distribu- 
tion of the velocity and temperature is determined 
within the layer by taking Dy =-0 (this value comes 


n[2u..°?/(1 — u,?)] [0Z/Oy,*]2  (8") 


n 
(n+1)/n 9 
= 5 38 
+ i (38) 


~ ZP — Pilp/s.)* 2V 2. VZ, — 201 — &) “xX 
(02,/0&:*) = @ (39) 


p. [vr/2evar gl/" Wn" 
a mm 1/2 d 9 Poa (40) 
Pe 0 (Z, — Z) “Vv i 


nearer to being the true one as m approaches unity). 
The procedure just outlined permits one to determine 
the values of (02z/0n*),*-9 if the above-mentioned as- 
sumptions are employed. Once having obtained this 
temperature gradient, a second approximation for the 
values of Z and then of z may be deduced if necessary. 
It should be noted that the values of Do affect the in- 
ternal solution of Z only. It is not necessary, there- 
fore, to repeat the calculations leading to the external 
solution of Z, which, in fact, are the only laborious ones, 
in order to obtain a second approximation for 2 which 
will be satisfactory throughout the layer. 
One obtains finally that 


*1 _ 
AQ = y) ie (5) dé = 
0 On 0 


” 


Ps ‘7 =) 
a Ke - 1é (42) 
Cpbpt p n Vv f Zp (= uit ( 


corresponding to the selected value of 7;. 

If the system is thermally insulated (AQ = 0), a 
first approximation determination of 7), may be made 
by taking the value of the latter to be that correspond- 
ing to the condition for flow along a flat plate 
[(Op/d&) = O]. In this case, the following equation 
describing the law of variation of the temperature, holds: 


(T — T,)/T, = —(P,/2) U2/(e>T>) 
[P,/(S0y2)] (1 — P,)[U2/(¢pT>) (U/Us)2 +... (43) 


Therefore, it appears that, in the case now being ex- 
amined, a parabolic law of variation of 7/7, in terms of 
U/Vc,T, is valid throughout almost the whole extent 


of the layer. If one retains only the term 


(T — T,)/T, = —(P,/2) U2/(GpT>) — (48) 


and imposes the condition that Eq. (43’) should blend, 


for U’ = U., into the relationship that holds for the 
external variation of 7, which is, of course, 7./7» = 
1 — [U.2/(2c,7»)], wherein 7) is the stagnation tem- 
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TABLE | 
Complete Velocity Distribution Through Boundary Layer for Three Selected Locations 


,* 
= ¢ 2V¢ n* ¢ Vv 2¢ n* ¢ 2Vt n* 
av i* £ = 0.1 £ = 0.1 & = 0.] ¢=03 t= 0.3 ¢= 03 —& = 0.5 & = 0.5 &= 0.5 
0.00 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 0.000 
0.05 0.027 0.230 0.420 0.024 0.211 0.712 0.021 0.190 0.921 
0.10 0.053 0.327 0.596 0.048 0.304 1.010 0.042 0.277 1.304 
0.15 0.080 0.400 0.726 0.072 0.376 1.227 0.064 0.347 1.580 
0.20 0.106 0.461 0.833 0.096 0.439 1.404 0.086 0.405 1.805 
0.25 0.182 0.514 0.926 0.120 0.491 1.537 0.108 0.466 1.996 
0.30 0.158 0.562 1.008 0.144 0.537 1.693 0.130 0.511 2.163 
0.35 0.182 0.604 1.082 0.167 0.579 1.816 0.152 0.552 2.314 
0.40 0.206 0.642 1.151 0.191 0.617 1.928 0.173 0.589 2.454 
0.45 0.229 0.677 1.214 0.212 0.651 2.032 0.195 0.624 2.583 
0.50 0.251 0.709 1.273 0.233 0.683 2.130 0.215 0.655 2.703 
0.55 0.272 0.735 1.329 0.253 0.712 2.221 0.234 0.685 2.816 
0.60 0.292 0.765 1.381 0.273 0.739 2.308 0.253 0.711 2.923 
0.65 0.311 0.789 1.43 0.291 0.763 2.390 0.271 0.736 3.024 
0.70 0.329 0.811 1.479 0.308 0.785 2.469 0.287 0.758 3.12) 
0.75 0.345 0.83 1.525 0.324 0.805 2.545 0.303 0.775 3.195 
0.80 0.361 0.849 1.569 0.340 0.824 2.617 0.318 0.797 3.284 
0.85 0.375 0.866 1.612 0.353 0.841 2.687 0.332 0.814 3.370 
0.90 0.388 0.881 1.653 0.366 0.856 2.755 0.344 0.830 3.453 
0.95 0.400 0.894 1.693 0.378 0.869 2.821 0.356 0.844 3.534 
1.00 0.411 0.906 1.732 0.389 0.882 2.885 0.367 0.856 3.6138 
1.10 0.429 0.926 1.807 0.407 0.903 3.009 0.385 0.877 3.765 
1.20 0.444 0.943 1.880 0.422 0.919 3.128 0.400 0.894 3.911 
1.30 0.456 0.955 1.950 0.484 0.932 3.244 0.412 0.908 4.053 
1.40 0.465 0.965 2.019 0.444 0.942 3.356 0.421 0.918 $.191 
1.50 0.472 0.972 2.085 0.450 0.949 3.466 0.428 0.925 4.326 
1.60 0.477 0.977 2.151 0.456 0.955 3.574 0.4384 0.981 4.458 
1.70 0.481 0.981 2.216 0.459 0.959 3.681 0.437 0.953 4.589 
1.80 0.484 0.984 2.281 0.462 0.961 3.787 0.440 0.938 4.719 
1.90 0.486 0.986 2.345 0.464 0.963 3.892 0.442 0.940 4.848 
2.00 0.487 0.987 2.408 0.465 0.965 3.996 0.448 0.942 4.976 

perature of the potential flow external to the layer, wherein 
then it is found that 
x= W*/(2Vi*) 
T,/To = 1 — (1 — P,) 4? (44) 
(10) Influence of the Boundary Layer on the External 
Gross BOUNDARY-LAYER EFFECTS Pressure Distribution 
9) Thickness of the Boundary Layer From the continuity equation, one has 


The thickness of the boundary layer may be obtained 
eS : (0/OE) (pu) + [(0/O~1) (pv) | [pu/(p.u.)] — 


from Eq. (40) by inserting in it, for each £, the value 

, : ; Jah: ‘is (0/01) (pu) | [pv/(p.u.)| = O 
of W*/(2Vé,*), for which ¢”? = 0.997Z,"*. This | v1) (ou) | [o0/(o | 
value varies but slightly with ¢, and in the numerical — or 


example that follows, it is seen to be close to 2. 
(0/0E) [(p/p.) (u/u.)| + 


Thus, it may be presumed that 
(p?/pa*) (u?/u..”) (0/O~) (v/u) = O 


"3 1/) 


6* = V D¢,# Ps / = =e (45) 
0 ¢ 7 


Des whence, it results that, within the layer, 


v v “Wp u? O “‘@ pth a ' 
() (0), = LE a(n LR Ean 
u/; u/s o 0 p’u OF \p..u 0 OE pu 


If the pressure variation due to the change in direction of the velocity is assumed to be that predicted on the basis 


of linear theory, (pb — p.)/(p2U.*) = (1/V M..2 — 1)Band if 6B denotes the angle between LU’, and U.,, then the 


change in pressure is given by 


Ap, | v v\_ l "W5D Ug 
-meel)-Ol-wia Rt» 
7 VM? -1 u/s u/ VM.? —1J0 Of \ pu 


poU 
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wherein Ap, = p. — 


the latter is taken into account. But itis true that 


-JANUARY, 


1952 


be, and p, denotes the pressure on the outer edge of the boundary layer when the presence of 


“W5O Pott. 2V/t* MePelte ” (x) ge O Datla 
dy, = dx 
0 Of \ pu "Keo Baodattea J0 0f:* \ pu 


wherein (x)5« is the value of x that corresponds to n* = 6 


Ap. l 


It follows that 


2V/t1* Uepelle | i ra) (t= ) 
Se a i: , : dx 
PoU wo" V M..* — 1 WR babatl 0 O£,* pu 


It appears that the influence of the boundary layer 
on the pressure distribution in the external stream is 
dependent upon the Reynolds Number; 
while in natural scale the influence may be small, under 


therefore, 
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Fic. 1. Successive steps in the numerical procedure of estab- 


lishing the correspondence between the velocity-related param- 
eter 2;,, and the adjusted location along the constraining wall, 
t,** 

$l 





experimental conditions, which usually correspond to 
smaller Reynolds Numbers, such influence may be ap- 
preciable. 


CALCULATED DATA 
(11) Numerical Application 


The procedures described in the preceding para- 
graphs are now to be applied to the following example. 


whence u,.7 = 4/9: Z, = 0.5(1 — 0.2é,*); » = 0.75: 
and FP, = 0.75. 


Only the results of the calculations are presented be- 
low. The subsidiary steps are to be found ‘‘in extenso” 
in a Cornell Aeronautical Laboratory report’ that is 
available upon request from the Johns Hopkins Uni- 
versity, Applied Physics Laboratory. 
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Complete velocity distribution through the boundary 
layer at three selected locations. 
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Fic. 3. Comparison of first and second approximations to the 
temperature profiles encountered at three selected locations 
spaced along the wetted wall. 
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: TABLE 2 
First and Second Approximations to Temperature Factor 


yi" ps p - . in 
“I “II “ “II 7 “11 
20/t* &= 0.1 &é= 0.1 &é=Q3 &=03 & = 0.5 &= 0.5 
0.00 1.405 1.405 1.370 1.370 1.337 1.337 
0.05 1.386 1.463 1.354 1.353 1.324 1.305 
0.10 1.367 1.483 1.338 1.345 1.311 1.291 
0.15 1.348 1.493 1.322 1.337 1.297 1.281 
0.20 1.329 1.497 1.306 1.330 1.283 1.271 
0.25 1.310 1.496 1.289 1.322 1.269 1.262 
0.30 1.292 1.492 1.276 1.313 1.255 1.253 
0.35 1.274 1.485 1.257 1.304 1.241 1.244 
0.40 1.256 1.475 1.241 1.295 1.226 1.235 
0.45 1.239 1.464 1.225 1.285 1.212 1.226 
0.50 1.222 1.453 1.210 1.275 1.199 1.217 
0.55 1.206 1.438 1.196 1.264 1.185 1.208 
0.60 1.191 1.422 1.181 1.253 1.172 1.199 
0.65 1.176 1.404 1.168 1.242 1.159 1.189 
0.70 1.162 1.386 1.154 1.231 1.147 1.180 
0.75 1.148 1.368 1.142 1.219 1.135 1.172 
0.80 1.135 1.349 1.130 1.208 1.124 1.163 
0.85 1.123 1.330 1.118 1.197 1.114 1.154 
0.90 1.112 1.310 1.108 1.186 1.104 1.145 
0.95 1.101 1.291 1.098 1.174 1.094 1.136 
1.00 1.091 | .272 1.088 1.164 1.085 1.127 
1.10 1.074 1.235 1.071 1.142 1.069 1.111 
1.20 1.059 1.200 1.057 1.123 1.056 1.096 
1.30 1.046 1.167 1.045 1.104 1.044 1.081 
1.40 1.036 1.137 1.035 1.086 1.034 1.068 
1.50 1.027 1.109 1.027 1.070 1.026 1.055 
1.60 1.021 1.083 1.020 1.054 1.020 1.042 
1.70 1.016 1.060 1.015 1.040 1.015 1.031 
1.80 1.011 1.038 1.011 1.026 1.011 1.020 
1.90 1.008 1.018 1.008 1.012 1.008 1.010 
2.00 1.006 1.000 1.006 1.000 1.006 1.000 


Proceeding upon the basis of the above-assumed 
values, one deduces: 


1/2 


T/T. = 18 — 16Z,; U./Ue = (1 — 0.2&*) 
: - / (1 + 0.16&*)~* (1 — 0.2¢:*)~°° d&* 
/0 


?,/p. = (1.8 — 1.6 Z,)* 


Fig. | shows the diagrams giving p,/p., &*, &** in 
terms of &. 
From Eq. (44) it follows that 7,/7, = 8/9, whence 


T,/T, = [1.6/(1.8 — 1.6Z,)]; >» = 
1 — [1.6/(1.8 — 1.6Z,)]°-7 


From this relationship for 2, ,, one deduces that 2, » 
may be fairly well approximated, in terms of £**, by 
the following equation: 


— 0.4226 + 0.12808,** — 


21, >= 2. A ati” a 
0.01344&**? + 0.001479£,**? 


m 


The complete velocity distribution within the bound- 
ary layer, throughout the three cross sections selected 
as the ones corresponding to locations = 0.1, € = 0.3, 
and § = 0.5, is plotted to give the diagrams of Fig. 2, 
and the values are listed in Table 1. 

Fig. 3 and Table 2 give the values of z as a function, 
respectively, of n* (distance normal to the constraining 
walladjusted for Reynolds Number) and of ¥:*/(2+/ &*), 
which is a parameter representing the stream function 
(nondimensionalized by reference to a viscosity coeffi- 


cient) that has been ‘‘reduced”’ to account for the value 
of the local velocity external to the layer. Hence, 
this curve essentially describes the temperature dis- 
tribution within the boundary layer as would be cal- 
culated as a result of applying the first and then the 
second approximation procedures to this determination 
of the temperatures. 


It should be remarked that, particularly for cases 
where & = 0.3 and — = 0.5, the values of z obtained as 
a second approximation are fairly close to those ob- 
tained as the first approximation. Consequently, one 
may safely draw the conclusion that the law of dis- 
tribution of velocity within the layer as well as the law 
of variation of the temperature therein, as now ob- 
tained in second approximation, is sufficiently accu- 
rate. At this point, it should be noted that departure 
of the true value of z from those assumed as a first 
approximation has but little influence on the graph of ¢ 
versus y,*/(2+/£,*), so that a continued process of cor- 
rection would involve solely the redetermination of y*, 
which is not a laborious task. 
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Computation of the Laminar Boundary 
Layer with Suction’ 


FRIEDRICH O. RINGLEBt 
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SUMMARY 


In this paper, the incompressible laminar boundary layer is 
computed using, as approximate profile for the velocity u tan- 
gential to the boundary, an exponential function of a polynomial 
of the fourth order in y with coefficients that are functions of x. 
If the velocity outside the boundary layer and the suction (blow- 
ing) velocity are given functions of x, the first of those coefficients 
is determined by an ordinary differential equation of the first 
order which can be solved by numerical methods in the general 
case and exactly by simple analytical expressions in many special 
cases. The other coefficients are determined by elementary 
functions of the first coefficient and its first derivative. This 
method, which is convenient and combines high accuracy with 
uniformity of the procedure, is applied here to the flat plate with 
and without suction, the stagnation point, the asymptotic suction 
profile, and some other examples of interest. 


NOTATION 


x,y = coordinates within the boundary layer tangential and 
normal to the boundary 

u,v = velocity components in the x- and y-direction within 
the boundary layer 

p(x) = pressure 

U(x) = velocity outside the boundary layer 

v(x) = suction (blowing) velocity 

7o(x) = shear at the boundary 

p = mass density of the fluid 

v = kinematic viscosity of the fluid 

a(x) 

b(x) = coefficients of the velocity profile 

Az) u = U(1 — e% t+ 9% + cx? + dy) 

d(x) 


INTRODUCTION 


— PRACTICAL VALUE OF boundary-layer suction, 
especially of continuous (porous) suction in order 
to diminish the profile drag and to increase the maxi- 
mum lift, proved by many experiments, led during the 
last years to a considerable number of investigations 
about methods for the computation of the laminar 
boundary layer with and without suction conditions. 
The literature up to 1938 is discussed in detail in S. 
Goldstein’s Modern Developments in Fluid Dynamics, 
Volume I, Chapter IV. Since that time, a considerable 
number of papers about the same problem have been 
published by Schlichting, Falkner, Piercy, Whitehead, 
Tyler, Thwaites, Trilling, and others.'~* Most of these 
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authors are using the momentum equation of the 
boundary layer. Trilling solved recently the equation 
in Crocco’s form, though in a somewhat complicated 
manner, by a series of powers of the velocity u with co- 
efficients that are functions of x. 


The present paper starts from the boundary-layer 
equations in the original form (Prandtl) and assumes an 
approximate profile of the type 


UG = ety + b(x)y? + c(x)y3 + d(x) y*) 


This choice of the velocity profile is suggested by the 
following general consideration. It can be assumed that 
a solution u of the boundary-layer equations can be ex- 


panded in the form 
u = U| 1 — exp ( > an(v)y") 
n=1 


at least within a neighborhood of y = 0 for any x. Re- 
placing now the power series in y by a polynomial in y, 
satisfying the boundary-layer equations near y = 0) as 
far as possible for the chosen degree of the polynomial, 
one obtains an approximation near y = 0. However, 
one can expect to obtain in this way, because of the 
general simplicity of the shape of a velocity profile, 
even a good approximation for any value y > 0 if the 
coefficient of the highest power in y proves to be nega- 
tive, because u approaches then the asymptotic value 
U(x) in an exponential character as it must be expected 
out of physical reasons and as it is suggested also by the 
few known correct solutions. Using as degree of the 
polynomial in y the value 4, one obtains a high ac- 
curacy in all well-known examples of the boundary- 
layer theory not only near y = O but for every y > 0. 
Another reason for the choice of this velocity profile was 
the fact that it does not offer particular mathematical 
difficulties for the computation as normally connected 
with suction problems especially. Therefore, this 
method seems to be well suited for theoretical and 
practical use, combining high accuracy with conceptional 
simplicity. 


(1) DETERMINATION OF THE COEFFICIENTS OF THE 
ASSUMED VELOCITY PROFILE 


The boundary-layer equations in the original form 
are 
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LAMINAR BOUNDARY 
Ou Ou p'(x) O7u 
u +v = , 
Ox Oy p oy" (1) 


(Ou/Ox) + (Ov/Oy) = O 


where u, v are the velocity components of the flow at the 
point x, y within the boundary layer. p(x) is the given 
pressure distribution of the flow outside the boundary 
layer, p is the mass density, and » is the kinematic 
viscosity of the incompressible fluid. p’(x) means the 
derivative of p(x) with respect to x. 
Outside the boundary layer, Bernoulli's equation 
UU’ + (p’/p) = 0 (2) 
is valid where L’ is the flow velocity at the border of the 
boundary layer. If the suction (blowing) velocity 
along the surface y = 0 is denoted by %(x), the second 


of Eqs. (1) yields 


Y Ou 
v = u(x) — dy 
J 0 Ox 


and the elimination of p and v from the Eqs. (1) and (2) 
leads, therefore, to the condition 


Ou 4 ( (x) “y Ou ; >" yur Ou . 
— tt / ax > Jay "ay? 


(3) 


for u, which may be denoted as the boundary-layer 
equation. 

With p(x) is U(x), a given function. 
solution of this equation with the form 


An approximate 


om U1 — 7 toh tr ter (4) 
shall now be determined. Here are a(x), d(x), c(x), d(x) 
unknown functions of x which are supposed to have 
continuous derivatives of the first order. It is further 
supposed that U’ has a continuous derivative of the 
second order* and wv a continuous derivative of the 
first order. 

The boundary condition u = 0 for y = 0 is always 
satisfied by the assumed function [Eq. (4)]. The 
boundary condition « = U for y = © is satisfied if 
d < Oor, generally, if the coefficient of the lowest power 
of y which does not vanish is negative. It is then, also, 


every derivative 

(0"u/dy"),-0 = 0 
and, further, 

(0u/Ox),-9 = O 


The functions a, b, c, d can be determined by satisfy- 
ing the boundary-layer equation [Eq. (3)] along y = 0 
as far as possible. The best way to do this is to dif- 
ferentiate Eq. (3) three times with respect to y. One 
obtains then, for y = 0, inclusive { Eq. (3) |, the following 


four conditions: 


* This supposition is necessary for the theoretical derivation 
only. The final result will not contain the second derivative. 
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Ou ae O7u 
Vo —- l l —_= / * =() 
oy oy? 
O*u O*u 
Vo aes v - = 0 
Oy" oy* 
Ou Ou O*u O'u 
Vo =, = 
Oy Oxdy - Oy* oy! 
Ou O*u O*u Ou 
z rs Vo [= | - = 0 
Ov Oxdy" Oy' oy? 


The derivatives that appear in these equations have to 
be calculated from Eq. (4) for y = 0 and substituted 
here. Thus, the following four conditions for a, b, c, d 
are obtained: 


a + 2b = (1 V) (avo + U”) 


a* + 6ab + 6c = (v/v?)(avy + U’) 


a‘ + 12a*b + 12b? + 24ac + 24d = 
— (1/v)(U’a? + Uaa’) + (x?/v*)(avo + U") 
2a[L’’(a? + 2b) + U(2b’ + 2aa’)] — u(at + 12a*b + 
12b? + 24ac + 24d) = vi[a(a*t + 12a*b + 
12b2 + 24ac + 24d) + Sb(a* + Gab + 6c) + 
36c(a? + 2b) + 96ad] 


The first three of these equations determine the func- 
tions }, c, and d as functions of a: 


2b = (1/v)(am + U’) — a? | 


v 3a oj ; 
6c = & —_ ) (avy + U’) + 2a’ | 


y? v 
(U’a? + Uaa'’) — 6a* + 
y 


) 


Uo" Vo a’ , 0” rI\ 9 
(* — da = + 12 ) (avo + l | - (avo + l )- 
py Vv v* 


If these values for b, c, d are substituted into the 
fourth equation, the following condition for a is found: 


9 
Zoey 


10v» a ee oma a“VUp 
— — 28 (avo + U’)? + [ 55a* — 21 
y= vy Vv 


) 


a aU" U° — 2 f d - 
5— — —] Qn t U’) + -al (avy + U’) + 
py? py? v dx 
Uo - aoe ; ao is 
— 5a) Uaa’ —24va®* + d5a4tvyg — a® = QQ (6) 
Vv v 


Thus we have the results: The coefficient a(x) is 
connected with the outside velocity U(x) of the flow and the 
suction velocity v(x) by the relation (6). The coefficients 
h(x), c(x), and d(x) are functions of a(x), U(x), and v(x 
according to fq. (5). 

Eq. (6) is an ordinary differential equation of the first 
order for a if U and vw are given functions. It is also an 
ordinary differential equation of the first order for a if 
U and a are considered as given functions. 

Now the shear along the boundary y = O is given by 


tT) = vp(Ou/Oy),-9 = —vpla (7) 
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Thus, if U(x) and the shear 7o(x) are given functions, 
a(x) is determined by Eq. (7) and the ordinary differen- 
tial equation of first order [Eq. (6)] for v solves ap- 
proximately the problem of the determination of the 
suction (or blowing) velocity that is required in order to 
produce a given shear 79(x). 

A dimensionless form of Eq. (6) can be obtained re- 
placing the variables x, v1, U, and a by the dimension- 


less variables 


Uo / va ] 


(1/lva?) (dU /dE) § 


g = x/I, t 
U/lva? ¢ 


} (8) 


oOo = 
where / is any significant length. The differential equa- 


tion [Eq. (6)| then becomes 


i 
26 — (+f) = elt) + et) + walt)? + 


dé 
449 i Sees do. . 
o3(C)t? + 26 4 4) \ae é (9) 
with 
golf) = 2867 — 55¢ 24 
gilt) = —10f? + 775 — 60 (10) 
g(f) = —25¢ + 50 | 
g(t) = ¢ — 15 


However, this general dimensionless representation is 
not always convenient. Different dimensionless repre- 
sentations will be applied, therefore, in the following 


from case to case. 


(2) THe LAMINAR BOUNDARY LAYER WITHOUT 
SUCTION 


Without suction follows from Eq. (6) because of 
YU = O 
Sa ane Ti 5 __ da’ 

2+ 55U'a? — 24va! — — l = 0 
y v 2 as 


(11) 


Therefore, a? is a solution of this Riccati’s differential 
equation if there is no suction. a* being determined, the 


coefficients 5, c, d are obtained from 


2b = U'/v — a’ 
6c = —3(U'/v)a + 2a 
(12) 
2 UU" 13 /U'’\? 6 
iio « — —— 4 — a gil 
” ad 2) Vv 2) 
according to Eq. (5). 
For the flat plate without suction (U = const. = Uo) 


Eq. (11) reduces to 
(5/2) Uo(da?/dx) + 24va* = 0 
which has the general solution 


l a’ = (48 Dd) (vx Uo) + Cc 


where C is a constant of integration. From Eq. (12) 


there follows 
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— a*/2, C= 6/3, d = —a‘*/20 


b= 


The velocity profile {Eq. (4)] of the boundary layer has, 
therefore, the form 


(1/2) (ay)? + 1/3) (ay)? — (1/20) (ay)4 


u = U(1 — e” ~ ) 
If now x = Ois the beginning of the flat plate, the initia] 
condition u — U> for x > 0 and every y > 0 shows that 
it has to be C = 0. Thus, 


a= am ¥/ (5 48) (U5 Vx) 


The shear along the surface of the plate follows from 
Eq. (7) and has the value 
> 0.323 Uo? Vy Uox 
Blasius’ correct solution yields the numerical constant 
0.332, the error therefore being 2.7 per cent. 
Let us apply Eq. (11) to the boundary layer of the 
stagnation point where 


U = Uj(x/l) 
Eq. (11) then gives 
da* 56 Uo 4. 9942 1S at 
=— 22a? — 3 
dinx 5 5 U/W = 


Approaching the stagnation point with « — 0, a? would 
become infinite if not da?/dx —~ 0. From Eqs. (12) and 
(4), it would follow that u = 0 for x = Oand y > 0, 
which would be a contradiction to the velocity distribu- 
tion of a stagnation point. Only a constant solution a’ 
of Eq. (13) therefore comes into question. The two con- 


stant solutions of this equation are 


(ly/Up)a? = 1.529, (lv/Us)a? = 0.763 


Now there has to be, however, 


13 (/Uo\? 6 
24d = — — -a'<0 
o ly » 
Thus only 
a = —1.237 V U>/lv (14 
The shear 
«tone IUul (‘) 
7 er y \ 


The correct value of the constant calculated by Hie- 
menz" is 1.233, the error of the approximation being 
With a, the coefficients b, c, d of the 
From Eq. (12) it 


0.52 per cent. 
velocity profile are also constants. 
follows, with Eq. (14), that 


b = — 0.264( Uo ly) 
c = —0.0115(U)/ly) ”’ (15 
d = —0.00855( U//v)? 


Setting now 


9 = yV U/W 
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the velocity profile of the boundary layer becomes 


1.2377 0.2642 0.01153 0.008554 
” ” ” 


—@ 


I| 


u/U(x/1) 


This profile is plotted in Fig. | as solid curve. Values 
that Hiemenz obtained by exact computation are indi- 


cated by small circles. 





The general equation [Eq. (11)] can be simplified 


considerably. It can be transformed especially into 
such a type that the second derivative of U’ disappears, 
which is convenient for practical applications. This can 
be done replacing a® by the function \ according to 


a? = (d/v) + (4/5)(U’/v) (16) 


Eq. (11) then yields, for A, the condition 
UN = 0.256U’? + 6.64U’’ — 9.6? 


By solving this equation, a? and the shear 7» [from Eq. 
(7)| are computable without knowledge of U”. The 
last coefficient d of the velocity profile, however, re- 
quires L’", as follows, from Eq. (12). 

This Riccati’s equation for \ can be transformed 
further by the usual substitution 


A = (U/9.6)(d In 2/dx) (17) 


into a linear differential equation of the second order. 
Replacing the differentiation with respect to x by a dif- 
ferentiation with respect to x//, and U’ by U/W», this 
differential equation gets the dimensionless form 


din(U/U,) dz 
d(x//) d(x/1) 


ne eee. = 
2.4576 z:=0 
L d(x 1) 


A solution z of this equation yields, because of Eqs. (16) 


dz 
d(x/l)2 


(18) 


and (17), 


lva? nS ding 4 din (U/Us) - 
U  96d(x/l) | 5 d(x/l) aia 
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If especially 
U/Uo = [a + B(x/l)]™ 


(a, 8, m are constants), the general solution of Eq. (18) 


becomes 


s= K, (« + 3*) + Ke (« + a*) | (20) 


A + 1\° 
( a“ ) + Bm? 


where 


and 


A = B = 2.4576 


5.64, 


K,, AK» are arbitrary constants. For /ya*/U’, their ratio 
only is important and this has to be determined from 
the initial condition of the particular problem. 


In the special case 
U £P = e* + B(x/l 
the general solution of Eq. (18) is 


= Kye" + Koe™%" (21) 


with 


Ni,2 = (8/2)(A + V A? + 4B) 


Both cases are convenient for the practical computa- 
tion of the boundary layer if these two velocity func- 
tions are used for a stepwise approximation of any 
given velocity distribution. Both cases are character- 
ized by the fact that the velocity profiles of boundary 
layer become similar (Thwaites’). 

If, finally, x// in Eq. (18) is replaced by 

E=k thi f (U/Us)**4d(x/1) 


with two arbitrary constants ky and ), this equation 


a d In (U Up) 2 
= 2.4576 je Z 
dé 


Thus, the investigation of the laminar boundary layer 
without suction is reduced to a differential equation of the 


gets the form 
d*z 


dé? 


type 


(3) APPROXIMATION OF THE VELOCITY PROFILE BY A 
DISPLACEMENT AND MOMENTUM 
THICKNESS 


GAUSS-CURVE; 


The main purpose of the following approximation of 
the velocity profile by a Gauss-curve is to simplify the 
calculation of the displacement and the momentum 
thickness of the boundary layer. 

If one tries to solve the boundary-layer equations in 
the same way as before by the simpler expression 
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ay + BM) (22) 


u= U(1l—e 


instead of Eq. (4), the accuracy of the result becomes in- 
sufficient. This is not too surprising because the deter- 
mination of the velocity profile by the given method 
starts essentially from conditions along y = O only. 
An accurate result can be expected, therefore, only if 
these conditions are sufficiently strong. This was the 
case with Eq. (4). However, it is possible now to deter- 
mine a sufficiently accurate Gauss-curve [| Eq. (22) | using 
the results obtained for the exponential solution with 
four coefficients. 

To this purpose, the fourth-order solution |Eq. (4) | 
shall be replaced by a Gauss-Curve [Eq. (22)] which 
has the same tangent as Eq. (4) for y = 0 and which 
passes through a certain point of Eq. (4) inside the 
boundary layer. If 

| eo + by® + cy* + dy* 
it is (OH/Oy),~9 = a and the tangent of this curve E(y) 
at y = 0 passes through the point y, = —1/a at the y- 
axis. The point P; of E(y) which belongs to this y; has 
the £-value 
Lag Sr see eee 
If this point P; of E(y) is chosen as the point through 


which the curve 


has to pass, one obtains 
B= b — (c/a) + (d/a?) (24) 


The solution (4) with the coefficients a, b, c, d determined 
by Eqs. (6) and (7) is in this way approximated by the 
Gauss-curve 

u = U(1 — et) 
with B given by Eq. (24). 

In the case of the stagnation point, because of Eqs. 

(14) and (15), 


B = —0.2789(U)/lv) 


If the corresponding Gauss-curve would be plotted in 
Fig. 1, there would be no perceptible difference with the 
already plotted curve | Eq. (4) J. 

Since a, b, c, d have to be found anyway, this simpli- 
fication to a Gauss-curve as such would not be essential. 
However, so it is for the determination of the displace- 
ment and the momentum thickness of the boundary layer. 

The displacement thickness 6*(x) is defined by 


si u 
6*(x) = | ( — r) dy (25) 
0 
’ u\ Uu 
6*(x) = f (: _ py dy (26) 
0 


With 
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é*(x) = / Edy (27 
/0 
6*(x) = 6*(x) — Z FE? dy (28 
0 


If £ is replaced now by the Gauss-curve [Eq. (23) 
the integration can be performed with the aid of Gauss 


») ae: ; 
e” @s (29 
Vit /o : 
and one obtains 


/ a*/4g 
é a 
é* (x) = == 1 — 6(- ) (30 
2 v-B 2-8 


vo 


The integral | i? dyin Eq. (28) follows from the same 
0 


error integral 


formula with 2a, 2b, 2c, 2d instead of a, b,c, d. If we de- 
note, therefore, 
B = 2b — (c/a) + (d/2a?) (31) 


it is 


2 x ‘ a* B 

? l 

| E? dy = vrs - [ = @ (- | )| (32 
0 2 vV-8 Vv ~-i 


The case of the flat plate without suction becomes 
B= —O.SS83a" 


with a = —0.323 V Uo/vx, and, because of Eq. (30), 
a” = 1.75 V ve/ Up. 


is 1.72. For the same case, 8 = —1.35Sa” and, because 


The correct value of the constant 


of Eqs. (32) and (28), 6* = 0.641 V vx/ Uo, the correct 
value (Blasius) of the constant being 0.664. 
(4) THe Fiat PLATE WITH CONSTANT SUCTION AND 
THE ASYMPTOTIC SUCTION PROFILE 
The flat plate with constant suction was computed by 
Iglisch'! in 1944. In order to treat this problem with 


the present method, the differential equation {Eq. (6) 
for a has to be solved under the conditions 


U = UU, = const., Vo = const. < 0 


which reduces this equation to 


. Vo = , ™ 2 
Up (: — sa) a’ = 24va*! — 6O0upa*? + 


_. Uo" Uo° Vo 
00 : 
v p~ is 


It can be solved by direct integration. 
dimensionless variables 


E = (u7/vUp)x, t = %/ra (34) 


one obtains 


Introducing the 
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LAMINAR BOUNDARY 


Flat plate with constant suction 


40\-— { T 





Approximation 








r) 10 ; : 45 
Ve j/ vx 
Fig. & “VI si 
dé 5t — 3f? 
= aes = (35) 
dt 24 — 60¢ + 500 — 15f7 + #4 - 
Here also, because of Eq. (7), 
t= — pU vo To (36) 


If the flat plate extends along the positive x-axis and 
the fluid flows in the direction of increasing x, the shear 
~ has to be infinite at x = 0, and from Eq. (33) it 


follows that 


9 bt — 3f 
é = / a a dt 
J, 24 — 0+ DP — ir + Ff 

or 
t = 0.4642 — 21n (1 — ¢) + 0.25 In (2 — #t) — 

0.351 In (10.899 — ¢t) + 2.101 In (1.101 — ¢) (37) 
The coefficients of the velocity profile of the boundary 
layer as functions of ¢ are 


= (Uo v) (1 t) 


a 
b= (*) () 7 ) 
2 \v t iy 
1 /w\?/1 3, 2 
ni ( ) ( ee :) > (38) 
6 \v t - g° | 
1] 21 18 6 
—2+ a a a 
~)' t t? Pe # 
d = a a 
24 \y 5 — 3 J 
The results following from these formulas have been 
plotted in Fig. For various values of t = —pU vo/ 70, 


= (uv?/vUo)x can be calculated directly 
from Eq. (37). ¢ is plotted versus V/¢ in Fig. 3. For 
the same values of ¢, the coefficients of the velocity pro- 
file follow from Eqs. (38). With a, 6, c, d is 8 known 
from Eq. (24) and, therefore, also 6* as function of ¢ and 
of VY With a, b, c,d is also known B from Eq. (31) 
and, therefore, 6* from Eqs. (32) and (28). In Fig. 2 are 
plotted —(v/v) 6* and —(u/v)0* versus Vé The 
small circles in this figure indicate the values calculated 
The accordance with the results of the 


2 
the coordinate 


by Iglisch. 


L 


AYER WITH SUCTION 5 
present theory is excellent, and the simplicity of the 
necessary numerical calculation involved by this method 
is obvious. 

Ihe preceding example includes the asymptotic suc 
tion profile. With — + is ¢ — 1, and therefore, be 
cause of Eqs. (38), 


b=cxd=Q0 


a = Uo/),; 


Thus, indeed, 


(v0/¥) ¥ 


u/Up = 1-—e 
With ¢ — | there follows, in this special case, 
—pU w/o = | 
A generalization can be obtained if one asks for the con- 
ditions under which, generally, 
—pUy/m = 1 (39) 


Because of 7) = —vpUa, the 


The general differen- 


for variable U and w. 

equation has then to be a = w/v. 

tial equation [Eq. (6)]| then yields 
UU" — 9U"? + (u?/rv) U’ = 0 (40) 


and the coefficients of the velocity profile become 


a = %/v 
b = (1/2)(U’/r) 
c = —(1/3)(v/v?) UL’ (41) 
few (- heey. = v) 
24 y? y* y® 


For any given not constant L’ exists a suction velocity 
v determined by Eq. (40), for which Eq. (39) holds 
under the condition that d < 0. If 


U = UU, = const. 


Eq. (40) is satisfied for every v%, and from Eq. (41) there 


follows 
Uo Uo dv? 
a= b=c=0,,d= —- = 
v {Xv* dx 
Thus we obtain 
»(v0/v) ¥ (U0 /48v*) (dv02/dx) y4 ( $2) 


u Uo = !—¢ 
as a suction profile of the flat plate for which 7 = 
—plUov, under the condition only that 


dvy” dx > 0 


This condition includes that the result has a meaning 
only for an infinite flat plate in both directions with ar- 
bitrarily increasing suction velocity in the flow direc- 
tion. The asymptotic suction profile appears again for 
Yo = const. 
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(5) Exact SOLUTIONS OF Eg. (9) CONSTANT 
Q 


SHEAR 


WITH 


The general equation in the dimensionless form [Eq. 


(9)] can be integrated exactly if ¢ and do/dé are con- 
stants. In this case there follows either 
U/Uo = (a& + B)”\ (43) 
ao = y(aé + B) f ; 
or 
U U - et + B 
c (44) 
c= y¥Y f 
where m, a, 8, and y are arbitrary constants. Accord- 
ingly, 
™m™ = pV ly(aé 7s By)? 2)m (1/2) (45) 
or 
Tm = pV v/ly eS@/2at + 6) (46) 


In the special case m = 1/3, the first 7) is constant. 


If, moreover, 


¢ = may = 5/4 (47) 
(9) has the solution ¢ = 0, 
which means that no suction is necessary. With a = |, 


8 = 0, it follows that for the velocity distribution 


and therefore ay = 15/4, Eq. 


q? Uo = gt (48) 
outside the boundary layer the shear along the boundary 


2 vUo | Ul 
o/is” l y 


% = (49) 


ts constant without suction. This special case corresponds 
to the symmetrical flow against a corner with plane sur- 
faces under 90° angle. 

The other special value a = 0 for which 7» is constant 
leads in both cases to 

U = Uo = const., o = op = const. 

and represents, therefore, the flat plate in longitudinal 
flow with constant shear due to suction. Eq. (9) yields in 
this case 


Sek ee dt 
aw ” te 24 — 68 + De — 


if suction is applied for —§ => &. 


Bf? (50) 


If ro is the shear with- 
out suction at the point & and at the same time the pre- 
scribed constant value of the shear for — > £, it follows 
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{= —pl ‘oo TOO 
é =~ £ Too U(x aia = \ol 
ou pU,? y 


because of Eq. (8). In Fig. 3, ¢ is plotted versus 


(E — &)/oo. This result shows that the suction ve- 
locity necessary to keep the shear constant increases 
rapidly from zero to a constant value corresponding to 


the asymptotic suction. 
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Tur ‘bulent Boundary Layer on a Cone in a 
Supersonic Flow at Zero Angle of Attack 


E.R. VAN DRIEST* 


North American Aviation, Ine. 


ABSTRACT 


A simple rule is given for transforming local heat-transfer 
coeficients from flat plates to cones at zero angle of attack for 
fully turbulent boundary layers in supersonic flow. The 
is that the cone heat-transfer coefficient is the flat-plate coefficient 
for one-half the Reynolds Number the the Mach 


Number and wall-to-free-stream temperature ratio remaining the 


rule 


on cone, 


same. 


INTRODUCTION AND SUMMARY 


_ CALCULATION OF THE TEMPERATURE of a cone 
in supersonic flight requires a knowledge of the 
heat-transfer coefficient of the boundary layer of the 
cone. Hantzsche and Wendt! solved the problem for 
laminar boundary layers. They found that the local 
laminar heat-transfer coefficient on a cone at zero angle 
of attack can be obtained by merely multiplying by 
V3 the local laminar heat-transfer coefficient 
flat plate. Since heat transfer and skin friction 
proportional, this rule applies also to local skin-friction 
Because laminar heat transfer 


on a 
are 


coefficients for cones. 
and skin friction vary inversely with the square root 
of the Reynolds Number, another way of formulating 
the transformation rule for conical laminar layers is 
to state that the cone solution is the flat-plate solu- 
the Reynolds Number on the 


tion for one-third of 


cone. 


The present report offers a simple transformation 
rule for fully turbulent boundary layers on cones, 
similar to that for laminar layers. It is found that, 
for turbulent boundary layers, the cone solution for 
local heat transfer is the flat-plate solution for one-half 
of the Reynolds Number on the cone, the Mach 
Number wall-to-free-stream temperature ratio 
remaining the same. The effect is an increase over the 
local turbulent heat-transfer or skin-friction coefficient 
of flat plates of about 10 to 15 per cent. The work 
is an extension to compressible flow of an analysis by 
Gazley,? who neglected the effect of compressibility 
and assumed a power law for velocity distribution; 
his results approach those obtained when the above 
general rule for turbulent boundary layers is ap- 
plied to incompressible flow at high Reynolds Num- 


and 


bers. 
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BOUNDARY LAYER ON 
ATTACK 


MOMENTUM EQUATION FOR THE 
A CONE AT ZERO ANGLE OF 


The von Karman momentum equation for a boundary 
layer in steady state on a body of revolution at zero 


angle of attack is 


a "5 . Op 6 
rpu(u.. — u)dy = — rdy — Put» (1) 
Ox 0 Ox 0 . 


is a coordinate distance measured 
along the body from the forward most point, y is the 
other coordinate measured from the surface along 

u is the flow velocity relative 
, p is the 


In this relation, x 


normal to the surface, 
and parallel to the surface at the point (x, 
density of the fluid at (x, y), 7 is the normal distance 
from the body axis to (x, y), p is the pressure at (x), 
T» is the shear stress at the surface, and 6 is the bound- 
ary-layer thickness measured normal to the surface. 
Subscripts w and © refer to the wall and outer edge of 
the boundary layer, respectively. 

With a cone, when the flow is supersonic and the 
shock wave is attached, the pressure along the rays 
is constant for inviscid fluids. Therefore, since the 
boundary layer is thin (6<vr), it is assumed that the 
pressure is also constant along the boundary layer 
reduces to 


i.e., OpP/Ox = 0. Hence, Eq. (1) for cones 
ra) 6 
rpu(u. — u) dy = —Puyty (2) 
Ox 0 
Furthermore, as a consequence of 6<r, Eq. (2) be- 
comes 
) ” 
Te pu(u. — ua) dy = —Pyty (3) 
Ox Al 
or 


o = 
/ pu(u.. — u) dy + 


Ox 
in f° 
pu(u.. —u)dy=—Ty (4) 
Ve Ox J % 
(5) 


But for the cone 
‘t0) (Ofw/Ox) = 1/x 


so that, finally, from Eq. (4), the momentum equation 
for thin boundary layers on a cone in supersonic flight 


at zero angle of attack is 
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: f ( )dy + 
u(u. — u) dy 
Ox 0 P . 
| 6 
f pu(u.. —u)dy (6) 
x 0 


For turbulent flow, the fluid properties, p and wu, in 
this equation take on mean values. 
HEAT 


FRICTION AND TRANSFER 


ON A CONE 


TURBULENT SKIN 


It is assumed that the same flow mechanism holds 
locally on a cone as on a flat plate. Hence, the same 
procedure is followed for skin friction on cones as was 
used in the writer’s previous report*® on turbulent 
boundary layers on flat plates. According to that 
procedure, it was at first shown that the apparent 
shear stress was given by 


r = —(pv)'u’ (7) 


whence, when @ ~ O for thin boundary layers and the 


, 


triple correlation p’u’v’ is neglected, 


T= —pu'v’ (8) 


where p is the fluid density and uw and v are the velocity 
components parallel to the free stream and normal to 


the plate, respectively. The primes indicate fluctua- 


tion from the mean. Introduction of the Prandtl 
mixing length / then gave 
tr = pl*(du/dy)? (9) 


The shear was next assumed constant and equal to the 


wall value, and / was taken as Ky. The mean local 


density was related to the mean local velocity by using 
r Te (> ) ul 4 
le Tw ie tics 


— ] Uu u 
——— (1 - ) (10) 
UU « U « 


bo 


in 
P/ Pu T/T (11) 
to obtain 
a 1 +B - a:(*)'| (12) 
Pw U. Un 
where 
A? = [(y — 1)/2]M.?/(To/T ) 
and 
B= ({1 + [(y — 1)/2]M.2}/(To/Ts)) — 1 


and where 7° is temperature, y is the ratio of the specific 
heats, M is the Mach Number, and w and = refer to 
wall- and free-stream conditions, respectively. Inte 


gration of Eq. (9) then yielded 
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- 2A?(u/u.) — B 1 . B 
sin! z a ae eo ieee = 
A (B? + 4A?) ”* A (B? + 4A?*)” 


l he, . l tj. ¥ 
(r+ In . *) (13 
Ue Vv... K a Was 
where F is a constant and » is the fluid kinematic vis. 
cosity. 

With flat Karman 
relation constitutes only the first term on the right-hand 
Substitution of Eqs. (12) and (13 
in that term then gave 


plates, the von momentum 


side of Eq. (6). 


Divttna Af. 5 , B | 
Tr = as { Ja? exp a a 
kK dx \ A (B2 + 442)'7"]{ 
(14 
where 
D=e** 
a= Ku./V te/ feo 
1 ° —- 
fan i |. 
o (1+ Bz — A?Y?)” 
a. 2A*z — B 
exp a> ——- — 7, ae 
A (B? + 4A?)” 
2 = Uu/u. 


In these equations, coordinate x is the distance 
measured from the leading edge of the plate in the 
free-stream direction. The integral J was expanded 
in a series and terms of higher order than 1/a? were 


neglected, since a was large, resulting in 
| a . 2A? — B 
iow >. ep an? — 
a?(1+ B— A?) A a?(B2 + 4A?)” 
(15 
After rearrangement, Eq. (14) became 


Pull « D l 


= Ht . i). . a’ . 
Mu K* 


(1 + B — A?) 


d f | a , 8) } 
ay ‘ 1 oa Ss 1 . ( h) 
lx \S xp (sin a in ( 


where 
a = (2A? — B)/(B? + 4A?)' 
and 
8 = B/(B? + 4A?)” 
The lower limit on a for x = 0 was taken as zero, since 


the shear at the leading edge of the plate was assumed 
to be infinite. Integration of Eq. (16) then yielded, 
for large values of a and negligible variation of wall 
temperature with x, 

Pull aX D l 


= ° ws 6g? 


Mn K* (1+ B — A?)” 


a 
exp F (sin—! a@ + sin a)| (17 
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from which came the form of the local skin-friction 
(heat-transfer) formula 
1 


9 2 
“s 1a + sin 


~! B) = const. + 


ere 
Atjo (7 a/7 oo) : 


l 1 us 
K [in Rojo — ( ao s) In =| (18) 


since 


Po = po(T, wo/ 5 


te = te(T./T a)” 

Cre = 20/pata® 
ee, Se 

- « RATT.) °* 
Ro = PaotkaX/ Mba 


A similar development can be traced for the cone in supersonic flow with attached shock wave. 


Eq. (16), one obtains, starting with Eq. (6), 
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Fic. 1. Method of obtaining the local cone heat-transfer 


coefficient from the flat-plate solution. 


Instead of 


rua D I “ [<4 |S i ieee ‘alba lf * (si iives ) 
= . =e 9 Jexp sin sin exp sin sin 
Bw K* (1 + B — A?)” dx \ ; A = x | ” f 
(19) 
Now, with the assumption that a is large, Eq. (19) can be written as 
Pull cx D l E os 5 ee | 4 e* |; 5 Glee + eb 
= . — 4a? exp sin sin + ex sin S 
nw K* (@+B—A)” Lax 0 “PLA - i *se"um” «= =e 
(20) 
Rearrangement gives 
Proll D I - |: ; \ 
-xdx = — - 7, + d \xa* exp (sin~! sin! B) |> (21 
bu ae eS ae ee es Yh — | 
which integrates to 
] Pull oX D l - : ( ss , 4 ‘ _ | (29 
, a — & 7, ° X%* 2° Gap sin sin ) 22) 
rs K* (1+ B-— A?) viewed - ' 
or 
1 pu lheoX D l : 5 ie aie a)| to: 
. = . 7. * a ex sin sin 2: 
[i «= © 6te ee oa 


Comparison of Eqs. (23) and (17) shows that, for 
the same cone Mach Number and wall-to-free-stream 
temperature ratio, the cone solution for local skin 
friction (and therefore heat-transfer coefficient, since 
they are proportional) is the same as the flat-plate 
solution when the cone Reynolds Number is divided 
by 2. The use of this rule is shown in Fig. 1 in which 
the local heat-transfer coefficient is plotted as a func- 
tion of Reynolds Number and temperature ratio at a 
Mach Number of 5 for flat plates. For example, if the 


5, the wall-to-free- 
Reynolds 


cone free-stream Mach Number is 
stream temperature ratio 4, and the cone 
Number 10’, then the cone local heat-transfer coefficient 
is found in Fig. 1 on the temperature ratio 4 curve at 
Reynolds Number 5 X 10°. Thus, in this particular 
case, it is found that the heat-transfer coefficient for 
the cone is 1.140 times that of a flat plate under the 
same conditions. This ratio will vary with all of the 
above flow parameters. 


(Continued on page 72) 








Column Behavior Under Conditions of Impact 


GEORGE GERARD* ann HERBERT BECKER? 


New York Unwersity and Combustion I 


SUMMARY 


The magnitude of the stress introduced into a long slender bar 
by velocity impact depends only upon the conditions of impact. 
Thus, buckling of this bar under impact loading can only take 
place after the compressive stress wave has propagated over a 
critical’ length of the bar. Column behavior under impact is 
studied theoretically and by a novel experimental technique, and 
it isconcluded that a column can momentarily support a dynamic 
compressive stress of any magnitude which may be introduced by 


impact. 


INTRODUCTION 


HE INVESTIGATION OF COLUMN BEHAVIOR. under 
j paeeae of dynamic loading has been studied by 
several investigators. Meier! has studied the behavior 
of columns containing initial imperfections by investi- 
gating the motion of the center of the column for loads 
below the Euler load Although Meier 
speculates that the column can carry loads in excess of 
the Euler load, a study of his analytical results fails to 


and above. 


bear out this conclusion. 

Hoff? has studied a similar problem but has specified 
the displacement at the end of the column rather than 
the load. His investigation is concerned with an 
initially imperfect column in which one end of the 
column has been displaced at a constant velocity. 
Under such conditions it has been shown analytically 
that the column can carry loads in excess of the Euler 
load because of the lateral inertia of the imperfect 
column. 

In both the above investigations, it has been assumed 
that the period of the lateral motion of the column is 
small relative to the time required for a stress wave to 
propagate from one end of the bar to the other. Thus, 
it has been assumed that the stress is constant through- 
out the entire length of the column. In contrast, the 
investigation contained herein is concerned with the 
propagation of an elastic stress discontinuity in the 
column and accounts for the fact that the stress is in- 
stantaneously not constant over the entire length of the 
column. Under conditions of compressive stress wave 
propagation, it has been found useful to adopt the point 
of view of determining the critical length of a column 
subject to a given stress rather than finding the critical 
stress of a column of given length. 

It is the objective of this investigation to show that 
under certain conditions of impact loading, the com- 
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pressive stress momentarily supported by a perfect 
column at buckling may be of any magnitude in excess of 
its static Euler load. Under sufficiently high velocities 
of impact, the critical length of the column is deter. 
mined by the magnitude of the stress being propa 
gated, which is contrary to the usual situation of static 
buckling in which the critical stress is determined by 
the effective length of the column. 
due to the fact that the magnitude of the stress at the 
end of the bar which is introduced by velocity impact 


This phenomenon is 


depends only upon the velocity of impact and _ the 
modulus and density of the material. Initially, the 
stress is only at the struck end. Since there is no 
stress in the remainder of the bar, it is only after the 
stress has been transmitted over a critical length of the 
bar that buckling can possibly occur. 


SYMBOLS 


A = area 


c = velocity of propagation 

E = elastic modulus 

I = moment of inertia 

I = length 

Le = length traversed by stress wave in time ¢, Ly = ct 
Ler. = critical length in which buckling occurs 
m = mass 

t = time 

U = elastic strain energy 

v = velocity 

V = volume, V = LA 

u = lateral displacement 

x = coordinate along length 

o = stress 

oy, = ultimate tensile strength 

€ = strain 

p = radius of gyration 

¢g = density 


TEST TECHNIQUE 


Fig. 1 is a photograph of a long thin column in which 
buckling was caused by the propagation of a compres- 
sive stress wave. The impact loading was achieved 


experimentally in a simple manner. The specimen was 
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loaded in tension until fracture occurred. At fracture, 


the elastic strain energy stored in the specimen was 


U = (0,7/2E)V (1) 
This energy is available as impact energy. Thus, 
mv?/2 = o,’m/2E¢ (2) 


Noting that c = (E/g)'* for longitudinal wave propaga- 


tion, the effective impact velocity obtained is given by* 


E)o, (3) 


v= (Cc 


Since for a given material the effective impact ve- 
locity is a linear function of the stress in the bar, a 
simple technique for obtaining various impact velocities 
is suggested. By using different size notches at some 
point in the bar, the stress in the unnotched portion of 
the bar can be varied as desired. While the average 
fracture stress at the notch would be o,, the stress in 
the unnotched portion would be lower. According to 
Eq. (3), this would result in a lower effective impact 
vel city : 

The compressive stress wave that buckled the column 
shown in Fig. | was a reflection of the unloading wave 
(negative tension) which followed fracture at the rela- 
tively massive grips that held the specimen. The wave 
propagation following fracture is depicted in Fig. 2. 
The velocity of the unloading wave is that associated 
with the elastic modulus of the material, since the un- 


* For a rigorous development of Eq. (3), see Timoshenko, &., 
Theory of Elasticity, 1st Ed., pp. 381-384, McGraw-Hill Book 
Company, Inc., New York. 
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loading path is elastic. Upon reflection, however, the 
stress may reach the plastic portion of the dynamic 
stress-strain curve, and therefore the plastic stresses 
will propagate at lower velocities than the elastic 
stresses. With this technique of compressive impact 
loading, the effective origin of the dynamic stress-strain 
curve is essentially shifted to the fracture point of the 


static tension stress-strain curve. 


COLUMN BEHAVIOR 


In its gross features, column behavior under com- 
pressive stress wave propagation can be represented by 
the model shown in Fig. 3. To the long slender bar 
fixed at the lower end, a given weight is assumed 
attached at successively increasing distances from the 
lower end. When the weight has moved a sufficient 
distance from the fixed end, ‘buckling of the bar will 
occur in the length between the lower end and the 
weight. Under the corresponding conditions of com- 
pressive stress wave propagation, the magnitude of the 
stress wave is fixed by impact conditions. As the wave 
front moves away from the fixed end, buckling will 
occur when the “‘critical’’ length corresponding to this 
stress magnitude has been traversed. 

In the following analysis, it is assumed that the 
stresses are elastic, and therefore the stress wave is 
propagated as a single discontinuity and is constant 
over the length traversed by the wave. The equilib- 
rium differential equation for the column is 


EI0*w/0x' + cAd*w/dox* = 0 (4) 
This equation can be rewritten in the following form: 
p°0'w/Oxt + €0°w/Ox? = 0 (5) 


The assumed boundary conditions are clamping at 
the fixed end of the bar (x = 0) from which point the 
compressive stress wave originates. The additional 
boundary conditions must be prescribed not at the 
other end of the bar but at the stress wave front where 
the first influence of the disturbance is felt. Since there 
can be no lateral deflection in front of the stress wave, 
continuity conditions prescribed the boundary condi- 


tions at the wave front. Thus, 


w(0) = Ow/dx(0) = Ol 6) 
. ) 
w(L,) = Ow/dx(L,) = OF 
where L, = ct and is the distance the stress wave has 


traveled along the bar. A solution that satisfies Eq. (5), 
as well as the boundary conditions of Eqs. (6), is given 
by 


w = all — (cos2rx/L,)] forO <x Sk (7) 


By substituting the appropriate derivatives of Eq. (7) 
into Eq. (5), 


[16p2(4/L,)* — 4e(x/L,)*]a[1 — (cos 24x/L,)] = 0 (8) 








60 JOURNAL OF THE 
For a nontrivial solution, the first bracketed term must 


vanish. Consequently, 


Lee, = 2mpe~ (9) 


and the critical length of the bar in which buckling 
occurs depends upon the magnitude of the stress wave. 
If use is made of Eq. (3), then Eq. (9) becomes 


L.. /p = 2n(c/v)' (10) 
The dependence of the critical slenderness ratio upon 
the impact velocity given by Eq. (10) is shown in Fig. 4. 
In cases where v > c, the compressive stress wave would 
still be propagated at a velocity equal to c. Therefore, 
for L., /p < 2m, buckling can never occur. 

It is interesting to note that, if the usual method of 
running static column tests (in which a column of fixed 
slenderness ratio is subject to increasing stress until 
buckling occurs) is used for dynamic loading, then 
buckling always occurs over the entire length. of the 
column. This is due to the fact that the impact veloci- 
ties would be increased only until the proper stress 
magnitude for buckling in the entire length is reached 
and then not continued. If the buckling tests are run 
according to the model shown in Fig. 3, then buckling 
will take place in a critical length less than the actual 
length of the bar. 

The solution for the critical length given by Eq. (9) 
is based on the assumption that the stress wave is elas- 
tic. If the stress is of such a magnitude so as to be in the 
plastic region, then the plastic portions of the stress are 
propagated at lower velocities than the elastic portion 
as shown in Fig. 2. In this case, the elastic portion of 
the stress wave traverses the critical length of the bar 
first. Therefore, to a first approximation, it appears 
that the value of the critical length may be closely asso- 
ciated with the proportional limit stress, which in itself 
may be high because of the large strain rates associated 
The 


behavior of the column under conditions of plastic stress 


with the impact velocities under consideration. 


wave propagation requires detailed investigation, how- 
ever, and therefore the above remark is merely an 
intuitive suggestion. 


EXPERIMENTAL DATA 


Several column tests were run using the experimental 
technique outlined previously for obtaining compressive 
impact loadings of high velocity. A '/2-in. strip of 
0.010-in. 24S-T3 aluminum alloy of 10-in. length was 
held in a testing machine by Templin grips. The speci- 
mens were loaded to failure in tension, and it was ob- 
served (as shown in Fig. 1) that a series of short wave- 
length buckles formed immediately above the grip after 
reflection of the unloading wave. The wave length of 
the buckles of four similar specimens was approximately 
'/,in. The ultimate tensile strength of the specimens 
averaged 64,000 Ibs. per sq.in. It is to be noted that, in 
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stretching the bar to fracture, any initial imperfections 
in the column are effectively removed and the column 
as tested is essentially perfect. 

From Eq. (9), the estimated critical length based on 
the assumption of elastic buckling was 0.23 in., which 
certainly is of the order of magnitude of the wave 
lengths observed. 


CONCLUSIONS 


Both experiment and theory indicate that, under con- 
ditions of velocity impact in which the magnitude of the 
compressive stress wave may be far in excess of the load 
required to cause buckling of the column under static 
conditions, buckling may be confined to a length less 
than the actual length of the column. This phenomenon 
indicates that a column can momentarily support a 
dynamic compressive stress of any magnitude that may 
be introduced by impact. Therefore, dynamic buckling 
does not place an upper limit to the compressive stress 
that can be carried by a slender bar in the sense of 
static buckling. 


(Continued on page 65) 
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On Supersonic Flow of a Two-Dimensional 
Jet in Uniform Stream’ 


S. I. 


PAIt 


Unwersity of Maryland 


SUMMARY 


It is well known that a supersonic gas jet issuing from a reser- 
voir into a medium at rest has a periodic structure if the differ- 
ence of pressure of the jet from the medium is not large. This 
paper investigates a similar situation when a supersonic gas jet is 
issuing into a uniform stream. It is found that if the uniform 
stream is subsonic, the supersonic jet has a/most periodic struc- 
ture, and the approximate wave length increases with the Mach 
Number of the surrounding stream for a given Mach Number of 
If the uniform stream is supersonic, the supersonic jet 


the jet. 
The transmission and the 


does not have periodic structure. 
reflection of small disturbance at the boundary of the jet—i.e., 
Factors of transmission and 
They are 


a vortex sheet—are investigated. 
reflection of disturbances on the vortex sheet are found. 
functions of the Mach Numbers of the jet and that of uniform 
stream. The reflection wave may be of the same sign or of op- 
posite sign of the incident wave or zero according to the Mach 


Numbers of the two streams. 


List OF SYMBOLS 


A = constant 

B = constant 

f(z) = function of 

F(z = function of z 

g(z = function of 

M = Mach Number of basic flow 

R = reflection factor defined by Eq. (24) 

: = transmission factor defined by Eq. (23) 

u = x-wise velocity of component of the resultant 
flow 

U = velocity of the basic flow 

i = y-wise velocity component of the resultant flow 

= = distance along the axis of the jet 

y = distance perpendicular to the axis of the jet 

B, = VIM, —1| (1 = lor 2) 

6 = wave-length factor defined by Eq. (10) 

d = frequency or eigenvalue 

o = velocity potential of disturbance 

subscript 1 = values for quantity in the jet 


subscript 2 = values for quantity in the surrounding stream 


(I) INTRODUCTION 


N HIS EXPERIMENTAL INVESTIGATIONS Of a jet of gas 
I issuing from an orifice into the medium at rest, 
Emden! found that, if the velocity of the jet is super- 
sonic and the excess pressure of the jet over the sur- 
rounding medium is not large, the jet has periodic struc- 
ture. It was Prandtl? who first attempted to explain 
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this phenomenon from a theoretical point of view by the 
method of small perturbation and reasonable good agree- 
ment of the theory with experimental results was 
Several other authors had discussed this prob- 
However, the 


found. 
lem from a theoretical point of view.*~® 
behavior of the supersonic gas jet in a uniform stream 
has not yet been studied. This paper is to study the 
case of a supersonic gas jet in a uniform stream when the 
pressure of the jet differs slightly from that of surround- 
ing stream. Because the behaviors of the subsonic and 
the supersonic uniform stream are entirely different, we 
have to treat these two cases separately. 

Since the pressure difference between the jet and 
the uniform stream is not large, we may assume 
that the gas velocity normal to the jet axis is small—in 
other words, that the jet expands very little—hence, 
the method of small perturbation may be used in this 
problem. We also assume that there exist velocity po- 
tentials inside the jet, as well asin the uniform stream 
i.e., even if shock wave may happen in the jet, the 
strength of the shock is assumed to be small. Thus, the 
vorticity introduced is still negligible. Both viscosity 
and heat conduction are neglected. 

In this paper we shall consider only the two-dimen- 
sional flow. Here, besides the discussion of periodicity 
of the supersonic jet, some interesting phenomena of a 
small disturbance reflected from, and _ transmitted 
through, a surface of discontinuity have been obtained. 

The author also found that, as far as periodicity is 
concerned, the two-dimensional and the axially sym- 
metrical jets behave similarly. However, in case of 
reflection and transmission of disturbance on a surface 
of discontinuity, the axially symmetrical case is much 
more complicated. Hence, the results of the axially 
symmetrical supersonic jet in uniform stream will be 
reported in another paper elsewhere. 


(Il) FORMULATION OF THE PROBLEM 


The schematic diagram of the problem is shown in 
At the exit of the nozzle x = 0, the jet has the 


Fig. 1. 

following characteristics: 
u = U, + uly); v = wl(y)t (1) 
M= M, + Moy) { 


where uo < U1, Uo ~ v0; Myo < Mi, where Mis the Mach 
Number corresponding to U;, and U; is certain mean 
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Fic. 1. Supersonic jet in uniform stream. 
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velocity in the jet where the pressure in the jet is equal 
to the pressure of the undisturbed uniform stream. 


Inside the jet we have 


; Od _ Op 
u=l (1 } =) y = a. (2) 


where Uq¢, is the velocity potential of the disturbance 
in the jet. 

The uniform stream is of a velocity U» and Mach 
Number J/, when it is not disturbed. 
components in the surrounding stream after it meets 


The velocity 


the jet will be 


of) . Ode 
= OP l " v= U2 é 
u ( $= > (3) 


where U2¢» is the velocity potential of the disturbance 
in the surrounding stream. 
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We also have 
1 > 04,/dx, 9¢;/dy 


where 7z = | or 2. 
The differential equation that governs ¢; is 


0*¢; | 0°; 
ae p 


d— M? 
Ox Oy" 


= 0 (4) 
The boundary conditions of the problem are as fol. 
lows: 
(1) By symmetry, at the axis of the jet, the y-wise 
velocity component is always zero—i.e., 


Od: Oy = 0 (5) 


at y = Ofor all x. 

(2) On the surface of the jet y = +1, we have two 
boundary conditions: 

(a) The pressures are the same on both sides of the 


interface—i.e., 
0¢)/Ox = (M2?/M,") (0d2/0x) (6) 


aty = +1. 

(b) The inclinations of the flow are the same on 
both sides of the interface—.e., 

0¢1/Ov = O¢2/Oy (7) 

aty = +1. 

(3) The boundary conditions at infinity depend on 
whether the uniform stream is supersonic or subsonic. 

(a) Ifthe uniform stream is subsonic—t.e., MJ. < 1 
we have that the disturbance will die out at infinity 


1.€., 
Og2/Ox = Ode/Oy = 0 (Sa) 


at y = +o, 

(b) If the uniform stream is supersonic—i.e., M2 > 
1—-we assume that no disturbances come from infinity 
and that all disturbances are generated from the inter- 


face—i.e., 
g2 = F(x - Bo|y ) (Sb) 


where F is an arbitrary function that is going to be de- 
termined and where 8.” = JA/.” — 1. 

Since the nature of the subsonic and of the super- 
sonic flow of the surrounding stream differ entirely 
from each other, we have to treat the problem sepa- 
rately for these two cases. 


(III) Supersonic JET IN SUBSONIC UNIFORM STREAM 


Since the differential Eq. (4) is linear, the method of 
superposition is applicable. It is then sufficient for us 
to consider one typical fundamental solution only in the 
analysis. In order to satisfy the boundary conditions 
[Eqs. (5) and (8)], we have the fundamental solution 


of ¢; as follows: 
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SUPERSONIC FLOW OF A 


¢; = cos \fyy (Ai sin Ax + B, cos Ax) t 9) 
od = e *'¥\ (4. sin Ax + Be cosAx) § Y) 


where A; and B; are arbitrary constants, \ is the eigen- 
value of the problem, and 6,7 = M,? — 1 and £6,? = 
1 — M.’. 

In order to satisfy the boundary conditions, Eqs. (6) 
and (7), we obtain the eigenvalue equation as follows: 


M? |l — M.? 
tan Ap; = 2 - = tan@ (say) (10) 
M,? M,? — | 
where 0 < 60 < 7/2. 
The eigenvalue \ is then 
A = (6 + nm)/B; (11) 


where 2 = 0, 1,2, 3,.... 


The complete solution of the problem is 


: . 6+ nr 
>> cos (6 + um)y EE sin r x + 


n=0 1 


6+ n 
Bi, cos 4 x] (12) 
B 


. O+ nr 
eee | 4s, sin x + 


n=Q 1 


06+ ur 
Bs, Cos x] (13) 


Bi 


The coefficients A,, and B,, can be determined from the 
initial conditions Eq. (1). After A;, and B,, are 
known, A», and By, can be determined by either Eq. 
(6) or (8). 

It is interesting to know that ¢; and @¢» of Eqs. (12) 
and (13), respectively, are, in general, not periodic 
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functions but that they are almost periodic functions 
which were first introduced by Bohr.’ Hence, the flow 
pattern in the supersonic jet in uniform subsonic 
stream is only almost periodic. As a first approxima- 
tion, the wave length of this almost periodic structure is 


L= (23 @) By (14) 
The variations of 6 with 7; and JJ. are shown in Figs. 
It is seen that this wave length increases as 
As Ms = 0, L = 1;, 
for two-dimensional 


2 and 3. 
Moe increases for a given MM). 
Prandtl’s formula 


which is the 


supersonic jet. 
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so that the boundary conditions [Eqs. (5) and (8)] are 
satisfied. In order to satisfy the boundary conditions 
[Eqs. (6) and (7)], we have the eigenvalue equation as 
follows: 


Ee 82? . | , 
- + tan? A6,] [1 + cot? \B.] = 0 (16) 
M,! B;? 

It is evident that there is no real eigenvalue \ that 
satisfies Eq. (16). Hence, the supersonic jet will not 
have periodic or almost periodic structure. 

The solution of this case may be obtained by method 
of characteristics. It is also possible to find some ana. 
lytic solution for simple cases, as follows: 

The most general solution of this case is 


Fic. 6. Supersonic jet in uniform stream (MM, = 1.98: Mz = —— . 
as ss e id pie = f(x — 1) a(x y) » 
1.40). (Fig. 14 of reference 7; courtesy of J. G. Wilder, Jr., $1 It Bry + g(x + By} (17) 
Cornell Aeronautical Laboratory, Inc.) go: = Fix — fr y|) f 
where /, g, and F are arbitrary functions to be deter- 
(IV) SUPERSONIC JET IN SUPERSONIC UNIFORM mined from the initial conditions. In order to bring 
STREAM out the essential feature of the flow pattern, we consider 


First, we would like to know whether or not the super- 


the simple case that, at x = 0, 


sonic jet still has periodic or almost periodic structure. u(y) = AU, wl(y) = 0 (18) 
We put ; 
I for —1 < y < 1, where A isa small constant. 
¢: = cos AP [A; sin Ax + B, cos X x] (15) In order to satisfy the boundary conditions on the 
° | F o rs é “i 
dé: = Az sin A(x — B2\y|) + Bz. cos A(x — B2\y\)f interface y = + 1, we have 


di = fo(x — Biy) + file — By) +... +Shile — By) +... + gol(x + By) + 


where 


and 


where 


gilx + Bry) +... + gale + Bry) +... (19) 


A , 
n(x — Bry) = R"-— (x — Bry) for (2n — 1)Bi <x — By < (2n+ 1)B, 


> (2n + 1) B, 


= 0 f se sia dia 
; a : Px < (2m — 1)B; 
(20) 
A . 1 
g(x + By) = oie (x + By) for (2n — 1)Bi <x + By < (2n+1)B, 
. > (2n + 1)B; 
=(Q x , 
( forx + By < (2n — 1)6, 
»= F(x — Bs| y|) + Fi(x — Boly|) +... + Fia(x — Bolyl) +... (21) 
a A : ' 
= FR nan (x — Be|y|) for 2nBi — Bo<x — Boly| < 2(n + 1)Bi — Bo 
~ (22) 
: 1 >2(n+ 1)B; — Be | 
=> 0 f == 91V 
tics Bs , 2 2nB — Be 
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SUPERSONIC FLOW OF 
T is the transmission factor that is the ratio of the 
nondimensional magnitude of disturbance in the uni- 


form stream to that in the jet and which is given by the 
following formula: 


‘iat M,? M,? rans l 
T=2/(— +4/— 
M,? M,? - 1 
The values of 7 for various 14, and Ms, are shown in 
Fig. 4. 
R is the reflection factor that is the ratio of the mag- 
nitude of the reflected disturbance from the surface of 


(23) 


discontinuity to the magnitude of the incident disturb- 
ance and which is given by the following formula: 


R Ga M,* — , (eS ee ) 
| he M2 — 1// \M;? M,? — 1 


(24) 
Its value lies between —1 to +1. 
various \/, and MJ, are shown in Fig. 5. 


The values of R for 
It is seen that 
the reflected disturbance may be of the same sign, or of 
opposite sign of the incident disturbance, or zero ac- 
cording to the values of AJ, and Mp. 

There are few experimental data to check the present 


theory. However, the schlieren pictures (Fig. 6) taken 
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by Wiider’ indicate that there is no periodic structure 
in the supersonic jet in uniform supersonic stream. 
Since the difference between 1/4, and Mz is small in 
Wi'der’s pictures, R is extremely small; as a result, in 
most of his pictures, only the first transmitted dis- 
turbance F\(x — 2|y| ) is seen without reflection. In 
the reflected 
Si(x — Bry) and g;(x + B,y) are seen. 


only a few cases, first disturbances 
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Buckling of an N-Section Column 


G. Sri Ram and G. V. R. Rao 

Research Student and Assistant Professor, Respectively, Indian 
Institute of Science, Bangalore, India 

September 29, 1951 


yen A COLUMN BUILT IN at one end and free at the other, 
subject to an axial load ?. The failure of the column is 
assumed to take place due to lateral instability rather than by 
direct compression. When the column is made of n-sections along 
its length, an analytical solution for the critical load P,;. can be 
obtained as shown in Fig. 1. 





In tn 








Let the 7-th section of the column have a length /; = L; Li-t 
and moment of inertia J;. This section can be characterized by 
kili, where kj? = P/EI;. The governing differential equation for 
the deflection of the column in 


EI, (d?y;/dx?) = P(6 Vi) (1 


Thus for the n-sections along the length we have » such differen 
tial equations. One can assume 


66 


y = 6[1 — cos (xx/2L)], one obtains a critical load of P.,. = 


Mi = 6 + Aji cos ki x + B; sin kx 9 


as a solution of Eq. (1). From the continuity of the deflectioy 
curve of the column we can equate the deflection and the slope of 
the column at the i-th joint obtained from the i — 1 section anq 
i-th section. These together with the conditions at the ends of 
the column give equations to determine the constant coefficients, 
A; and B;. If these coefficients should have nonzero solutions, 
one must satisfy a determinantal equation. The expansion 
of the determinant leads to a transcendental equation, the soly. 
tion of which gives the critical load of the column. 

For example, in the case of a column built up in two sections 
the equation 


O = 1 — (K,/Ke) tan &l; tan Role 


gives the buckling load of the column. This condition is the sany 
as that given in references 1 and 2, whereas, for » = 3, there 


seems to be an error in reference 2. For n = 3, we must satisfy 


O = 1 — (K,/Ke) tan &/; tan kele 
(K,/K3) tan Rl; tan kl (K2/K;) tan kel, tan hl; 


For n = 4, we must satisfy 


O = 1 — 3(Ky/K,) tan kglg tan k,l, 4+ 


(K,/K2)(K3/K4) tan Rl; tan Role tan kyl; tan ky, 


where, in the summation, ry = 1, 2, 3, 4and gq < r. 
f 


In general, for n-sections we obtain an equation of the form 
O = 1 — DXXfi/K;) tan &jl; ie = 2 3... . 2) 3) 


Each f; in the above equation depends upon the functions /,, f,, 


fi = fi + ki, tan Ri_yli_y (1 — Xd(f-/k,-) tan Rl { 


summed for r = 2,3,...,2 2,and f, = K;, tan &,/ 

With the above formula one can write the characteristic equa 
tion for any number of sections in the length of the column. Th 
equation that gives the critical load can only be solved by trial 
and error. As ” increases, the work becomes laborious. In prac- 
tice, one can see that tan &;/; is small and that products of tan 


k,l; four times or more can be neglected. Eq. (3) will then reduce 


to 
0 = 1 — 3(K,/K,r) tan kgl, tan yl; 5 
where r = 1,2,...andq<~¥, 


One can easily show that the f; in Eq. (3) are all positive 
Further, Eq. (4) shows that the approximation made by neglect- 
ing higher powers of products of tan &;/; gives a higher value for 
f;. Considering the plot of the right-hand sides of Eqs. (3) and 
(5) against P, one can easily prove that the solution of the ap 
proximate Eq. (5) is always lower than the solution of the com 
plete Eq. (3) 

From the application of the energy method, one obtains 
approximate buckling load, always higher than the true buckling 
load Therefore, the method outlined above helps to obtain 
bounds on the buckling load of a column. 

A symmetrical column with both ends pinned can be reduced to 
the case treated above by considering only half the length of the 
column. For an unsymmetrical column pinned at both ends 
similar formulas can easily be obtained. A numerical example 1s 
carried out applying this method for a column built up in fout 
sections. The following data are taken for the example 


I, is the section of the column at the built-in end and 


I 
= 0.81, = ().64, = 0.36 


L = length of the column. The approximate solution obtained 
from Eq. (5) gives a critical load of P,, = 1.979 (EI,\/L? 
From the strain energy method, by assuming a deflection curv 
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9034 (EI,/L? The analytical solution as given from Eq. (3) 
2.019(EI,/L?) 


yields P 
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A Solution of the Surging Problem in Axial- 
Flow Compressors 


John Mitchell Stephenson * 
Farnborough, England 
August 21, 1951 

SUMMARY 


f an axial-flow compressor has many stages on a single rotor shaft, and 


there is only one combination of pressure ratio and mass flow 


adiabatic flow 
This is called the 


for which the efficiency of compression is a maximum 
“design point’ on the performance graph and usually corresponds to full load 
At any other point the efficiency is reduced, and surging 


or thereabouts 
This effect places an upper limit of 


may occur with catastrophic results 
about six on the design pressure ratio of such a compressor 
The reason for this limitation is analyzed, and the known ways of avoiding 


it areenumerated. A further solution is then described, which depends on ex 


traction of heat from the gas during compression It is shown that high 
pressure ratios can be achieved with good efficiency and stability and without 
mechanical complications 

(1) INTRODUCTION 


-— FLOW THROUGHOUT AN AXIAL-FLOW COMPRESSOR Can most 
conveniently be deduced in terms of the performance of each 
stage, by which is meant each pair of annular blade rows. To 
compress efficiently, none of the blades of the compressor must be 
stalled. But the chief limit to the operation of axial compressor 
is surging, which occurs when most of the stages are stalled. The 
mechanism of surging is not fully understood, but the effect of 
violent fluctuations of pressure is well known. Surging must al- 
ways be prevented, and the most important characteristics of a 
compressor is the surge line, which describes the maximum steady 
pressure ratio attainable with a given mass flow. 

We shall show that the line of most efficient running of a single 
shaft adiabatic compressor is close to the surge line. This fact 
limits the design pressure ratio of such a compressor to about six. 
Above this it is not usually possible to avoid surging at part 
loads. 

(2) THe EFFICIENCY OF AN ELEMENTARY COMPRESSION 


Consider a gas flowing steadily along a uniform tube, such that 
at one section there is an infinitesimal pressure rise dp. If the 
work put into the gas to effect this compression is dW per unit 
mass flow, we may define the efficiency » of the process as the 
tatio of pressure energy produced to the work applied. That is, 

7 = dp pdW 
where p is the density of the gas 
If the compression is adiabatic dW is equal to the change of 


enthalpy, Cpbd7. But if an amount of heat dQ per unit mass flow 
isextracted from the gas during the compression, we have 


dW = CpdT + dQ 
We now define a cooling factor ¢ as the ratio of heat extracted to 
the work done, or 
¢ = dQ/dW 
Combining these equations we have 
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Or 








0 ds , $e 7 : , ; 
06 07 08 09 4 1-0 


Fic. 1. 


n/(1 — §) = (dp/p)/(Cp dT) 


The differential law relating pressure and temperature is ap 


proximately 
dp/p = [n/(n 1)}(d7/7 


where » = y for an isentropic compression (i.e., 7 = 1 = ¢). 
Hence, using the perfect gas equation of state, 


y - 1 
. ) pCpl 
Y 


n (? l ( n ) 
oa (1) 
i-—¢ ¥ n l 


In the particular numerical case that we shall consider later, 


p = pRT = 


we have 


y = 7/5and n = 6/7, so that 
n = 3/(2 + ¢ (la) 
Note that, if ¢ = 0, the compression is adiabatic; while if 
¢ = 1, we have m = 1 and the compression is isothermal. We 


shall assume here that ¢ always lies between these limits so that 
no heating or refrigeration is applied. Clearly, for a given value 
of the efficiency (which is independent of the cooling), we can re 
duce n from the adiabatic value to any other value greater than 


one simply by cooling the gas 
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(3) THE FLow RANGE OF A SINGLE COMPRESSOR STAGE 


The dimensional analysis of a compressor stage shows that 
there are a small number of coefficients in terms of which the per- 
formance can be described. For a stage of given geometrical de- 


sign the only quantities that may be varied independently are 


va = the mean velocity of the gas in the axial direction 

u = the tangential velocity of the rotor blades at a reference 
radius 

l = a length giving the geonietrical scale of the stage 

v = the kinematic viscosity of the gas at same point 

a = the velocity of sound in the gas at inlet 


The variables can be combined to form three dimensionless 
ratios, each a measure of one aspect of the physical behavior of 
the gas. These are the Reynolds Number, Re = val/v; the 
Mach Number, WV = va/a; and the flow coefficient ¢ = va/u. 

The dependent quantities of interest are the pressure rise per 
unit mass, Ap/p, and the efficiency. The pressure rise of a single 
stage is as a rule sufficiently small to use the elementary efficiency 
» as defined in Section (2). The pressure rise can be expressed in 
the dimensionless form 


y = Ap/(1/2)pu? 


We may say, therefore, that the performance of a stage is given 
by ¥ and », which are unique but as yet unknown functions of the 
independent parameters ¢, Re, and M. The effect of the last two 
is known to be small, provided that Re > 3 X 105 and M < 0.5, 
approximately. Thus the performance can be expressed simply 
in terms of ¢, which determines the angle of attack of the gas on 
the blades. 

We shall not now derive y and n analytically in terms of ¢ and 
the blade shape but shall simply show the characteristics of a 
This stage is designed to operate at the 
point of best efficiency, where ¢ = ¢d. We see that there is only a 
small range of ¢ over which the stage works efficiently. And at 
values less than ¢s, Y has a negative slope, and the blades are 
If the stage is running alone, the flow surges at these 


typical stage in Fig. 1. 


stalled. 
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points. 
surging does not occur until most of the stages are stalled, 


(4) THE COMBINATION OF SEVERAL STAGES IN SERIES 

The problem of combining stages so as to produce a high pres- 
sure ratio is to ensure that all the stages operate at or near the 
design value of ¢ under all likely conditions of operation. 

Consider a compressor with many stages, each having the per. 
formance shown in Fig. 1. We take any two separated stages jn 
this compressor and refer to them with suffixes 1 and 2 The 
mass flow Q is given for any section of area A by 

Q = pAVa 


The law of compression is 


p2/p: = (po pi)! 7 = r! n 


2 = (2a) ( l ) 
91 - QA ote ils 7" 


A, is determined for a particular density 


Hence, 


Now the area ratio A, 
ratio and is fixed thereafter. Let us suppose that this corresponds 
to a design pressure ratio R, when there is no cooling (¢ = 0) and 
the efficiency is 7 = 6/7. Then, from Eq. (la), we have n = 3/2 

If also at the design point “; = w2/k where k is a constant and 
Qi = Qe, then 

A,/A. = RR*/* 

and Eq. (2) becomes 


HD) ‘bt = (Qeuy Qyu2) (RR® " r} *) (3 


A typical simple compressor has continuous adiabatic flow and 
Thus u, = w2/k, Q: = Qe, and m = constant 
It is clear from Eq. (3) that in this 


a single shaft. 
(approximately) for all r. 
case 
¢d2 > ¢ forallr < R 

Thus, the stages are never all working at maximum efficiency, 
except at the design pressure ratio and mass flow. The (r, Q, u 
diagram of such a compressor is shown in Fig. 2 with the surge 
line. 

Now ¢1 = Q:/:Aim, and therefore the condition that ¢; is con 
stant implies that Q,/m is constant. Similarly, ¢2 is constant im- 
plies that Q:/mr!/” is constant. The lines of constant ¢ are shown 
in Fig. 2 for the case where section 1 is at inlet and section 2 at 
outlet. Since as we have seen, ¢2 > ¢, neither of these lines corre- 
sponds to efficient compression, and in the second case the com- 
pressor surges when r < 0.9. The line of best efficiency actually 
lies somewhere between these two and corresponds to constant ¢ 
for a stage midway between inlet and outlet. It is parallel to the 
surge line, which approximately corresponds to stalling of the 
midway stage. 

To ensure stable running, especially when R is large, it is 
therefore necessary to design every stage, and in particular the 
first few, to have a large flow range. We shall not consider this 
subject here, but there are more general solutions of this problem, 


some of which are well known. 


(5) SOLUTIONS 


(a) The most simple and specialized solution is to keep the 
pressure ratio always at the design value R. When a compressor 
forms part of a closed system with a continuously circulating gas, 
it is sometimes possible to alter the inlet density level at will. In 
this way the mass flow and power may be altered over the whole 
range, while the rotational frequency remains constant. Thus 
the efficiency is always a maximum. This solution is most suit 
able in connection with constant speed electrical apparatus 

(b) A device that has long been used on airplane propellers is 


to rotate the blades about a radial axis. In this way the flow co- 


However, when several stages are combined in Series 
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TABLE 1 


efficient ¢¢ corresponding to maximum efficiency can be varied 
over a large range. Of course the twist of the blades cannot be 
altered, so that the efficiency falls off slightly when they are ro- 
tated. It is difficult to adjust the rotor blades, but, if only the 
stators are rotated, the efficiency can be improved and surging 
avoided. For example, when r < R, the entry stator blades and 
the stators of the first few stages should be rotated to a high 
stagger (fine pitch), while the stators of the last few stages are 
moved to low stagger (coarse pitch). 

From Eq. (3) however, we see that ¢; = ¢2 can be maintained 
for all pressure ratios less than R by reducing either (c) Q2/Q,, (d) 
u;/u2, (e) n. We discuss these in turn. 

(c) In a compressor with inlet at atmospheric pressure it is 
easy to reduce Q2/Q, below the usual value of one by blowing off 
gas between the stages. However, a great deal of power is lost in 
doing this, so that the efficiency of any plant incorporating such a 
device is low. Blowoff is only defensible therefore for transient 
conditions in starting or decelerating compressors. 

(d) “,/u2: can be appropriately reduced below the design value 
k if all the rotors are able to turn independently. This arrange- 
ment was suggested by Dr. A. A. Griffith in 1929. It is compli- 
cated but is frequently approached by dividing the stages into two 
groups on separate shafts. 

(e) The solution put forward here and which is in theory most 
promising for many applications is to reduce the exponent n by 
cooling the gas during the compression, such that m is given by 
Eq. (1). A device to effect this was patented by the author and 
Power Jets Ltd. in July, 1948. 


(6) CooLED COMPRESSION 


We consider a single shaft compressor without blowoff, so that 
Eq. (3) becomes 

¢:/¢1 = R’*/r 

and we wish to adjust so that ¢@2 = ¢, for all r. Hence, we re 


quire 
log r/log R = 2n/3 
Using Eq. (la) we have 
log r/log R = 2/(2 + $) 


If fis varied according to this formula, all the stages work at their 
design efficiency. But since we have assumed that ¢ < 1, there is 
a value 7) below which it is not possible to reduce m when n = 1 
and the flow is isothermal. Then 


log ro/log R = 2/3 

Values of ro and R are compared in Table 1. We see that, if a 
single-shaft isothermal (or adiabatic) compressor will run stably 
up to a design pressure ratio of five, then a cooled compressor 
having the cooling shut off at full load can be designed for a pres 
sure ratio of eleven, with the same stability. 

In practice it may not be desirable to vary ¢ and » continu- 
If, for example, full cooling is applied up to one-half the 
required pressure ratio and switched off completely above that 


ously. 


the compressor may be sufficiently stable. 

Suppose the cooling is achieved by spraying liquid into the 
stream, so that the liquid is at the temperature of the stream and 
does not evaporate. Then, if the specific heat of the liquid is C 
and the liquid-gas ratio m by weight, we have 


dQ = mC dT 


13 14 15 
6.09 


9 10 11 12 
4.35 4.66 4.96 5.24 5.53 5.79 


Hence 
¢ = mC/(Cp + mC) 


Note that, for isothermal compression, m — ©. For air and 


water, Cp = 0.24C and 
¢ = m/(0.24 + m) 

An approach to the solution given here is made by varying the 
thermal ratio of an intercooler placed between two multistage 
compressors that operate in series. This easy way to 
achieve freedom from surging troubles on am intercooled gas- 
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turbine engine. 


Plastic Buckling of a Simply Supported Plate 
in Compression 


Richard A. Pride 
Aeronautical Research Scientist, Langley Aeronautical Laboratory, 
Langley Air Force Base, Va. 


September 25, 1951 


ges APPLICATION OF THEORIES OF PLASTICITY to plate buckling 
problems has resulted in widespread controversy in recent 
years. In the case of long simply supported flanges, a large 
discrepancy exists between the predictions of buckling theories 
based on deformation or finite theories of plasticity’ ? and those 
based on flow or incremental theories of plasticity. The experi- 
mental data are in excellent agreement with deformation-type 
buckling theories and appear to contradict altogether the flow 
type of buckling theory. It has been both suggested‘ and denied® 
that the discrepancy between the experimental results and the 
predictions of flow-type buckling theories can be ascribed to the 
effects of the initial imperfections. A crucial question in such a 
debate is how large the initial imperfections actually were. 

In the case of simply supported plates the discrepancy between 
the two types of theory is generally much less pronounced. Here 
also experimental results® are in much better agreement with de- 
formation than flow-type buckling theories. Several investiga- 
tors have questioned the validity of this comparison, however, 
because two of the deformation type of buckling theories used in 
the comparison were based on the Shanley principle of buckling 
with increase in load, whereas the flow type of buckling 
theory was based on the assumption of constant loading during 
buckling. Pearson’ has recalculated the buckling load predicted 
by flow-type theories for simply supported plates taking the 
Shanley principle into account and found that this modification 
substantially improved the agreement with experiment in the 
lower stress range. Although Bijlaard* has pointed out that a 
large part of the improvement is due to the use of the simplifying 
assumption that Poisson’s ratio is one-half in the elastic, as well 
as the plastic, range, it seemed desirable to provide an experi- 
mental check in a region in which the discrepancy between the 
two types of theories is much more pronounced. 

The discrepancy between the theories becomes greater at high 
buckling stresses—that is, for plates having large thickness- 
width ratios. However, it would not be satisfactory to go far in 
this direction because the thin plate equations on which all of 
the plastic buckling theories have been based would cease to be 


applicable. A better way to provide the desired check is to use 
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Fic. 1. Plastic buckling theories compared with tests. 


a material with a low yield stress. With this purpose in mind, 
square tubes of 145-T6 aluminum alloy similar to those with the 
highest buckling stress in a previous experimental program® were 
annealed before being buckled, producing the results shown in 
Fig. 1. 
per sq.in., and the buckling stress was reduced from an average 
(three tests each) of 60,400 to 19,700 Ibs. per sq.in. as a result 


The yield stress was reduced from 62,400 to 12,300 Ibs. 


of annealing the material. 

Deformation-type plastic buckling theories continue to be in 
good agreement with these tests. Flow-type theories, however, 
are completely inadequate even with Pearson’s modification. 
The dashed portions of the flow theory curves in Fig. 1 indicate 
the trend of the theory if the stress-strain curve is assumed to 
rise to a correspondingly high stress. On this basis Pearson's 
analysis gives the same buckling stress for the low-strength an 
nealed material as for the high-strength material. Since for 
reasonable strains the stress-strain curve of the material does not 
rise indefinitely, the flow theory curves probably should be cut 
off in the vicinity of the ends of the solid lines, Fig. 1. This 
would seem to imply that buckling could not occur at all for the 
annealed tubes tested 

The small scatter of the test results makes it seem doubtful 
that the large discrepancy between these results and the predic 
tions of flow-type buckling theories can be ascribed to the effects 
of initial imperfections. However, this evidence alone cannot be 
considered conclusive. Additional data bearing on this possible 
explanation are provided by the representative stress-deflection 
curves shown in Fig. 2, from which the order of magnitude of the 
initial imperfections can be inferred. 

The results just discussed taken in conjunction with the pre- 
viously established results for buckling of long simply supported 
flanges confirm the view that for practical purposes the buckling 
stresses of well-made plates are in good agreement with the pre 
dictions of buckling theories based on deformation theories of 


plasticity without modification for effects of initial eccentricity 
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Fic. 2. Stress-deflection diagrams from buckling tests 
(Plate thickness, ¢ = 0.154 in.) 
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A Method for Calculating Airfoils with 
Prescribed Pressure Distribution 
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ijn PROBLEM OF CALCULATING AIRFOILS with prescribed pres- 


sure distribution is important for laminar airfoils and suction 
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Fic. 1. Source distribution and airfoil shape in the x, y-plane. 
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Fic. 2. Diagram for evaluating the coefficients a; and a» 


airfoils. Furthermore, the problem of cavitation in fluid flow 
and the problem of critical Mach Number at subsonic flow are re 
A solution by means of the method of singulari 
and Allen.’ 


extension of the method by Riegels suited for practical 


lated herewith. 
ties has been given by Riegels' In this note we give an 
applica 
tion. 

he starting point of our calculations is the relation given by 


Riegels* between the velocity of the potential flow on the surface 


of the airfoil, U., and on its mean line, Um, see Fig. 1, 

Us(x) = Um(x)/V 1 + (dy/dx)? (1) 
y(x) being the shape of the profile. Between the source distribu 
tion g(x), which induces velocities on the mean line, and the air 


foil shape the following relation holds: 


dy/dx = (1/2) [q(x)/U.] (2 
Therefore one gets from Eq. (1) for the velocity at the surface of 
the airfoil 
\ l cow ae. alae } 
l 4 
Ux) = I w J x’ 4 3 
V1 + (dy/dx)? 
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For the airfoil shape y(x), at first supposed to be symmetrical, 
we give an expression, which fulfills identically the condition of a 
0 for x = +c/2), 
trigonometric variable, 


closed contour (y = ¢ being the well-known 


x/(c/2) = —cos¢ 
y/(c/2) = y dy sin (ng) 
n=1 (4) 
dy cos (Ng) 
= > Nan - 
dx as sin ¢ 


with 0 < ¢ < a for —c 
then becomes 


2<x< +c/2 
(compare also reference 2) 


The integral in Eq. ( 


+c/2 , "x 
1 dy dx 1 
. Ix'x — x’ a x 
TJ xt me me fQak x x T J olnd 
. _ cos (ng’) sin ¢’ 
ps» n(Qn) . : > dy’ 
nant sin ¢ cos ¢ — cos¢ 


> sin (ng) - 
iu N(dy) . (oO) 
oe sin ¢ 
Herewith one obtains finally from Eq. (3) after some simple 
calculation 


} n(a,) sin (ng) = 


n=1 


sin g X 


en 


cos (ng) \? 
ea! + > x(a, ) : - \ fie) (6) 
} 


sm} sin ¢ 
For slender airfoils the square root at the right-hand side of this 
equation is only a little different from 1 for nearly the whole range 
taken as a correction term 
With 


can be regarded as a 


of values of ¢. Therefore, it can be 
that can be approximated in a simple way as given below 


Us(¢ 


known function f(¢), 


) given, the right-hand side of Eq. (6) 
from which the coefficients a, can be calcu 
lated by Fourier analysis 


9 " 
an = fig) sin (ng) deg (7) 
mn) Joao 


Then the shape of the airfoil is obtained from Eq. (4 By the 


aid of the Prandtl-Glauert rule, the effect of compressibility can 


simply be introduced by putting 


f(o) 


¥)incomrpr 


a V1 M 2 f(¢) compre 3) 


M... being the free-stream Mach Number 


For the first approximation of f(¢), only the coefficients a, and 


ad, are taken into account. These can simply be obtained from 
Fig. 2 dependent on the following two parameters of the pre 
scribed velocity distribution: 

(1) Maximum velocity of the potential flow, | l 

(2) Position of maximum velocity, [«/(c/2)}; war, correspond 


ing tO Yma 
This diagram 
Joukowsky 


tained from Eq 


has been obtained from the “generalized” 


airfoils (with rounded trailing edge), which are ob 
e zero, except for a, and 
Joukowsky 


approximation for f(¢) ob- 


1), if all coefficients ar 


ad». (The case dz = a;/2 gives the “‘special”’ airfoils 


The 
tained in this way yields the coefficients a, 


with sharp trailing edge first 


according to Eq. (7) 


and herewith a good approximation for the shape of the airfoil 


From these coefficients, a», a second approximation for f(¢) can 


be obtained, which in most cases then yields the shape of the air 


foil with fully sufficient accuracy. It is sufficient to calculate 


in the first approximation the coefficients until » = 4, and in the 


second approximation until » = 8. For the evaluation of 


the coefficients a, and the coordinates of the airfoil, a simple 


general scheme has been set up which reduces the whole numerical 
Fig 


work to forming simple products and their sums 3 gives the 
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1. approximation 
—-— 2.approximation 
3. approximation 
——— Method of Allen [NACA Rep. 833) 
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Fic. 3. Prescribed velocity distribution and calculated airfoil 
shapes. (Airfoil thickness is enlarged to five times the actual 
thickness. ) 


results of the first, second, and third approximation of our calcu- 


lation compared with the result of Allen? for one of his examples 


(base airfoil velocity for ‘‘semi-low-drag-airfoil’’). Up to now, 
only symmetrical airfoils with the same velocity distribution on 
both sides of the airfoil were dealt with. But it can easily be 
shown that the same procedure can be applied if the prescribed 
velocity distribution is different on the lower and upper surface 
of the airfoil. The procedure then has to be applied separately 
for the lower and the upper surface and yields the airfoil coordi- 
nates for both sides. 

A pure source distribution g(x) yields the velocity distribution 
(Uup + Uiow)/2, a pure vorticity distribution y(x) yields the veloc- 
ity distribution (Uu, — Uiow). Both together above the chord 
yield a velocity normal to the chord 


2r 
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which might also be obtained from a pure source distribution gf 
the strength gup = 2vup. Then Eq. (2) holds for the slope of the 
upper surface with the required source distribution gup, and in the 
same way for the lower surface. Therefore, all relations given 
here hold also for asymmetrical airfoils. 

If, however, the airfoil has an angle of attack of some magnitude, 
Eq. (4) is inconvenient for the shape of the airfoil, because the 
Fourier series with only sine terms is incomplete, as it gives al. 
For small angles of incidence, as they 


ways y = Oatx = —c/2. 
mostly occur in practice, this failure is of minor importance, Fog 
larger angles of incidence, it can be overcome by splitting off that 
part of the induced velocity which gives the angle of incidenge 
From Glauert’s theory of thin airfoils and with regard to Eq. (1), 


one obtains for the angle of incidence, 


Viow | dy >| — v 
ra f sin 9 de (10) 


T JU * 
— +. 


where the braces mean the distribution of vorticity y of the 
The reduced velocity distribution having no angle of in 4 


a= 


airfoil. 
cidence is, then, 
. ; ay ctg (¢/2) 
Uw * = Ui ; (11) 
low V1 + (dy/dx)? 


In Eqs. (10) and (11) the slope dy/dx has to be taken from the 


symmetrical basic airfoil. 
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Turbulent Boundary Layer on a Cone in a Supersonic Flow at Zero 


Angle of Attack 
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